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DR. K. R. RAMANATHAN ON HIS SIXTIETH 
BIRTHDAY 


THE present issue of the “‘ Proceedings” of the Academy commemorates 
the sixtieth Birthday of Dr. Ramanathan which coincides with the date of 
its publication. Such commemoration appears appropriate in view of the 
eminent position he has attained as a man of science, the notable conttfi- 
butions which he has made to his chosen fields of knowledge and the leading 
part he has played in the development of an Indian School of Meteorology. 
Ramanathan is a Foundation-Fellow of the Academy and served conti- 
nuously for fifteen years as a member of its Council, besides being a Vice- 
President for the triennium 1943-46. His constant presence and active 
participation in the annual meetings of the Academy greatly promoted the 
usefulness of those meetings. Significant contributions by him and his 
collaborators have appeared in the “‘ Proceedings” of the Academy. 


Ramanathan began his scientific career as a Lecturer in Physics and 
Director of the Observatory at Trivandrum. He left this position towards 
the end of the year 1921 to join the group working in the research laboratory 
of the present writer at Calcutta, and collaborated in the studies of the 
molecular scattering of light there in progress. His work in this field proved 
highly fruitful and secured for him the distinction of being the first Doctor 
of Science of the Madras University. He did not long remain at Calcutta, 
but left to join a teaching appointment at Rangoon, where, however, he 
continued his investigations. In the year 1925, he resigned his appointment 
at Rangoon and accepted a call to the position of Meteorologist under the 
Government of India. At various periods in his service, he was Director of 
the Agra Observatory, Director of the Colaba and Alibag Observatories at 
Bombay, Director of the Solar Physics Observatory at Kodaikanal, Super- 
intending Meteorologist at Poona and Officer on Special Duty at Delhi. 
Having joined the service at a somewhat later age than his contemporaries, 
his retirement became due before his turn to be Director-General of the 
Department came up. He, however, officiated as Director-General on several 
occasions. Following his retirement from service in February 1948, he 
joined the Physical Research Laboratory at Ahmedabad then being organised. 
He is there now actively one his investigations on the physics of the 
upper atmosphere. 


The wide range of subjects Ramanathan handled during his official 
career in the Indian Meteorological Service is reflected in the variety of topics 
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dealt with by him in his scientific publications. A selected list of the papers 
issued during his period of service is appended to this short summary of his 
career. It will be seen that terrestrial magnetism, seismology, and the 
spectrum of the night sky, besides meteorological optics and acoustics were 
all of interest to him. His principal preoccupation, however, has been the 
study of the thermal structure and movements of the upper air and it is in 
this field that his contributions are outstanding. A paper by him on the 
distribution of ozone in the upper atmosphere is appropriately included 
in the present issue of the “* Proceedings”. Ramanathan’s recent election 
to the Presidentship of the International Meteorological Association is a 


gratifying recognition of the position he has won for himself by his contri- 
butions to knowledge. 


It remains for me as President of the Academy to express my grateful 
appreciation of the prompt and kindly response of the eminent meteoro- 
logists in several countries who were invited to contribute papers to this issue 
of the “ Proceedings” and thus to join in honouring Dr. K. R. Ramanathan 
on his sixtieth Birthday. 

CORDIAL messages have been received from the following who, 


on account of shortness of time, have not been able to send in papers for 
publication in this number: 


Dr. J. Bjerknes, Los Angeles; Professor Dr. Rudoff Geiger, Munchen; 
Professor G. M. B. Dobson, Oxford; Professor Dr. F. W. Paul Gotz, Arosa; 
Professor Dr. E. Regener, Ravensburg; Professor Joseph Kaplan, Los 


Angeles; Professor Dr. H. J. Bhabha, Bombay; Dr. S. Chandrasekhar, 
Wisconsin. 


C. V. RAMAN. 
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Dr. K. R. Ramanathan on his Sixtieth Birthday 
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Dr. K. R. Ramanathan on his Sixtieth Birthday 
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POLAR AND TROPICAL AUROR&: AND THE 
ISOAURORAL DIAGRAM 


By SYDNEY CHAPMAN 
(The Queen’s College, Oxford, and 
The Geophysical Institute, College, Alaska) 


Received December 2, 1952 


It is a pleasure to be among those who pay tribute to the long, distinguished 
and devoted scientific career of Dr. K. R. Ramanathan, President of the 
Association of Meteorology (part of the International Union of Geodesy 
and Geophysics), on the occasion of his sixtieth Birthday. The topic here 
discussed is closely related to the subject (Terrestrial Magnetism) chosen by 
him for his Presidential Address to the Physics Section of the Indian Science 
Congress in 1939—one of the many branches of geophysics to which he has 
given fruitful attention. 


INTRODUCTION 


The aurora polaris or polar light is observed commonly in high lati- 
tudes in both hemispheres, as the northern lights or aurora borealis, and 
the southern lights or aurora australis. Its panoramic beauty and changing 
form, colour and position much enhance the interest of the night sky in 
those regions, even to primitive human minds; and to those who know 
something of its origins and its widespread associations—with the distant 
sun, and with the geomagnetic field proceeding from deep within the earth 
—its delights are intellectual as well as visual. 


India is far from the usual regions of its appearance; but on very rare 
occasions it has been seen there, at least in the more northerly parts; and 
I am told that there are passages in the sacred literatures of India that are 
interpreted as referring to the aurora borealis. Though the aurora is 
rightly named polar, it is a phenomenon not wholly withheld from tropical 
gaze, either north or south of the equator. When visible from the tropics 
the aurora may justly be accorded the name tropical; in general, if not 

invariably, on the rare occasions of its occurrence the aurora tropicalis and 
the aurora polaris are one phenomenon. It is one to which, on account of 
its rarity in India, little attention has been given even by the atmospheric 
physicists of India; but part of the purpose of this paper is to propose that 
records of its observation from India should be searched for among the 
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dusty records of the past—from newspaper files and archives of events of 
former decades—because of their value and interest to-day. 


2. THE IsOcHASMIC DIAGRAM (Fritz) 


Herman Fritz, Professor at the Eidgenéssische Polytechnic at Ziirich, 
was one of the founders of the science of auroral morphology—the study 
of the geographical distribution of auroral appearance over the earth, 
and of its changes in the course of the day, the year and the sunspot cycle. 
In 1873 the Vienna Academy of Sciences (Wissenschaften) published! his 
great catalogue (Verzeichniss) of aurora borealis (and australis), the result 
of immense labours of compilation of auroral reports derived from very 
widespread sources, indicated therein. Already in 1866 he had found from 
his studies* that the frequeacy of visibility of aurore does not increase 
steadily towards higher latitudes, but that there is a line or zone of maximum 
frequency of visibility, along which the aurora can be seen on every clear 
dark night; this is called the auroral zone. It is a nearly circular oval of 
radius about 23° angular or over 1,500 miles linear, centred not at the geo- 
graphical pole, nor at the north pole of magnetic dip (76° N., 102° W.), but 
at a point that can be identified as the pole of the geomagnetic axis of the 
earth (79° N., 70° W.); this point may therefore also be called the auroral 
centre. The frequency of visibility of aurora decreases on proceeding either 
north or south from this zone. As Fritz remarked on p. 67 of his treatise 
‘** Das Polarlicht ’’ mentioned below, the existence and rough position of the 
auroral zone had already been realized many years earlier by Muncke, who 
wrote an article “* Nordlicht” in the second edition of J. S. T. Gehler’s 
Physikalisches Worterbuch, 1825-1845; and in 1860, E. Loomis, in 
** Silliman’s ’’ American Journal of Science and Arts, 30, 79-100, clearly 
indicated that the frequency of auroral visibility has a zone of maximum, 
and published a tentative diagram showing a maximum zone in which the 
average yearly number of visible auroral nights is 80, and inner and outer 
bordering zones in which the number is 40; Loomis’s position for the 
auroral centre was about 83° N., 60° E. or about 5° distant from the more 
accurate position found by Fritz. Added 14 Jan. 1953. 


In 1874° he gave a map showing isochasms or lines of equal frequency 
of (nights of) auroral visibility; these oval lines lay outside the auroral zone, 
and were approximately parallel to it; they indicated the places along which 
aurore could be seen 100, 30, 10, 5, 1 and 0-1 times per year (in the average 
over many years). Inside the auroral zone he drew only one line (not an 
isochasm), along which the aurora is visible as often to the north as to the 
south—the line of ‘ neutral’ directivity of auroral visibility (Fig. 1). 
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The distribution of isochasms, or lines of equal auroral frequency, in the northern 
hemisphere, according to Fritz. 


In 1881 he published his famous treatise* “ Das Polarlicht ’, describing 
these and allied researches by himself and others, and summarizing the 
auroral knowledge and speculations of that time. 


3. MODERN ISOCHASMIC DIAGRAMS (Vestine, Snyder) 


Though Fritz’s diagram was constructed by a simple method lacking 
in some refinements, it was a most valuable contribution to geophysics; and 
it was not superseded until, after 80 years, Vestine® produced a modern 
isochasmic diagram, constructed by improved methods, and based not only 
on the data collected by Fritz—extending from 503 B.C. to 1872 August— 
but also on the new data furnished by the two International Polar Years 
(1882/3 and 1932/3), and by many other modern sources. The auroral zone 
and the six outer isochasms given by Fritz needed some but not much amend- 
ment; and Vestine not only gave 9 instead of 6 outer isochasms, but also 
4 inner isochasms; the innermost has a radius rather less than half that of 


\ 
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theauroral zone, and corresponds to 50 nights of auroral visibility per 
average year. 


The outermost isochasm given by Vestine, as by Fritz, was still for 0-1 


visible occurrence per average year. In the longitudes of India it lies nearly 
along the circle of latitude 47°, far to the north of India. 


A year later, Vestine and Snyder,® improving on pioneer efforts by 
Davies,’ and by Geddes and White,® published a valuable isochasmic map 
for the aurora australis. Besides showing the southern auroral zone, two 
inner and three outer isochasms were drawn; the outermost corresponds to 
5 visible occurrences per average year; the material on which this diagram 
is based is far less than that which underlies the northern diagram. 


4. NEED FOR ISOAURORAL DIAGRAMS 


Valuable as are these important diagrams, they need to be supplemented 
by a new kind of auroral diagram giving isoauroral lines or isoaurores. These 
are defined as lines at all points along which there is the same frequency of 
(nights of) auroral occurrence above the point, instead of the same frequency 
of visibility from the point. Though to the general public such diagrams 
might not be of such direct interest as the isochasmic diagrams, they would 
represent more precisely the geographical distribution of occurrence, and 
would be of greater value to the auroral theorist, by indicating the extreme 
limits of latitude of occurrence which auroral theory, however satisfactory 


it may become for the ordinary high latitude aurore, must also be able to 
explain. 


The conception of the isoauroral diagram is obvious when once pre- 
sented, but (so far as I know) it was a novel conception when it first occurred 
to me in the summer of 1950, while I was Research Associate at the California 
Institute of Technology, working on a U.S. Signal Corps contract. 
Dr. Vestine, with whom I corresponded on the matter, while concurring in 


the importance of such a diagram, emphasized the magnitude and difficulty 
of the task of its preparation. 


Since then I have referred to the isoauroral diagram in several lectures 
on the aurora, and also in correspondence with Mr. J. Paton of the Uni- 
versity of Edinburgh (Director of the Auroral Section of the British Astro- 
nomical Association), in connection with his plans for the extension of auroral 
observation in Britain, and also elsewhere, by British air pilots. I have 
particularly stressed the importance of observers recording occasions when 
they see aurore lying or passing overhead; in the case of such zenithal occur- 
rences, the latitude of the aurora is unambiguously determined; when the 
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elevation of the aurora above the horizon is Jess than 90°, its determination 
is of great value for use with any similar determinations made by observers 
situated to the north or south, though it may be difficult always to be sure 
that the measured elevations refer to the same time and same point of the 
aurora; without such simultaneous observation from more than one place, 
the measurement of the elevation (and azimuth) of auroral points does not 
determine the location of the point unless its height is known. The Nor- 
wegian observations of auroral heights—particularly those by Stérmer and 
Harang—show that though there is a normal statistical distribution of 
height (varying, however, with latitude, e.g., as between Northern and 
Southern Norway), there are special types of aurore—not necessarily 
recognizable at sight—whose height distribution departs considerably from 
the normal. 


Next in definiteness to observations of overhead aurore are pairs of 
observations of the same aurora from two places, from one of which it is 
seen to the south, and from the other to the north, indicating that the true 
latitude of the aurora lies between those of the stations: the respective values 
of the elevations above the southern and northern horizons will determine 
the true latitude at least approximately, as well as giving (a perhaps less 
accurate value of) the height. 


Mr. Paton has emphasized some of these points in his circulars to 
observers. It may be hoped that through his efforts, and those of leaders 
in continental Europe, America and the Soviet Union, a more widespread 
network of auroral observers will be developed; the observations of air 
pilots will be of special value, because they are free from the limits imposed 
in the past by cloud cover. In this work the new techniques of radio observa- 
tion of aurore will also play a useful part. In a decade or so, with such 
more continuous and more quantitative observations, it should be possible 
to accumulate a body of auroral statistics surpassing in value all the long 
records of the past—so far as concerns geomagnetic latitudes above about 55°. 
This should permit the construction of the northern isoaurores from the 
auroral zone down to this geomagnetic latitude; and the development of 
civil and military aviation over the northern polar cap will assist in the 
delineation of the inner isoaurores, within the zone. 


The delineation of the southern isoaurores, around the south geo- 
magnetic axis pole or auroral centre, is bound to lag behind, owing to the 
less favourable disposition there of habitable land, and the sparsity of land, 
sea and air traffic south of geomagnetic latitude — 55°. 


F 
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The isochasm approximately coinciding with the circle of geomagnetic 
latitude 55° corresponds nearly to the visibility of about ten aurora in an 
average year, according to Fritz and Vestine. 


THE ENUMERATION OF THE ISOAURORES 


The general nature of the isoauroral diagram can easily be inferred. 
Our present knowledge, especially that embodied in the isochasms, suggests 
that the frequency of auroral occurrence is a function mainly of geomagnetic 
latitude, though the apparent fact that the auroral zone and isochasms are 
not truly circular signifies that geomagnetic longitude is a minor independent 
variable of the auroral frequency function. Hence the isoaurores may be 
expected to have almost if not exactly the same form as the isochasms; to 
this extent the patterns of the two diagrams will be closely similar. The 
difference will lie in the enumeration of the frequencies, isochasmic and 
isoauroral, associated with the same curve. Outside the north auroral zone, 
aurore boreales are most commonly seen to the north; they differ among 
themselves in height (and in consequent range of visibility). At the frequently 
occurring lower limit of height of 100 km. an auroral point can be seen (if 
the sky is clear and the view is open down to the horizon) from a distance of 
over 1,120 km. or about 700 miles, corresponding in latitude to about 11°; 
even though an aurora may not be noticed unless it is elevated at least 5° 
above the horizon, yet, owing to the fact that it is likely to extend at least 
a few tens of km. above 100 km., it will contribute to the isochasmic fre- 
quency over a belt of latitude, to north and south of itself, at least 22° wide- 
If aurore occurred at one definite height, it would be possible by a fairly 
simple mathematical transformation to derive the isoauroral numeration of 
the isochasmic curves from their isochasmic numeration; the actual height- 
frequency distribution of aurore makes this a far more difficult task,— 
indeed, at present it is not feasible, because of the sparsity of reliable height 
determinations for aurore in geomagnetic latitudes below about 60°. 


It is at least clear that for geomagnetic latitudes less than about 65°, 
most of the aurore seen will lie to the north, and the isoauroral frequency, 
of those lying actually above the places along a given isochasm, will be less 
than the isochasmic frequency. The ratio of the isoauroral to the isochasmic 
frequency is likely to decrease from nearly one, for the auroral zone itself, 
to zero at the geomagnetic equator; and similarly it is likely to decrease 
within the auroral zone, with increasing geomagnetic latitude, to some value 


as yet unknown, but perhaps significantly different from zero, at the auroral 
centre. 
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This implies that the number N associated with the isoaurores will 
decrease more rapidly, in passing both northwards and southwards from the 
auroral zone, than the isochasmic number; or in other words, for equal 
increments of N the isoaurores will lie more closely towards the zone than 
do the isochasms. 


6. ISOAURORES AND ISOCORONES 


The isoaurore numbered N is the locus of places above which some 
auroral luminosity is observable on N nights per year over the average of a 
long period of years. This simple conception can of course be elaborated ; 
for example, an isoauroral diagram may be drawn for a particular local 
time, or a particular season, or a particular sunspot epoch; and one can 
imagine a continuous series of diagrams, permitting cinematographic repre- 
sentation of the variation of (average) frequency of overhead aurore, with 
respect to local time, season or sunspot epoch. 


The idea can also be extended in association not with time but with 
auroral character; for example, N may refer to the average frequency of 
auroral coronas observed at the observing stations; in this case the auroral 
location considered is not in the celestial but in the magnetic zenith, whose 
celestial elevation varies with latitude: such lines may be called isocoronal 
or isocorones. Similarly isoaurores may be delineated for red aurore, or 
for auroral arcs crossing the zenith. 


The N corone associated with an isocorone N outside the northern 
auroral zone will be overhead at places somewhat to the south of the iso- 
corone. At these places other (non-coronal) aurore such as arcs may also 
occur; hence at these places the isoauroral number will exceed N. As the 
isoauroral number increases with latitude (up to the auroral zone), its value 
for the isocorone N will a fortiori exceed N. Hence the isoaurore N must 
lie outside the isocorone N. 


It is an interesting subject of enquiry to determine the lowest geomagnetic 
latitude at which auroral corone appear. 


7. IsOAURORES BELOW GEOMAGNETIC LATITUDE 55° 


As suggested in § 4, energetic widespread observation of aurore above 
geomagnetic latitude 55° may in a decade or so enable reliable isoaurores 
to be drawn for the region between this latitude and the auroral zone. 
Below 55°, however, where aurore occur only rarely, a longer period of 
observation is necessary. A future of unknown duration may be expected 
to provide the data for the complete construction of the isoauroral diagram ; 
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but if it is desired to construct a provisional diagram now or in the near 
future, recourse must be had to past records. 


From these I hope soon to evaluate the isoaurores over England down 
to its southern limits, in geomagnetic latitude about 53°. Even so limited 
a task involves much searching among old periodicals and other records. 
My experience leads me to think there may be a wealth of auroral informa- 
tion scattered throughout such sources, including old newspapers which in 
some cases may have long lapsed and been forgotten. This material may 
become more numerous the higher the latitude, up to the limits where settle- 
ment and civilization becomes sparse. Where aurore are common it is 
not worthwhile delving into such sources, but with decreasing latitude and 
auroral frequency, the task becomes more worthwhile, and also more limited, 
in the sense that the dates for which a search is likely to be fruitful are already 
known from existing catalogues of aurore and magnetic storms. 


Of course much of the old records brought again to mind by such a 
search will be vague and unquantitative, and therefore of little value. Here 
and there however, reliable definite items of information will accrue, which 
will help in reconstructing at least some ‘ picture’ of the development in time 
of the position, form and colour of the past aurore, for use in the construc- 
tion of a provisional isoauroral diagram. 


INDIAN AURORZ 


Particular interest attaches to the most extreme aurore, whose range 
of visibility has extended even into the tropics. Owing to the obliquity of 
the geomagnetic axis, such extension is easier in and near the longitudes of 
the auroral centres (70° W. in the north, and 110° E. in the south) than mid- 
way between. India, with a range of longitude from about 70° to 90°E., 
is rather unfavourably placed for auroral observation; the geomagnetic 
(gm.) latitude there is about 10° less than the geographic (gg.), e.g., for 
Bombay, 72°-8 (E.) it is 9°-5 (gm.) as against 18°-9 (gg.), and for Madras, 
further east (80°-2), it is 3°-1 (gm.) as against 13°-1 (gg.). In the north, 
India and Pakistan, including Kashmir, extend to about 25° geomagnetic 
latitude, nearly the same as Tokyo; at the best, therefore, the expectation 
of aurore visible from India is very low. This however, imparts exceptional 
importance to any that have been seen there. 


Fritz’s catalogue lists only two such, and one of these (Verzeichnis, 
p. 160) refers to ‘‘ an aurora-like appearance seen near the western horizon 
from Benares on 1847 November 11; I have now looked up his sources® 
which give the additional information (though without an indication whence 
the report came) that the appearance was like a comet’s tail, and later took 
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a snake-like form and disappeared in 10 minutes. In his Polarlicht, (p. 61) 
Fritz terms this observation ‘‘ somewhat questionable’, and this view is 
supported (in the absence of clearly favourable evidence) oy the fact that in 
general, when aurore appear beyond their usual range of latitude, they are 
exceptionally brilliant and widespread in that range, and are accompanied 
by magnetic storms and strong earth currents. These associated phenomena 
were apparently not observed on 1847 Nov. 11, though on that date (as on 
very many others) aurora was observed at New Haven, Conn., U.S.A., and 
at Jakobshavn (N. Greenland). The Benares report probably refers to 
a somewhat long-lasting meteor trail. 


The other case noted by Fritz (Verzeichnis, p. 244 for aurora borealis, 
p. 255 for aurora australis) was that of 1872 Feb. 4, discussed in §10. This 
case cannot be doubted, and is of great interest. 


Prof. D. S. Kothari in 1949 kindly drew my attention to another report 
of an aurora seen from India, and recorded by Joseph Dalton Hooker.?° 
This report is reproduced in § 8 and discussed in § 9. 


It seems to me somewhat unlikely that the aurora has been altogether 
beyond visibility, at least in northern India, before and since 1872. I think 
it is much to be desired that a future auroral watch be kept, by professional 
and amateur astronomers and naturalists in northern India, and that 
warning should be broadcast when solar disturbance makes it somewhat 
probable that a magnetic storm and a great auroral manifestation may 
occur: so that the occasions when there is some chance of observing another 
Indian aurora may not be missed. Further, I would suggest that the news- 
papers and other periodical literature and records in Lahore and other 
northern centres should be searched for auroral references, particularly 


relative to the dates of the greatest aurore and magnetic storms of the past 
century and a half. 


8. HOOKER’s AURORAL REPORT (pointed out by Dr. D. S. Kothari) 


On the night of 1848 Feb. 14, J. D. Hooker with three companions 
(Messrs. Williams, Haddon and Theobald) were on the east bank of the 
Soane River, at 24° 52’ N., 84° 22’E., altitude 345 feet, in geomagnetic 


latitude 13° 6’ N. The sky was calm and clear, the moon three quarters full, 
and bright. 


At 9 P.M. they saw about 30 “ lancet beams rising in the north-west, 
from a low luminous arch, whose extremes bore W. 20° S., and N. 50° E.; 
altitude of upper limb of arch 20°, of the lower 8°. The beams crossed 
he zenith, and converged towards S. 15°E. The extremity of the centre 
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one bore S. 50° E., and was 45° above the horizon. The western beams 
approached nearest the southern horizon. All the beams moved and flashed 
slowly, occasionally splitting and forking, fading and _ brightening; they 
were brightly defined, though the milky way and zodiacal light could not be 
discerned, and the stars and planets, though clearly discernible, were very 
pale.” 


“At 10 P.M., the luminous appearance was more diffused, the upper 
limb of the arch less defined; no beams crossed the zenith, but occasional 
beams appeared here and faded away.” 


** Between 10 and 11 P.M., the beams continued to move and replace one 
another, as usual in aurore, but disappeared from the south-east quarter, 
and became broader in the northern hemisphere; the longest beams were 
near the north and north-east horizon.” 


** At 10-50 p.m., the dark band had increased so much in breadth that 
the arch was broken up in the north-west, and no beams appeared there. 
Eighteen linear beams rose from the eastern part of the arch, and bore from 
north to N. 20° E.” 

“ Towards 11 p.M., the dark band appeared to have replaced the 
luminous arch, the beams were all but gone, a few fragments appearing in 


the N.E. A southerly wind sprang up, and a diffused light extended along 
the horizon.” 


‘“* At midnight, I saw two faint beams to the north east, and two well- 
defined parallel ones in the south-west.” 


9. DISCUSSION OF HOOKER’s REPORT 


Though Hooker’s book was published in 1854, its mention of observa- 
tion of an aurora clearly remained unknown to Fritz, whose Catalogue does 
not include it; but on p. 161 it notes that on 1848 Feb. 14, aurore were 
observed at Christiania (now Oslo), where Hansteen saw a strong arc between 
9 and 10 p.m. local time (about 5 hours later than Hooker’s first observa- 
tion), and at Montreal and Toronto, Canada. 


But one would expect an aurora; thus visible in India to be a very out- 
standing one, visible at many other places down to and even beyond the 
same low latitude, and associated with a very great magnetic storm. En- 
quiries on the latter point elicited from Dr. S. L. Malurkar, Director of the 
Bombay Observatory, and from the Astronomer Royal, the information 
that 1848 Feb. 14 was not a day of great magnetic disturbance at Bombay 
or Greenwich (at that time, some years before automatic photographic 
recording was introduced, the magnetic instruments were read by eye only, 
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thus no continuous record was available). The Catalogue of Fritz shows 
that there was no worldwide auroral display on that night. 


What, then, should one think of Hooker’s account? He mentions in 
his book (1854) that “ the appearances are as noted in my journal at the 
time”. He added: ‘ They so entirely resembled auroral beams, that I 
had no hesitation at the time in pronouncing them to be such. This opinion 
has, however, been dissented from by some meteorologists, who consider 
that certain facts connected with the geographical distribution of auroras 
(if | may use the term), are opposed to it. I am well aware of the force of 
these arguments, which I shall not attempt to controvert; but for the informa- 
tion of those who may be interested in the matter, | may remark that I am 
very familiar with the aurora borealis in the northern temperate zone, 
and during the Antarctic expedition was in the habit of recording in the 
logbook the appearance presented by the aurora australis. (My compa- 
nions), who were also witnesses of the appearances on this occasion, consi- 
dered it a brilliant display of the aurora”. 


These statements, and the excellent objective quantitative record of the’ 
appearance, so like the best descriptions of undoubted aurora, indicate that 
what Hooker and his companions observed was an objective phenomenon, 
and give some confidence that it was in fact an aurora. On this inter- 
pretation, it was one which was unusual in that the area over which it seems 
to have appeared was limited. As far as | know, there are few if any other 
such well authenticated “ limited-area ’ low-latitude aurore. However, as 
the arch rose to an elevation of 20°, and as the ‘ beams’ (more commonly 
called rays) crossed the zenith into the southern half of the sky, it seems 
almost impossible that they should have escaped all other observation 
(despite the bright moon) over the great area within whose range of visibility 
they must have been presumed to lie. This would almost certainly have 
included Allahabad, Benares and Patna, not improbably Cawnpore, 
Lucknow and Calcutta, and perhaps Delhi. I would suggest that the news- 
papers published shortly after 1848 Feb. 14 in these cities should be 
examined to find any references they may contain to this phenomenon. 


The year 1848 was certainly a year of exceptional auroral frequency, 
and a sunspot maximum year of decidedly high sunspot number; it included 
six (an unusually large number of) crimson aurore, which rarely occur at 
other sunspot epochs (Fritz, Polarlicht, p. 135). This adds to the probability 
that what Hooker observed was a real aurora. 


Incidentally this reference to red aurore calls to mind the one sur- 
prising omission in Hooker’s report, the absence of reference to the colour 
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of the beams and arches; this suggests that his aurora was not crimson, or 
he could scarcely have failed to say so. 


However, this Section raises anew the question of the appropriate 
meaning of the word aurora. The late Lord Rayleigh introduced the term 
nonpolar aurora to signify what is now more often called the night airglow—or 
light of the night sky as seen at times of geomagnetic calm, or uncomplicated 
by the polar aurora associated with the entry of particles into our atmosphere 
from outside. Though patches of special luminosity—sometimes emitting 
the sodium ‘lines "—are at times seen in the night sky, the night airglow 
usually differs greatly in appearance and distribution from the aurora polaris, 
which Hooker’s aurora so strongly resembled, though its apparently isolated 
regional character seems to distinguish it sharply from almost all other 
aurore seen in low latitudes. The only other low-latitude account I have 
come across at all similar to Hooker’s, is one by A. Poey, Director of the 
Havana Observatory, Cuba, 23° 9’ N., 76° 5’ W.; 34°-3 gm. N. On the night 
of 24/5 March 1860 he saw from Havana an aurora to the east, of arc form, 
showing marked undulations and fluctuations, and appearing and disappear- 
ing over an interval of 5 hours—11 P.M. to 4 A.M.:—as recorded in the 
Comptes Rendus, 1860. Fritz mentions this in his Verzeichniss, but records 
to other aurora seen elsewhere on that night.—Dr. C. T. Elvey informs me 
that he will shortly publish an account of some unusually definite forms of 
the night airglow seen from Texas.—Added 14 Jan. 1953. 


10. THE INDIAN AURORA OF 1872 FEB. 4 


The aurora of 1872 Feb. 4 was certainly one of the greatest on record; 
Fritz in his Catalogue (p. 244) mentions two places in India where it was 
observed, Bombay (gm. lat. 9°-5) and Lahore (31° 34’N.; 74° 21'E,; 
gm. lat. 22°-1). He adds that Jannsen was observing in Ceylon (7° N. gg.) 
that night, and does not mention seeing an aurora. In Polarlicht (pp. 32, 
151) Fritz mentions also Raikote (30° 39’ N.; 75° 35’ E.; gm. lat. 21°-1) 
with Bombay and Lahore. He does not give his sources for these Indian 
observations; Nature (1872, 5, 323) quotes from a letter of Feb. 5 by 
George Draper, of the British Indian Submarine Telegraph Company, that 
the aurora was visible at Bombay, but without details of the Bombay auroral 
observations; and I have been unable to find (as yet) any useful details about 
the Indian observations at the above three places, in any literature accessible 
to me in England. 


However, at my request, Dr. S. L. Malurkar kindly made a search in 
the great Bombay newspaper The Times of India for some days after 1872 
Feb. 4, and supplied me with four extracts from the issues for (a) Feb. 6, 


187 


(b) Feb. 12, and (c,d), Feb. 15; for Feb. 15 there were two reports, one 
of them (c) being taken from the Delhi Gazette for Feb. 12. I have no 
space to reproduce these extracts here; I hope to use them in a future study 
of this aurora, together with a great amount of material from other parts 
of the globe. But I may mention that all four reports contained useful 
information relating to observations from Bombay, Rawal Pindee, Madhopore 
and Jacobabad. They also mentioned that the aurora had been visible 
from the following other Indian places, Lahore, Mooltan and Sukkur, and, 
outside India, from Aden. (The geomagnetic latitude of Aden, 9°, seems to 
be slightly less than that of Bombay, and the Aden observation, thus noted 
from a newspaper, constitutes, as far as I yet know, the record extension 
of auroral visibility towards the geomagnetic equator.) 


Professor S. K. Mitra has at my request kindly sought for Calcutta 
references to this aurora; the Journal of the Asiatic Society for 1872 and 
1873 contains no mention of this aurora. The Statesman files in Calcutta 
are incomplete before 1875, but I hope may be searched from London files. 
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{ would suggest to scientists of India that it is worthwhile to continue 
the search for Indian observations of this aurora, south as well as north of 
Bombay, with the aim of determining the southern limit of its visibility over 
India, and all quantitative details as to geometry and times, that may contri- 
bute to our knowledge of this remarkable event. I hope similar searches 
will be made in the records of other countries, and in nautical records, for 
low-latitude observations of this aurora (although a number of such are 
already in hand). 


In addition, similar Indian data for other great aurore might be looked 
for, with particular reference to the following :— 
1859 Aug. 28-Sep. 4. 
1882 Nov. 17-21. 
1903 Oct. 31-Nov. 1. 


1909 Sep. 25. 
1921 May 13-16. 
1938 Apr. 16. 


I hope also that a lookout for future auroral appearances visible in India 
may be kept, as suggested in § 7. 


SUMMARY 


(1) After describing the isochasmic diagrams of Fritz and Vestine, which 
map the distribution of frequency of auroral visibility from any point on 
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the earth, it is pointed out that there is a need for a new type of isoauroral 
diagram, showing isoaurores or lines of equal frequency of occurrence above 


any terrestrial point. The nature and “‘ cnumeration”’’ of the isoaurores 
is discussed. 


(2) The few aurore which have been reported as visible from India 
are discussed: namely (1) one reported from Benares on 1847 November 11, 
(2) one reported by J. D. Hooker as seen on 1848 Feb. 14, from a point on 
the Soane River at 24° 52’ N., 84° 22’ E..—the validity of both these reports 
has been questioned—and (3) the undoubted case of the aurora of 1872 
Feb. 4, seen as far south as Bombay, and from many more northerly points 
in India and elsewhere. 


(3) An appeal is made for a search by Indian scientists for material 
relating to Indian observations of these and possibly other aurore. 
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THE importance of the vertical eddy flux of momentum for the dynamics 
of the Trades has recently been emphasized by Riehl, et a/.t and by Sheppard 
and Omar.’ These investigations show that the wind veer with height in 
the lowest layer is much smaller than what could be expected from the 
classical conception of the frictional boundary layer. Thus the clockwise 
angular turning of wind with height in the Trade-wind belt of the eastern 
Pacific Ocean is, according to Riehl, et a/., only about 4 degrees up to the 
900-mb. surface, whereas the corresponding turning according to Sheppard 
and Omar amounts to about 7 degrees from surface to 2,000 ft. at Johnson 
Island (16° 44’ N.; 169° 31’ W.). Since the total wind velocity increases 
with height at the rate 2-3 m. sec.-! in the layer in question the cross-isobaric 
wind component remains practically constant in the frictional layer. 


If one applies these results on hemispheric scale the cross-isobaric wind 
component corresponds to the mean meridional component v, whereas the 
component parallel to the isobars corresponds to the mean zonal component u. 
The variation of the zonal component of the shear stress +, with height is 
determined by the equation 


po, = 2 sin + (1) 


Here p is the density of air, w is the angular velocity of the earth’s rotation 
and ¢ is the latitude. In this equation the stress due to horizontal shear 
has been neglected. This latter is important when considering large-scale 
disturbances. According to computations by Widger’ and Mintz! the 
horizontal flux of momentum due to large-scale eddies, however, seems to 
be relatively unimportant in the lowest part of the atmosphere discussed here. 


Integrating equation (1) between the earth’s surface and an arbitrary 
level H one gets: 


H H 
f p vdz + f pat dz (2) 
0 
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This equation determines the flux of westerly momentum into the layer 
between the surface and the level H. Since the vertical flux of momentum 
is determined by 


the flux of westerly momentum is directed upwards in the layer where the 
east wind increases with height and downwards where it decreases. At the 


level where be == 0 there is no vertical eddy flux of momentum. If this level 


is denoted by H, the surface stress can be determined from the equation 
Hy 
Tx9 = 2w sing f p vdz, (4) 
0 


if the acceleration term in equation (2) can be neglected. 
The surface stress can also be computed from the empirical formula 
Tx,0 = kpocy” COS a, (5) 


Here k is the drag coefficient, cy is the velocity of the wind at the surface 
(‘“‘ anemometer height”) and a is the angle between wind direction and 
isobars. The negative sign indicates an upward transport of westerly 
momentum from the earth. Thus there are two possibilities to determine 
the surface stress; it can be computed from equation (4), if the height H, 
and the value of v between surface and that level are known, or it can be 
computed from formula (5), if the drag coefficient .k and the surface wind 
are known. 


At the latitude of no zonal wind, corresponding approximately to the 
latitude of highest surface pressure, the level H, sinks down to the surface. 
Here the surface stress vanishes both according to formule (4) and (5). This, 
however, does not necessarily mean that v vanishes in the surface layer. For 
the Trade-wind zone it is not possible to fix a constant H,. According to 
Sheppard and Omar the level varies between about 300 and 800 m., whereas 
according to Riehl, et al., it could be estimated at around 900 mb. in the 
Trade-wind zone of the eastern Pacific Ocean. It is obvious that the level 
of maximum east wind varies both with latitude and wind velocity. In the 
following a probable mean value of 700 m. will be used for the latitude of 
strongest east wind. 


The zone of maximum east wind during winter lies somewhere in the 
latitude belt 10°-15° N. In this belt the acceleration term in equation (2) 
is small and can be neglected when the equation is applied to the lowest 
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part of the atmosphere. The maximum southward component can accord- 
ing to Riehl and Yeh* be observed around the same latitudes. Assuming 
with Riehl and Yeh that v has the value — 2:5 m. sec. at latitude 13° N. 
and that the angle between surface wind and isobars, or parallel circles, is 


km 


Fic. 1. Schematic distribution of zonal wind component, u, and meridional wind compo- 
nent, v, with height in the zone of strongest east wind in the Trades. H, denotes the height of 
strongest east wind, whereas H, is the level where the meridional component vanishes. 


25° one gets profiles of the wind components of the type presented in the 
figure. Referring to the previous discussion and especially to the wind 
profiles given by Sheppard and Omar the level of maximum east wind, H,, 
is put at 700 m., whereas the level of no meridional wind, H., is assumed 
to be 3,000 m. This latter value seems to be, at least partly, supported by 
results presented by Starr and White. It might be pointed out that the 
level of maximum zonal wind must be somewhat higher than the level of 
maximum wind component in the direction of the surface wind. For this 
latter Sheppard and Omar gave a mean value of about 500m. In the figure 
it was further assumed that the increase of west wind with height in the free 
atmosphere amounts to 3 m. sec.-! per km., corresponding to a meridional 


temperature gradient of about 0°-3C. per 100 km. 


Using the above values one gets for the surface stress 7,9, according 
to equation (4), the value 0-64 dyne cm.-*, corresponding to a value for the 


coefficient k in (5) of 0:0017. This value should be compared with the 
A3 


3H, 
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H, 
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somewhat lower value of 0-0012 presented by Sheppard and Omar. The 


difference depends upon the value used for H,. With H, = 500 m. the drag 
coefficient would have been 0-0012. 


It is important to point out that the drag coefficient depends upon the 
roughness of the ocean surface. Computations of the coefficient k by use 
of other methods (Sverdrup’), e.g., over the Baltic Sea have given values 
around 0-:0025. It seems quite possible that the roughness of the sea surface 
in the Trade-wind belt with its long, fast moving waves is much less than 
in the Baltic area with its shorter and more irregular waves. Since the actual 
surface stress over the oceans is of greatest importance for several oceano- 


graphic and meteorological problems one must regret that the drag coefficient 
is not better known. 


Returning to equation (2) it is now possible to estimate the shear stress 
in the free atmosphere. If we still neglect the acceleration term the shear 
stress increases with height up to the level where v becomes positive. The 
increase of the shear stress in the 700 m. deep layer between surface and H, 
would be 0-64 dyne cm.-* at the latitude 13° N. in January. The follow- 
ing table gives the zonal shear stress at different heights computed from the 
figure. 


TABLE 
Zonal Stress at Different Heights at Latitude 13° N. in January 


Height m, | Sarface 500 | 1000 1500 | 2000 | 2500 3000 3500 | 
Zonal stress, dyne cm.~? | ~—0°64 | —0-18 | 0-3! 0-66 0-88 | 0-99 | 1-02 0-98 | 0-89 
| | 


The level H, in the figure, where v = 0, is the level of maximum shear 
stress or maximum downward transport of momentum. In the figure H, 
is assumed to be 3 km., or about 700 mb. The shear stress here is according 
to the table about 1 dyne cm.-? According to this result the flux of momen- 
tum downward through the 700-mb. level thus would be larger than the 
upward flux from the earth’s surface. 


The shear stress can be expressed by a coefficient of eddy viscosity 
according to equation (3). If the increase of west wind with height at the 
level 3 km. is 3m. sec.-! km.—', or = 3-10-* sec-1, as assumed in the 


figure, the corresponding value of the coefficient of eddy viscosity p is 340 g. 
cm.-!sec.-! For the eastern part of the Pacific Ocean, between the latitudes 
20° and 30°N., Riehl, et a/., got values between 370 to 680 g. cm.— sec.-}, 
with the highest value around the 800-mb. level. It seems probable that 
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the level of maximum » should increase somewhat southwards due to the 
increase of the height of the cloud layer. 


In the computation of the vertical transfer of momentum at a given level | 
H the acceleration term in equation (2) was neglected. That is permissible 


in the lowest part of the atmosphere outside the equatorial zone since Fa 


here generally is much smaller than |2w sin dv|. This can easily be seen, 
if we write equation (2) in the form 


(2e sin dv — v  — w dp, (6) 

pH 
Here g is the acceleration of gravity, p is the pressure and w denotes the 
vertical wind component. Considering a relatively shallow layer one has 
only to compare the possible magnitude of the three terms in the integral. 
The two last terms are of about the same order of magnitude. The hori- 


zontal shear < 5x10-* sec!, or of the same order of magnitude as 
the Coriolis parameter in the belt between the equator and 4°N. At lati- 
tude 20° N. E 5 | < one-tenth of the Coriolis parameter. If the integration 


is extended higher up the acceleration cannot any more be neglected, since 
the relative importance of the first term gradually decreases with height, 
whereas that is not generally true with the other terms. 


In the upper part of the atmosphere the horizontal eddy flux of momen- 
tum, which has been neglected in the previous discussion, becomes important, 
as, e.g., Mintz? and Starr and White* have shown. This type of eddy flux 
of momentum is, however, primarily connected with very large-scale dis- 
turbances not considered here. 


In a paper by E. Palmén and M. A. Alaka recently submitted for publi- 
cation in Tellus? a different method was used in estimating a complete 
budget of angular momentum in the zone 20-30°N. The stress exerted by 
the atmosphere at 700 mb. was here computed to about 0-5 dyne cm.~2, or 
approximately half of the value estimated above for the latitude 13° N. The 
difference could partly be explained as the result of difference in latitude, 


but depends essentially upon the values!for the surface stress used in the 
computation. 


The data available are still much too unsatisfactory to permit quanti- 
tative computation of the shear stress on a hemispheric scale. All our 
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knowledge, however, indicates that there is a considerable downward eddy 
flux of momentum in the Trade-wind zone and that this vertical flux of 
momentum should not be neglected in any study of the dynamics of the 
Trades. 
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Dés que furent connus les premiers résultats sur la répartition verticale de 
Yozone!: 2 indiquant la présence de ce gaz dans les régions trés froides de 
latmosphére, on se préoccupa de déceler l’effet de cette basse température 
sur le spectre d’absorption de l’ozone atmosphérique; les mesures de 
laboratoire avaient en effet montré une modification de certains coefficients 
d@absorption dans la région des bandes de Huggins lorsque l’ozone est 
refroidi. En 1935, Barbier, Chalonge et E. Vassy,* 4 partir de spectres 
stellaires obtenus 4 la station de Jungfraujoch (Suisse), mirent en évidence 
leffet de la basse température de la stratosphére sur le spectre d’absorption 
de l’ozone atmosphérique. Leurs résultats montrérent la nécessité de 
nouvelles mesures de laboratoire sur l’influence de la température sur les 
coefficients d’absorption de l’ozone. Ce travail effectué par l’un de nous, 
le conduisit 4 définir la notion de température moyenne de l’ozone atmo- 
sphérique*; il montra que les spectres de l’ozone atmosphérique (sous réserve 
de certaines conditions expérimentales) permettent de déterminer cette 
température moyenne 4 quelques degrés prés, ce qui est suffisant pour un 
bon nombre d’études géophysiques ou météorologiques. 


I. DEFINITION 


Cette définition de la température moyenne de l’ozone atmosphérique 
repose sur le fait expérimental suivant: pour les longueurs d’onde des 
minima d’absorption des bandes de Huggins, les coefficients varient suivant 
une loi trés approximativement linéaire, et en tous cas d’une facon univoque, 
entre — 80°C. et + 20°C., c’est-a-dire dans les limites de température que 
lon peut rencontrer dans la stratosphére et la mésosphére. Les maxima 
par contre, 4 quelques exceptions prés, ne varient pas, circonstance 
heureuse qui permet, quelle que soit la température, de déterminer sans 
erreur l’épaisseur réduite. La faible variation des maxima observée par 
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Barbier et Chalonge® puis par Vigroux® ne changerait en rien notre méthode 
de détermination de la température. 


Considérons une longueur d’onde A, le coefficient d’absorption fonction 
de la température est k,>. 


Supposons latmosphére divisée en couches minces d’épaisseur e, et de 
température T,,. 


L’absorption due 4 |’ozone total est 


appelons E l’épaisseur toiale Xe; on peut définir un coefficient d’absorption 
moyen par la relation: 


Sek.’ = Ek,* (2) 
k, est la valeur que l'on obtient expérimentalement 4 partir des spectres 
d’absorption de l’ozone atmosphérique (voir Fig. 1). 


Mais tenant compte de la loi de variation trouvée, nous pouvons écrire, 
avec une trés bonne approximation: 


kp=a+6T (3) 
et par suite: 
2ek, = [e (a+ 
=aE+ b&eT (4) 
Or on a de méme: 


Ek, = E(a+ 
En substituant (4) et (5) dans !’équation (2), nous en tirons: 
ZeT = 


(6) 
Cette équation définit une température moyenne de ozone atmosphérique, 8. 


Cette température s’obtient facilement a partir du coefficient d’absorption 
moyen kp. 


Sa définition par l’équation (6) permet de la calculer simplement si |’on 
connait les valeurs de e et T aux différentes altitudes dans l’atmosphére, ou 
de déduire l'une de ces données si l’on connait les autres; nous en verrons 
un exemple plus loin. 


II. CONDITIONS ET DIFFICULTES DE LA MESURE 


La mesure expérimentale de cette température moyenne exige un certain 
nombre de précautions et de restrictions. Elles sont dues aux deux fait, 
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suivants: premicrement, la variation du coefficient d’absorption est faible 
ou nulle dans certaines régions spectrales; deuxiémement, les bandes de 
Huggins ont une structure assez serrée, deux maxima étant séparés par 
25 4 30 Aen moyenne. Ceci entraine un certain nombre de limitations et 
de difficultés dans l’utilisation de la méthode. L’expérience a montré que 
le coefficient d’absorption varie en fonction de la température plus ou moins 
rapidement suivant la longueur d’onde; donc certains minima sont plus 
sensibles 4 cet effet, et par suite devront, dans la mesure du possible, étre 
utilisés de préférence aux autres pour la détermination de la température 
moyenne. Ainsi tous les minima de longueur d’onde plus courte que 3151 A 
ont une variation trop faible, le rapport des coefficients 4 + 20° C. et — 80°C. 
étant compris entre 0,85 et 1,00, alors que pour 3359 par exemple il est de 
0,40. Ces courtes longueurs d’onde ne peuvent étre utilisées. 


Pour les longueurs d’onde plus grandes que 3295 A, l’effet de tempé- 
rature est intéressant, mais le coefficient d’absorption est faible, ce qui 
entraine une faible variation en valeur absolue avec la température; comme 
la pratique des mesures d’ozone atmosphériques exige souvent, par 
exemple pour éviter de rencontrer l’effet Gdtz’ lorsque lon opére sur la 
lumiére solaire diffusée, l'utilisation de masses d’air assez petites (sec Z 


inférieur 4 3) la précision de la mesure de température sera trés faible dans 
ces cas au-dessous de 3300 A. 


D’une facon générale, pour la méme raison, l'utilisation de masses 
d’air trop faibles au moment de la mesure conduira 4 une erreur importante 
dans la détermination de la température. 


Les températures élevées de la mésosphére soulévent une difficulté d’un 
autre genre: les mesures de coefficients d’absorption 4 haute température 
effectuées par l’un de nous, puis par Vigroux,® font apparaitre une variation 
trés importante 4 la fois des maxima et des minima au-dessus de 18° C. 
et lorsque l’on rencontre des températures moyennes supérieures 4 18°C. 
et nous en avons observé, la détermination en devient délicate. 


Enfin la question du pouvoir séparateur de l’instrument employé est des 
plus importantes; l'emploi d’un spectrographe de dispersion suffisante est 
nécessaire (nous utilisons couramment un instrument qui a, dans les bandes 
de Huggins, une dispersion de 28 A au mm.). En effet, les maxima ne 
présentant pas d’effet de température, on congoit facilement que les varia- 
tions du coefficient d’absorption sont limitées au minimum lui-méme et 4 
une étroite région voisine; il faut donc que la mesure porte exactement sur 
le minimum, et un pouvoir séparateur de l’ordre de l’angstrém est indis- 
pensable; l'utilisation d’un domaine spectral plus large ne peut pas conduire 
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4 observation de l’effet de température. R. H. Kay® a tenté d’étudier cet 
effet en utilisant un spectrophotométre Dobson et il a abouti 4 un échec, 
Cet échec est incontestablement di 4 la largeur des bandes spectrales adoptées 
dans cet instrument; Normand et Kay® ont indiqué que les fentes corres- 
pondent 4 des largeurs de 10 A vers 3100 A (S,) et 30A vers 3300 A (S,); 
cette fente S, laisse passer exactement 2 maxima et la variation du minimum 
intermédiaire se trouve masquée par son étalement sur une trop grande 
largeur spectrale. A l’heure actuelle, seules les mesures avec un spectrographe 
assez dispersif peuvent permettre de déterminer la température. 


Encore se trouve-t-on souvent devant des circonstances aggravantes. 
Etant données les faibles différences d’absorption 4 mesurer, autour de 
0,05, nous améliorons la précision en utilisant plusieurs spectres; mais on 
observe, comme nous le verrons plus loin, des variations accidentelles ou 
systématiques au cours d’une journée. A moins d’avoir multiplié le nombre 
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Examples de détermination de l’épaisseur d’ozone; on remarquera que les points correspondant 
aux minima se placent en-dessous de la droite. 
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Examples de détermination de la température moyenne pour les spectres correspondant a la 
figure 1; chaque point correspond a une longueur d’onde. 


de spectres, la dispersion des résultats sera trop grande pour pouvoir donner 
une valeur de la température. Dans les mesures courantes, il faut compter 
sur une incertitude de 5° C. environ. Méme si cette erreur peut paraitre 
importante, il y a l4 mieux qu’un ordre de grandeur. Et de plus les varia- 
tions observées de la température moyenne de l’ozone atmosphérique sont 
telles que ce renseignement présente beaucoup d’intérét en géophysique. 
Nous allons le montrer par quelques exemples et nous verrons que les 
résultats des mesures se trouvent en bon accord avec les autres données 
météorologiques dont on peut quelquefois disposer en liaison avec les mesures 
dozone, ou les mesures spectroscopiques effectuées dans l’infrarouge. 


II]. QUELQUES EXEMPLES DE MESURES 


Les premiéres évaluations de température moyenne de !’ozone, E. Vassy,* 
Déjardin,!° furent données 4 partir de spectres pris aux latitudes moyennes 
(Jungfraujoch, Mont-Blanc, Montpellier); ces deux auteurs ont trouvé 
— 30°C. et — 35°C. Ce résultat, une fois acquis, a permis 4 E. Vassy* de 
donner un ordre de grandeur de la température de la stratosphére au-dessus 
du plafond des ballons-sondes; en effet on connaissait la répartition verticale 
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de Pozone!s* et les sondages fournissent celle de la température jusqu’a 
30 Kim environ. On a pu déduire que la température de l’atmosphére entre 
30 et 50 Km s’élevait fortement. Ce résultat, déja annoncé par Lindemann 


et Dobson en 1922 n’est plus contesté aujourd’hui et a été précisé 4 l'aide 
d’engins autopropulsés. 


En 1938, nous avons effectué au Maroc, A Ifrane (Lat. 33° 31’ N.) une 
série de mesures systématiques de la température moyenne de l’ozone 
atmosphérique!! échelonnées sur une période de 1 mois d’été; nous avons 
mis en évidence des variations journaliéres de la température moyenne, avec 
comme valeurs extremes — 7°C. et — 49°C. En possession de l’analyse 
de la situation météorologique pour les journées correspondant a nos mesures, 
nous avons pu tirer les conclusions suivantes: 


1° la température moyenne de |l’ozone est liée 4 l’origine de lair; on 
avait déja établi que l’épaisseur réduite est en liaison avec l’origine des masses 
d’air; la température vient confirmer cette permanence des caractéres d’une 
masse d’air pendant son déplacement au cours d’une invasion. L’ozone 
a montré de plus que ces déplacements de masses d’air, observés dans la 
troposphére par les météorologistes, s’étendent 4 la stratosphére. Ainsi 4 
Ifrane situé 4 une latitude de 33°, nous avons observé des invasions d’air 


polaire (épaisseur forte, température relativement élevée) et d’air saharien 
(faible épaisseur, basse température). 


2° les mesures d’épaisseur réduite et de température moyenne de l’ozone 
ont montré comment, pendant sa stagnation, l’air provenant d’une autre 
latitude dégénére et perd réguliérement et lentement ses caractéres d’origine. 


3° Grace au profil en coin de la masse d’air qui s’approche, l’ozone 
permet de déceler a l’avance les mouvements de cette masse d’air; ainsi les 
invasions d’air saharien affectent la stratosphére, c’est-a-dire l’ozone, 2 jours 
avant l’arrivée au sol. L’ozone s’est révélé dans cette étude un excellent 
indicateur des mouvements des masses d’air stratosphériques. 


Si nous nous plagons a une échelle plus grande, l’ozone nous apporte 
des renseignements sur un certain nombre de phénoménes géophysiques et 
en particulier sur l’influence du rayonnement solaire sur la stratosphere Ainsi 
nous avons montré!* comment a Abisko, lat. 68° 20’, variait la température 
moyenne de l’ozone atmosphérique: entre les mois de janvier et d’aout cette 
température passe de — 80°C. a + 20°C. et elle est en liaison avec le 
rayonnement solaire regu 4 haute altitude. Une telle variation annuelle de 
température, qui parait propre aux régions polaires, et doit étre attribuée 4 
Vabsence des alternances réguliéres du jour et de la nuit, doit certainement 
influer sur la circulation générale. 
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Revenant aux basses latitudes, nous avons fait en collaboration avec 
Khalek?* une année de mesures de |’épaisseur réduite et de la température 
moyenne de l’ozone 4 Kaboul (34° N.) qui montrent quelles masses d’air 
régnent en altitude au-dessus de cette région ot alternent les invasions d’air 
sibérien et d’air tropical. La comparison des résultats de Kaboul avec 
ceux d’ifrane dont la latitude est trés voisine promet des résultats 
intéressants. 


Une autre étude faite en méme temps 4 Kaboul a permis de compléter 
nos connaissances sur les variations journaliéres de l’ozone et la dépression 
de midi observée par Karandikar."* Nous avons retrouvé cette diminution 
de l’épaisseur autour de midi et nous avons pu étudier sa variation saisonniére, 
ainsi que l’élévation de température qui l’accompagne. 


1V. INTERET DES MESURES DE TEMPERATURE MOYENNE DE L’OZONE 


Nous avons choisi des exemples aussi variés que possible pour montrer 
combien l’étude de la température moyenne de l’ozone est fertile en découver- 
tes et en renseignements. 


Elle est en elleeméme une donnée intéressante puisque sa variabilité 


dun jour a l’autre correspond aux changements de nature des masses 
d’air. 


Elle permet d’atteindre la température des couches les plus élevées de 
ozone lorsque l’on posséde en méme temps sa répartition verticale, car la 
connaissance de la température ne dépasse pas 30 Km, plafond des ballons- 
sonde, dans les mesures courantes. Nous avons dit que les premiers 


résultats obtenus avec les fusées avaient confirmé les résultats fournis par 
lozone4 


Profitant des connaissances acquises grace aux fusées, qui ont fourni 
une distribution moyenne de la température” et plusieurs répartitions verti- 
cales de l’ozone!® au-dessus de New-Mexico, nous avons calculé la tempé- 
rature moyenne de l’ozone, 9 de l’équation (6), 4 partir de ces données; le 
calcul donne — 43° C. pour le 14 Juin 1949, épaisseur 1,9 mm., et — 50°C. 
pour le 10 Octobre, épaisseur 2, 38mm. La latitude de New Mexico étant 
32° 24’ N., ces valeurs, du moins la premiére qui est relative a l’été, doivert 
ctre comparées 4 nos mesures d’Ifrane (30° 31’ N.), et nous voyons que, 
mises 4 part les journées d’invasion d’air polaire, cette valeur de — 43°C. 
correspond aux résultats que nous avons trouvés. 


Les mesures spectroscopiques sur l’ozone sont donc un moyen simple 
davoir une donnée météorologique intéressante. N’oublions pas que les 
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bandes de Huggins ne sont pas la seule région spectrale permettant d’atteindre 
la température moyenne de l’ozone. La bande d’absorption infrarouge 3 
9, 6u utilisée d’abord par J. Devaux,'’ puis par Adel,'* a donné a cet auteur 
des températures moyennes de — 44°C. pour le 5 Septembre, — 53°C, 
pour le 3 Novembre; ces valeurs obtenues par 35° de latitude Nord, sont 
elles aussi tout 4 fait comparables aux valeurs rappelées ci-dessus. 


CONCLUSION 


Au cours de nos mesures de la température moyenne de !’ozone, nous 
avons été amenés 4 attirer l’attention sur l’insuffisance de nos renseignements 
sur la distribution verticale de l’ozone. II serait intéressant d’avoir des 
données plus nombreuses qui feraient apparaitre des variations avec la lati- 
tude et la saison. Ainsi le travail déjéa commencé par Ramanathan’ a 
fourni plusieurs distributions différentes. De méme, les différents tirs de 
fusées!® ont donné aussi des courbes différentes; on remarquera en parti- 
culier la courbe du 14 Juin qui indique une présence notable d’ozone jusqu’a 
70 Km; jusqu’A maintenant, on tendait 4 admettre que la quantité d’ozone 
au-dessus de 50 Km. était négligeable; cette idée est sans doute a réviser 
et les valeurs élevées de la température que nous avons obtenues tendraient 
aussi 4 faire admettre l’existence d’ozone dans les régions de température 
relativement haute de la mésosphére. 


SOMMAIRE 


Dés que l’on eut reconnu sur les spectres d’absorption de 1’ozone 
atmosphérique les modifications que fait subir 4 ce spectre une basse tempé- 
rature, on s’efforca d’en tirer parti. E. Vassy en 1936 donna une défi- 
nition de la température de l’ozone qui est sommairement rappelée ici. On 
insiste sur les conditions nécessaires d’une mesure correcte de cette tempé- 
rature moyenne a partir des spectres: choix de masses d’air moyennes, 
utilisation des longueurs d’onde les plus sensibles 4 l’effet de température, 
pouvoir séparateur de l’instrument suffisamment élevé. 


On rappelle quelques applications de la méthode qui ont déja donné des 
résultats intéressants: évaluation de la température moyenne aux latitudes 
des régions tempérées, comparaisons avec la méme mesure dans I’infrarouge, 
variations journaliéres et saisonniéres aux basses latitudes, variations annuelle 
dans les régions polaires. 


On termine par une comparaison avec les mesures directes fournies pat 
les fusées et lon montre l’intérét des mesures de température moyenne de 
Y’ozone atmosphérique dont la simplicité permet un emploi étendu. 


AY Pen 


Temperature Moyenne de L’ozone Atmospherique 


Gotz, Meetham ect Dobson 
E. et V. A. Regener 


Barbier, Chalonge et E. Vassy 


E. Vassy 

Barbier et Chalonge 
Vigroux 

Gotz 

Kay 


. Normand et Kay 
. Déjardin 
. A. et E. Vassy 


. Khalegq et A. Vassy 

. Karandikar 

. Havens, Koll et La Gow 

. Johnson, Purcell, Tousey et 


Watanabe 


. Devaux 
. Adel 
. Ramanathan et Karandikar 


BIBLIOGRAPHIE 


Proc. Roy. Soc., 1934, 145, 416. 

Physik Zeits., 1934, 35, 768. 

Rev. Opt., 1935, 14, 425. 

Ann. Phys., 1937, 8, 679. 

Ibid., 1942, 17, 280. 

C. R. Acad. Sc., Paris, 1950, 230, 2170. 
Zeits. Astrophysik, 1934, 8, 267. 


“Commission on atmospheric ozone, U.G.G.I.,”’ 
Bruxelles, 1951. 


Jour. Sc. Instr., 1952, 29, 33. 
Jour. Phys., 1937, 8, 1509. 


C. R. Acad. Sc., Paris, 1938, 207, 1232; La Météo- 
rologie 1939, 3. 


Tellus, 1950, 2, 69. 

C. R. Acad. Sc., Paris, 1952, 235, 737. 
Proc. Ind. Acad. Sci., 1948, 28, 63. 

J. Geophys. Res., 1952, 57, 59. * 
Ibid., 1952, 57, 157. 


C. R. Acad. Sc., Paris, 1935, 201, 1500. 
Astrophysic Jour., 1947, 105, 406. 


A.T.M.E., Transactions of the Oslo Meeting, 1948, 
492. 


203 

i 
a 
le 13 
i- 14 
! 
18 

ne 
é- 
fi- 
yn 
2S, 
re, 
es 
es 
lle 
ar 
de 


COMMENTS ON THE VORTICITY EQUATION 


By JACQUES-M. VAN MIEGHEM 
(University of Brussels) 
Received December 17, 1952 


1. IN order to avoid confusion in the attempts made to explain the occur- 
rence of large values of the absolute vorticity (twice the Coriolis parameter 
Or more) in the atmosphere, a clear distinction between individuai and local 
changes of the absolute vorticity should be made. 


Using a moving system of curvilinear co-ordinates, we shall establish 
analytical expressions for the rate of individual and local changes of the 
absolute vorticity, suitable for application to synoptic meteorology. 

The equation of balance of the absolute vorticity,’ giving explicitly the 
local variation of the absolute vorticity per unit volume, in terms of a vorticity 
flux and a vorticity production, seems to be the most convenient form of 
the vorticity equation for synoptic purposes. 

2. Let (x', x®, x*) be the general co-ordinates in a frame moving arbi- 
trarily with respect to the absolute frame, y, = y;; (i, j = 1, 2, 3) the coefii- 
cients of the metric form in the variables x’, y” the algebraic minor of y, 


dx’ 
divided by the determinant y = || = 


pectively the contravariant and covariant components of the air velocity y 
with respect to the relative frame, w the velocity of the relative frame with 
respect to the absolute frame, V = w + v the absolute velocity of the air 
(absolute wind) defined by the contravariant V’ and covariant V; =y,V' 
components in the moving frame x'x*x3, u the air velocity with seed to 
the earth (wind), w the angular velocity of the carth, & the contravariant 
components of the absolute vorticity curl V = 2+ curl u, p the pressure, 
a the specific volume of the air, s* the contravariant components of the 
solenoidal vector s = Vp x Ya and /* the contravariant components of 
the curl of the frictional force F (f = curl F). The system of co-ordinates 
(x!, x, x*) is supposed to be right-handed, so that, in the northern hemi- 
sphere, a positive vorticity is cyclonic and a negative one anticyclonic. We 
recall the formulas** 


Vy 
dV, dV, da op da 
dx” dx” dx’? Dax?” 
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where the set (i, j, K) is a cyclic permutation of the indices (1, 2, 3). 


It should be noted that the functions y,;, »/ and y are independent of 
the time t when the movement of the relative frame x!x?x* is a solid motion 
(for instance, when the frame x!x*x* is at rest with respect to the earth). 

In the case of general co-ordinates such as (x!, x®, x*), the vorticity 


equations are easily deduced from V. Bjerknes” circulation (or vorticity) 
theorem, we found?’ : 


where a. yx’ Tepresents the individual derivative along the air 


motion. 


We recall that & : \/, represents the orthogonal projection of the 
absolute vorticity curl V on the normal to the co-ordinate surface x* = con- 
stant at the point and at the time considered.* 


3. The transformation of the vorticity equations (1) by substitution of 
the continuity equation 


Ido _ 1 frvy (vw) 
a dt div vv= | 


converts these equations in a form appropriate to the study of individual 
vorticity changes, namely, 


Let us now suppose that at each time ¢, the equiscalar surfaces 4 = con- 
stant constitute a one parameter family of surfaces and let us identify these 
surfaces with the co-ordinate surfaces x* = constant and the orthogonal 
trajectories of + = constant with the co-ordinate lines x* = variable. Then, 
we have: 


= = = = 0, = 1, v = [¥11 Yee — 


v3 
V 733 curl, V = $*, (3) 
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where v, and curl, V are the orthogonal projections respectively of the 
relative velocity v and of the absolute vorticity curl V on the normal n to 
the equiscalar surfaces x* = = constant at the time ¢ and at the point 
considered. The positive » axis is chosen as to point in the direction of 
increasing values of #. 


Substituting (3) in (2), for k = 3, we obtain: 


1d 
a at 


where s, and /,, are respectively the number of isobaric-isosteric solenoids 
and the number of vortex unit-tubes of the frictional force per unit area in 
the surfaces 4 = constant. 


[a| curl, VJ = Ve.curl V + | 74|-(s, +F,,), (4) 


Finally, it is easy to see that 
| curl, V=Ve.curl V, | s,=| Va)n = Va) 


and | = F; hence, Eriel’s vorticity equation*: 


4 (aVp-curl V) = Veé.curl V + Vé.[Vp + curl F].| (5) 


4. Taking, for instance, 4 = lg@, where @ is the potential temperature, 
one has: \/7%.(Vp x Ya) =0. In addition, one makes the assumptions 
that frictional and non-adiabatic effects may be ignored. With these simplifi- 


cations, the equation (5) is found to express the conservation (FG =0) 


of the potential vorticity Z = as curl, V, where s = | = designates the 
static stability. 


Neglecting the slope of the isentropic surfaces, the n axis may be identi- 
fied with the ascending vertical z and accepting further the hydrostatic 
hypothesis, Z assumes the approximate form: Z = — & ¢, where 
g is the acceleration of gravity and where ¢ designates the vertical compo- 
nent of the absolute vorticity; hence, Rossby’s vorticity theorem®: in the 
absence of friction, the ratio €¢: Ap, where Ap is the pressure difference between 
two isentropic surfaces, is a conservative quantity, provided the motion is 
adiabatic. Consequently, vertical stretching of an isentropic layer can 
increase its absolute vorticity ¢. The depth of an isentropic layer can vary 
considerably in space; for instance, it could be that the weight Ap of the 
“‘deep ” regions be twice as large (or even more®) as the weight of the 
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“shallow” regions. The air moving adiabatically remains in the isen- 
tropic layer and consequently, the air columns will be shrinked vertically 
in the “ shallow” regions and stretched vertically in the “‘ deep ” regions. 
The potential vorticity ¢: Ap being advected invariantly, the absolute vorticity 
t will increase in the “* deep ” regions and decrease in the “* shallow ” regions. 
If, ab initio, the value of the absolute vorticity is normal in the “* shallow ” 
regions and if potential vorticity is advected with the air from the “ shallow ” 
to the ‘“‘ deep” regions (along the isentropic gradient — V/s or approxi- 
mately the horizontal gradient — 7,5 of the static stability s), the absolute 
vorticity can reach in the “deep” regions a value twice as large (or even 
more) as the normal value.* 


5. The conservation law of the potential vorticity Z can be written out 
as follows: 


1 d(curl,V) , lds = 


where, in virtue of the continuity equation, 
| da 
a at 


For large-scale motion in the atmosphere, the following simplifications are 
justified in most cases?: 


= div V = div u. 


(a) curl, V = curl, V = 2w, + 
1 
{(b) S= (7) 


divue 


>0 or <0 according to w +0 or<0, 
where, (u. v, w) are the cartesian components of the air velocity u with res- 
pect to the earth (wind components) in a system of local cartesian co-ordi- 
nates (x, y, z), the axis x being drawn to the East, the axis y to the North 
and the axis z to the zenith of the point considered. 


It should be noted that the vertical component ¢ of the absolute vorti- 
city is practically always positive (cyclonic) and that changes of ¢ with lati- 
tude are roughly as those of the Coriolis parameter. The yearly mean value 
of the absolute vorticity ¢ varies from 5-10~ sec! in the low latitudes to 
15-10 sec.-! in high latitudes,® while the mean value of the static stability 
is very nearly equal to 10-* m.— in the troposphere and varies between 3-10-® 


and 5-10-° in the stratosphere.’ 
Ad 
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From equations (6) and simplifications (7) we can deduce the following 
conclusions : 


(A) The absolute vorticity ¢ of a sinking air parcel increases, per contra 
the absolute vorticity of a lifted air parcel decreases. 


(B) During the formation of a new tropopause at lower levels, in a 
large developing cyclonic vortex (with respect to the earth),* the air parcels 


of the upper layers sink (5 <, 0) and the cyclogenetic deformation 
field acts in such a way that the static stability of the sinked parcels decreases 
considerably ¢ < 0), entailing the creation of an important amount of 


cyclonic vorticity. During this process of transformation of “‘ stratospheric 
air” into “ tropospheric air” through lateral outflow, the decrease of static 
stability may be as large as twice (or even more) the mean value of the static 
stability in the troposphere, doubling (at least), in virtue of (6), the value of | 
the absolute vorticity curl, V ~ ¢. 


(C) During the formation of a new tropopause at higher levels, in a 
large developing anticyclonic vortex (with respect to the earth),' the air parcels 


of the upper layers are lifted (G> a > 0) and the anticyclogenetic deformation 


field acts in such a way that the static stability of the lifted parcels 


increases considerably e > 0), entailing the creation of anticyclonic vorti- 


city. 

6. In order to estimate the maximum order of magnitude of the increase 
of absolute vorticity ¢, we shall consider an extreme and oversimplified 
example, namely, the advection. of an air parcel at 8 km. above sea level, 
from latitude 70° to latitude 40° N. At this level, the absolute vorticity ¢ 
is of the order of 14.10-° sec.-! at 70° N, while the static stability s varies 
from 3.10 at 70°N. to at 40° N.° Assuming adiabaticity 
and neglecting the horizontal gradient of the specific volume a, the product 
s¢ remains constant during the advective process, so that the absolute 
vorticity of the parcel at 40° N. reaches the value 42.10~ sec.-!, that is to say 
4 to 5 times the Coriolis parameter at this latitude. This crude example 
shows that very large values of the absolute vorticity (twice the Coriolis 
ee or even more) can occur in the atmosphere. 


7. Substituting now the identity - = 0 in the general vorticity 


equations (1), these equations assume Be form of an equation of balance’: 


‘ 
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1 2 + 


aiv — ¥ curl V) = + (8 a) 


which is appropriate to the study of local vorticity changes. 


The k-component of the absolute vorticity curl V in an unit volume 
(fixed with respect to the relative frame) can change in consequence of: 


(a) a vorticity flux & vy — v¥ curl V across the boundary of the volume; 


(6) a production within the volume at the rate s* + f*, determined by 
the isosteric-isobaric solenoids and the vortex tubes of the frictional force F. 


The flux of the k-component of curl V is defined by the contravariant 
components —  — &, the k-component of which is equal to zero 
identically (C =0); hence, the flux of the k-component of the absolute 
vorticity has no component normal to the co-ordinate surface x* = constant. 
This vorticity flux can be decomposed ina convective flux (é*v’) and a.non- 
convective. flux (— v*é’) due to the air flow across the surface x* = constant. 
Outside this surface, the two partial fluxes compensate each other exactly, 
the resultant flux normal to the surface x* = constant being equal to zero. 
Moreover, when the surfaces x* = constant are substantial surfaces (v = 0), 
the k-component of the convective flux and the three components of the 
non-convective flux vanish. The convective flux results from the transport 
of the instantaneous absolute rotation about the normal to the surfaces 
x= constant by the air flow in these surfaces. The non-convective flux is 


due’to the instantaneous absolute rotation of the air satel across’ the same 
surfaces. 


8. Taking again the equiscalar surfaces 4 = constant with the assump- 
tions (3), equations (8) and (8a) become respectively’: 

1 2, 
and 

+ — Xp) = + So (9) 


where y, = 711/22 — (7x2) is the discriminant of the metric form in the 
surfaces y = constant, div, the corresponding two-dimensional divergence, 
vy and X, the tangential components at the surface ; =- constant, at the 
point and at the time considered, respectively of y and X = curl V. 
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When the ¥-surfaces are isentropic surfaces, we have evidently s, =0; 
moreover, above the friction layer, we can neglect f,. Hence, the equation 
of balance of the absolute vorticity X,, normal to the isentropic surfaces: 


Mt (Vy,Xn) + divs — v,, curlg V) = 0. (9 6) 


Of course this equation is valid only in the free atmosphere. 


The local increase, per unit area in the isentropic surfaces, of the absolute 
vorticity X,, normal to these surfaces, is equal to the convergence of the isen- 
tropic vorticity flux X,Nvg — v, curls V. The convective part of the flux is 
due to the transport by the isentropic flow (adiabatic flow) of the instanta- 
neous absolute rotation of the air around the normal to the isentropic surfaces 
and the non-convective part, to the instantaneous absolute rotation of the 
air crossing the isentropic surfaces (non-adiabatic flow) around the axis of 
the isentropic component of the absolute vorticity The non-convective 
vorticity flux results from non-conservative processes through which abso- 
lute vorticity normal to the isentropic surfaces can be generated.* It should 


be noted that © <0 or v, <0, when y, — 0 for large scale air motion in 
atmospheric disturbances. 


9. The equation of balance of the vertical component of the absolute 
vorticity curl V can be deduced immediately from (8) or (8 a), using spherical 
polar co-ordinates (A, 4, 7). Let the longitude A be denoted by x, the lati- 
tude ¢ by x* and the radius r by x*. The coefficients of the metric form in 
this system of co-ordinates are y,,; = r? cos? 4, yoo = r? and y33 = 1. Intro- 
ducing now in (8), for A = 3, the zonal, meridional and vertical components 
of the wind and of the absolute vorticity, instead of the corresponding contra- 
variant components, one obtains the equation of balance of the absolute 
vorticity 

— we) + 2 Ir cos $ (Lv — = 5, + (10) 
ox r COS dy 


— -+ div, (<v, — w curl, V) = s, +f, (10 a) 


ag 
dt 


u=rcos¢a, y = rd, 


ow ou 


or 
where® 
| 
wu 
r 
¢ = 2 sin 
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= F,ig¢ 


with " = the subindice / referring to the 
horizontal component of a vector or to an operator in the geopotential sur- 
faces (assumed here to be spherical). At the synoptic scale it is justified 
to assume that the horizontal baroclinicity is negligible (s, = 0) and above 
the boundary layer that the frictional term is inefficient ( f, ~ 0); hence, 


the approximate equation of balance 


= 


+ div, —w curl, V) = (10’ a) 


is valid in most cases. Consequently, the local changes of the vertical compo- 
nent ¢ of the absolute vorticity curl V (the so-called absolute vorticity ¢) are 
due chiefly to the convergence or divergence of the horizontal vorticity flux 


C, = tv, — w curl, V, 


the convective part of which results from the horizontal transport fv, by 
the horizontal air motion v, of the instantaneous absolute rotation ¢ around 
the vertical axis z and the non-convective part—w curl, V, from the instan- 
taneous absolute rotation of the air moving vertically in the plane normal 
to the horizontal component curl, V of the absolute vorticity. 


For the zonal (C,), the meridional (C,) and the vertical (C,) compo- 
nents of the vorticity flux, we have the following pingygpio 


= Wu — wi ~ u {Qu sin + (3° [we], 


C, = — wy = v{2u sin + > 


C.= 0. 


It should be emphasized that the terms in parenthesis cannot be deleted 
in the jet-stream region (certainly not in the vicinity of the inflection 
points of the jet-stream profile) neither the terms in brackets in those regions 


where simultaneously the vertical baroclinicity (=. sz) is very high and 


the vertical motion exceptionally important. Hitherto this last part of the 
vorticity flux has been generally neglected, 


211 
); 
and 
5) 
te 
is 
ve 
Ds 
ld 
te 
al 
i- 
in 
ts 
a- 
te 
2) 


212 


Charney, J. G. 

Ertel, H. 

Newton, C. W., 
Phillips, N. A. 

_ Carson, J. E. and 
Bradbury, D. L. 

Palmén, E. 

Petterssen, S. 

Rossby, C. G. 


Van Mieghem, J. 


—— et Vandenplas, A. 


JACQUES-M. VAN MIEGHEM 


REFERENCES 
‘On the scale of atmospheric motion,” Geof. Publ., Oslo, 1948, 
17, No. 2, 1-17. 
“‘Ein neuer hydrodynamischer Wirbelsatz,” Meteor. Zeits., 
Braunschweig, 1942, Bd. 59, 277-81. 


‘*Structure of shear lines near the tropopause in summer,” Tellus, 
Stockholm, 1951, 3, No.3, 154-71. 


“*Aerologische Untersuchungen der atmosphirischen Stéringen,” 
Soc. Scient. Fennica, Comm. Phys-Math., Helsinki, 1933, 7, 
6, 1-65. 

“Some aspects of the general circulation of the atmosphere,” 
Centenary Proc. Roy. Met. Soc., London, 1950, 120-55. 

“Planetary flow patterns in the atmosphere,” Quart. Journ. 
Roy. Met. Soc., Supplement, Toronto, 1940, 66, 68-87. 

‘*Le bilan de la rotationnelle absolue dans l’atmosphére,” Tellus, 
Stockholm, 1951, 3, No. 3, 297-300. 

‘*Les équations de la dynamique atmosphérique en co-ordonnée, 
généralisées,’” Mém. de I’Inst. Roy. Météor., de Belgiques 
Uccle, 1950, 41, 1-57. 


i. 
4 
&e 
6. 
— 
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SON ROLE DANS LA FORMATION DES CYCLONES 
PAR P. QUENEY 

(Professeur @ I’ Université de Paris, France) 


Received December 26, 1952 


Divers auteurs ont repris depuis quelques années les anciennes études 
théoriques de RAYLEIGH sur les perturbations des écoulements barotropes 
continus en vue de les appliquer aux ondes du jet-stream et aux ondes de 
cyclone, ou au probléme de l’origine de la turbulence. Les résultats obtenus 
par RossBy, aujourd’hui classiques, ont prouvé l’intérét qu’il y avait pour 
la météorologie 4 développer de telles études, mais jusqu’ici la nature réelle 
des ondes de cyclone ne semble pas avoir été bien éclaircie. Le but du 
présent article est de décrire d’abord une propriété dynamique particuliére 
de certains écoulements barotropes, pour laquelle nous proposons le: nom 
de “ résonance interne’, et ensuite d’appliquer la théorie de ce phénoméne 
aux ondes de cyclone. 


1. ONDES. ELEMENTAIRES D’UN SYSTEME DE DEUX COURANTS DE COUETTE 


Considérons dans le plan des x y un mouvement barotrope et sans 
divergence d’un fluide incompressible non visqueux. Si YW désigne sa 
fonction de courant, les composantes U, V de sa vitesse et son tourbillon Q 
sont donnés par 


et il résulte des hypothéses faites que Q est une grandeur conservative. 


Supposons maintenant le mouvement constitué par un écoulement de 
base paralléle 4 l’axe des x, déformé par une perturbation ondulatoire de 
petite amplitude et ne dépendant de la co-ordonnée x et du temps ¢ que 
par un facteur sinusoidal en x et f, dont nous désignerons la longueur d’onde 
par A, le nombre d’onde par k (k = 2m/A) et la vitesse de phase par c. 
Dans l’écoulement de base on a 


V=0; U=U(y); Q=Q(y) = — U'(y) (2) 
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et si ¥, u, v, g, désignent les déformations de ¥, U, V, Q respectivement, dues 
a la perturbation, on a 


(3) 


En écrivant ensuite q ue le tourbillon total Q + q est conservatif, on obtient 
pour u l’équation différentielle suivante si on assimile cette grandeur 4 un 
infiniment petit: 


(U—c) y+ (4) 
ot Q’ est la dérivée Q’ (y). 


Si en particulier le profil U = U (y) est une ligne brisée continue formée 
de segments de droite, "écoulement est un syst¢me de courants de Couette 
sans discontinuité de vitesse, et dans ce cas Q est une constante dans chaque 
courant mais varie de fagon discontinue sur leurs lignes de séparation 
(frontiéres). Dans chaque courant et partout ol u est continu l’équation (4) 
se réduit alors a 


oug=0 (5) 


tandis que sur une frontiére ou sur une ligne de discontinuité de u (ligne de 
glissement) on a, par intégration de (4) a travers cette ligne, sur laquelle ¢ 
doit étre continu, 


(U — c) Au= ¥AQ (6) 


en désignant par Au et AQ les accroissements discontinus de u et QO dans 
le sens des y croissants. 


Cette relation montre que toute frontiére est nécessairement une ligne 
de glissement, ol g a par suite une amplitude infinie, mais dont l’intensité 
tourbillonnaire par unité de longueur, définie par J = / q dy, a une valeur 
finie égale a — Au. La méme formule (6) montre aussi que si on a une 
ligne de glissement supplémentaire (autre qu’une frontiére) on a nécessaire- 
ment c = U sur cette ligne puisque AQ n’est pas nul: donc la perturbation 
se déplace avec la vitesse que posséde le mouvement de base sur la ligne de 
glissement supplémentaire, et par suite le nombre total de ces lignes est au 
plus égal au nombre des courants de Couette. 


En raison de ce résultat nous appellerons “ ondes convectives”’ les 
perturbations qui ont une ou plusieurs lignes de glissement supplémentaires, 
tandis que celles qui n’ont que les fronti¢res comme lignes de glissement 


I 
I 
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seront appelées ondes fondamentales”. L’ensemble des ondes fonda- 
mentales et des ondes convectives ayant une seule ligne de glissement supplé- 
mentaire forme ce que nous appellerons les ‘‘ ondes élémentaires’’, en 
raison du fait que toute perturbation du systéme peut étre obtanue en super- 
posant un nombre fini ou infini de telles ondes. 


Dans la suite nous ne considérerons que le cas le plus simple, celui d’un 
systeme de deux courants de Couette illimités, séparés par la ligne y = 0; 
si Q, et Q, sont les deux valeurs constantes de Q, et U, la vitesse sur la 
frontiére, on a alors 


U= U, pour y > 0 


U=U,— Q:y, poury< 0 (7) 
et nous poserons pour simplifier l’écriture 


o = (8) 


Nous supposerons d’ailleurs que o et Q, sont positifs moyennant des sens 
positifs convenables sur les axes co-ordonnés. Dans le cas ot Q, est 
négatif l’écoulement a un maximum de vitesse pour y = 0. 


(a) Ondes fondamentales 
Si on se donne l’intensité tourbillonnaire de la frontiére y = 0, soit 
Ig = Ag Cos k (x — ct — Xp) (9) 


A, et sont deux constantes, la répartition de se trouve complétement 
déterminée par le fait que cette grandeur doit étre solution de 1’équation 
de Laplace dans chaque courant de Couette, continue sur la frontiére, et 
nulle 4 Vinfini; on a ainsi 


st etm — cos k (x — ct — x») (10) 


Donc toute onde fondamentale est une onde de surface simple ayant 
la frontiére pour axe. En écrivant ensuite que la condition (6) est satisfaite 
sur la frontiére on obtient l’équation de dispersion des ondes fondamentales, 
a savoir 

@ 
La déformation transversale 9 de la frontiére y = 0 est donnée par 


no = — = — cos k (x ef — x4) (12) 
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et plus généralement la déformation » de la ligne fluide paralléle a l’écoule- 
ment et d’ ordonnée » a pour expression 


A 


Cette formule montre que, alors que #, u et v ont leur amplitude la plus 
grande sur la frontiére y = 0, la variation de 7 en fonction de y est toute 
différente, et en particulier l’amplitude de » devient infinie lorsque U =, 
ce qui se produit pour lunique valeur positive y = o/kQ, si Q, est positif, 
et en outre pour la valeur négative y = o/kQ, si Q, est négatif. Le fait de 


trouver une valeur infinie pour 7 entraine que les lignes fluides voisines sont 
infiniment resserrées. 


Inversement 4 toute valeur positive de y il correspond toujours une 
longueur d’onde et une seule telle que l’on ait c = U(y) pour les ondes 
fondamentales ayant cette longueur d’onde, et if en est de méme si y est 
négatif dans le cas ot Q, est également négatif: pour une raison qui apparaitra 
plus loin nous dirons que cette longueur d’onde particuliére est la “* longueur 
d’onde de résonance ” relative 4 la valeur considérée de y. 


(b) Ondes convectives 


Nous considérons ici une onde ayant, en plus de la ligne de glissement 
y = 0, une seule ligne de glissement supplémentaire y=h. Si Jp et J, 
désignent les intensités tourbillonnaires respectives de ces lignes on obtient 
alors, en expérimant que satisfait a l’équation de Laplace en 
dehors des lignes de glissement, est continu sur elles et s’annule a l’infini, 


Donc toute onde convective élémentaire est la superposition de deux ondes 
de surface simples ayant les lignes de glissement pour axes. II! nous reste 
maintenant a exprimer que la condition (6) est satisfaite sur chaque ligne 
de glissement. Pour qu’elle le soit sur y=A/ il faut que l’on ait c = U, 
quelle que soit la longueur d’onde, ainsi que nous l’avons déja dit; et si on 
écrit que la condition est satisfaite pour y = 0 on obtient une relation entre 
I, et I, que l'on peut mettre sous la forme 


I, ~ Up) 


ol cg désigne la vitesse de phase des ondes fondamentales de méme longu- 
eur d’onde, c’est-a-dire la valeur de ¢ donnée par la formule (11). Il en 
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résulte que l’onde est complétement déterminée si on se donne par exemple 
lexpression de /,,, laquelle doit étre de la forme 

I, = A, cos k (x — U,t — x4) (16) 
ou A, et x, sont deux constantes. 


La déformation 7, de la frontiére est ici donnée par la formule 


La formule (15) montre que, lorsque /, est donné, l’amplitude de Jy 
tend vers l’infini lorsque cy tend vers U,, ce qui se produit lorsque A tend 
vers la longueur d’onde de résonance A, relative a la ligne y = h (si cette 


longueur d’onde existe): c’est ce résultat qui est 4 la base du phénoméne. de 
résonance dont nous allons nous occuper maintenant. 


2, RESONANCE DANS UN SYSTEMR DE DEUX COURANTS DE COUBTTE 
(a) Résonance aigiie 
L’onde fondamentale définie par (9) est la perturbation qui résulte d’une 


répartition initiale de q consistant en une concentration infinie sur l’unique 
ligne y = 0, avec l’intensité tourbillonnaire 


la valeur de q étant nulle partout ailleurs. 
De méme l’onde convective définie par (16) est la perturbation qui 


résulte d’une répartition initiale de g consistant en une concentration infinie 


sur les deux lignes y = h et y = 0, avec les intensités tourbillonnaires res- 
pectives 


i, = A, cos k (x (19) 


De la on déduit que la perturbation due a une concentration initiale infinie 
de q sur l’unique ligne y = h, avec Vintensité tourbillonnaire (19), s’obtient 
en superposant les deux ondes précédentes 4 condition de prendre i, = — i,', 
c’est-a-dire 


Xo — A; (21) 


A un instant quelconque / cette derniére perturbation comporte donc les 
deux lignes de glissement y = h et y = 0, avec les intensités tourbillonnaires 


| ig = Ag cos k (x — XX) (18) 


218 


P. QUENEY 


respectives J, donnée par (16), et J, donnée par 


I, = A, Kc — U,) [cos k (x — U,t — x,) — cos k (x — egt—x;,)] (22) 


et on a encore n, = — /,/2c. 


On a ainsi au total la superposition de trois ondes de surface simples, 
et on voit que celles qui ont la ligne y = 0 pour axe donnent lieu 4 un phéno- 
méne de battements si cy est différent de U,. En mettant J, sous la forme 


n= ~ sin 


(23) 
on constate en outre que la période et l’amplitude des battements tendent 
toutes deux vers l’infini lorsque cy tend vers U;,, c’est-a-dire lorsque A tend 


vers A, (ou lorsque k tend vers k,= 27/A,), et en prenant la limite de J, on 
obtient, lorsque A = A,, 


I, = — A,oe~*n'*' t sin k, (x — U,t — x,) (24) 


On a d’ailleurs k, | | = 0/Q, si h est positif, et k,|h | | si h est 
négatif, dans le cas ot Q, est aussi négatif. 


Donc pour la longueur d’onde de résonance A, les deux ondes de surface 
relatives a la ligne y = 0 se fondent en une seule onde de surface dont l’ampli- 
tude croit indéfiniment avec le temps (proportionnellement a 1), de sorte 
que c’est une onde de surface instable, et cela ne se produit que pour cette 
longueur d’onde: il s’agit évidemment d’un phénoméne de résonance aigiie, 
ce qui justifie le terme “ longueur d’onde de résonance”’. L’onde instable 
elle-méme est ce que nous appellerons: |’ “‘ onde résonance”’ relative a la 
déformation initiale définie par (19). Au contraire la troisitme onde de 
surface, ayant la ligne y = A pour axe, a une amplitude constante, et on voit 
qu’elle est négligeable en comparaison de l’onde de résonance lorsqu’on a 

t > ty avec fy = (25) 


Ajoutons que. pour l’onde de résonance, la déformation transversale de la 
ligne fluide y =A est infinie, tout comme dans le cas d’une onde fonda- 
mentale. 


On doit admettre que l’énergie de l’onde de résonance est, comme dans 
toute perturbation instable, constamment empruntée 4 l’énergie mécanique 
du mouvement de base et non fournie par la perturbation initialement intro- 
duite sur la ligne y = h; mais c’est tout de méme cette perturbation initiale 


il 
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qui détermine le développement de l’onde de résonance et dont l’amplitude 
fixe la vitesse de croissance de cette derniére. 


(b) Résonance amortie 


(1°) La méthode de Fourier permet d’obtenir plus généralement la 
perturbation due 4 une concentration initiale infinie de Q sur la ligne y =/ 
avec l’intensité tourbillonnaire 


i, = f(x) (26) 


oi f(x) est une fonction quelconque de x: il suffit de décomposer cette 
fonction en ses éléments sinusoidaux, puis de superposer les perturbations 
élémentaires correspondantes. Or il est évident que, comme conséquence 
du phénoméne de résonance précédemment décrit, les perturbations élé- 
mentaires dont la longueur d’onde est voisine de A, finissent par prédominer 
sur toutes les autres, et en conséquence la perturbation prend un caractére 
ondulatoire de plus en plus marqué: lorsque la condition (25) est remplie 
cette perturbation consiste pratiquement en une seule onde de résonance 
qui est une onde de surface ayant la frontiére » = 0 pour axe, A, pour 
longueur d’onde, et U, pour vitesse de phase; mais comme ici il y a inter- 
férence entre toutes les perturbations composantes de longueurs d’onde 
voisines de A,, amplitude de cette onde de résonance ne peut plus croitre 
indéfiniment avec le temps: elle tend vers une limite finie au bout d’un temps 
comparable 4 f, aprés avoir éprouvé une série d’oscillations d’amplitude 
décroissante. De méme la déformation », de la ligne fluide y =A reste 
finie, mais les lignes fluides voisines sont fortement tassées les unes contre 
les autres. 


A titre d’exemple si on choisit pour f(x) la fonction 


ou S et a sont deux constantes, et qui est représentée graphiquement par une 
courbe en cloche symétrique d’aire S, la valeur asymptotique de J, pour 
U;t U, est 


Ip — 2S ky gin k, (x — (28) 


et en faisant a = 0 on obtient |’effet d’une concentration initiale ponctuelle 
de g (concentration infinie au point de coordonnées x = 9, y=h). On 
obtient ainsi un train d’ondes dont la longaeur croit proportionnellement 
au temps. Cet exemple suffit 4 montrer comment une perturbation initiale 
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qui n’a absolument rien de périodique crée cependent une onde sinusoidale 
grace au phéno méne de résonance. 


(2°) Une concentration infinie de g sur une ligne y = / n’est pas physique- 
ment concevable, mais si 4 |’instant initial la déformation du tourbillon est 
localisée dans une bande suffisamment étroite paralléle 4 I’écoulement il 
est clair que l’on aura encore un effet de résonance amorttie. 


La perturbation peut étre obtenue en décomposant la zone perturbée 
en bandes infiniment étroites paralléles 4 l’axe des x et que !’on assimile a 
des lignes sans épaisseur, puis en superposant toutes les perturbations 
partielles qu’elles produisent, ce qui conduit encore 4 une intégrale, mais 
ou la variable d’intégration est cette fois y au lieu de k. On trouve ainsi 
que !’effet de résonance demeure important tant que la largeur de la zone 
perturbée est petite en comparaison de sa distance moyenne | h/ | a la frontiére 
y = 0, la perturbation limite (onde de résonance) étant encore un train 
d’ondes d’amplitude bien déterminée lorsque la condition (25) est satis- 
faite; sil n’en est pas ainsi, la perturbation n’arrive jamais 4 prendre le 
caractére ondulatoire, 4 moins que la déformation initiale ne soit déja sinu- 
soidale et de longueur d’onde voisine de A,. 


3. 


APPLICATION AUX ONDES DE CYCLONE 


Le probleme des ondes élémentaires d’un systeme de courants de Couette 
a son équivalent mathématique exact lorsqu’on considére un mouvement 
barotrope et sans divergence qui s’effectue, non plus sur un plan fixe, mais 
sur une surface de niveau supposée également fixe et en outre de révolution 
autour de la ligne des pdles: léquivalent d’un courant de Couette est alors 
un courant zonal (paralléle a l’équateur) dans lequel la composante verticale 
du tourbillon absolu est indépendante de la latitude, ce que nous pouvons 
appeler un “courant de Rossby”, et l’équivalent d’une onde sinusoidale 
én x et f est une onde sinusoidale en L et t, ol. L désigne la longitude. II 
en résulte que le phénoméne de résonance que nous avons décrit existe aussi 
dans le cas d’un systeme de deux courants de Rossby séparés par une bast 
de discontinunuité du tourbillon vertical absolu. 


Or l’acuité du maximum de vitesse du jet-stream aux niveaux voisins 
de celui de la tropopause permet d’assimiler cette partie centrale de la circula- 
tion générale 4 un tel systeme, du moins tant que l’on cherche seulement, 
comme nous le faisons ici, 4 élucider la cause physique des ondes de cyclone; 
d’ailleurs il se trouve justement que dans beaucoup de cas le jet-stream n'est 
pas loin d’étre un courant de Rossby sur le versant polaire de son maximum 
de vitesse, ce qui renforce encore la légitimité de cette approximation. 


| | 
| 
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D’autre part tous les mouvements atmosphériques de grande échelle sont 
approximativement horizontaux et barotropes aux niveaux voisins de la 
tropopause, et comme c’est la que les ondes de cyclone ont normalement 
leur plus grande intensité et leur caractere ondulatoire le plus net, nous 
avons tout naturellement été conduit 4 la conclusion suivante: 


Une onde de cyclone ne serait pas autre chose que l’onde de résonance produite 
dans le jet-stream par une perturbation locale du tourbillon vertical. 


D’ aprés la théorie de onde de résonance, cette perturbation doit étre 
introduite dans une zone suffisamment éloignée de l’axe du jet-stream en 
comparaison de sa largeur, et il semble que l’on puisse admettre une distance 
égale 4 la moitié de cette largeur comme minimum extréme: si on adopte 
une largeur minimum de 500 km, cela fixe a la fois le maximum de vitesse 
des ondes de cyclone et le minimum de leur longueur d’onde, soient respec- 
tivement 30 m/sec et 2000km environ pour un jet-stream moyen ayant 
40 m/sec de vitesse maximum et pour lequel o = 3, 5x 10~ sec". 


D’un autre cété la durée d’établissement de l’onde résonance, qui est 
de l’ordre de la quantité f de la formule (25), ne doit pas excéder 2 a 3 jours. 
Or, contrairement au cas d’un courant de Couette, le gradient méridien de 
la vitesse d’un courant de Rossby varie avec la latitude, et pour un jet-stream 
assimilé 4 un systeme de deux courants de Rossby ce gradient décroit vers 
le Nord et vers le Sud a partir de la frontiére des deux courants. De la il 
résulte que le produit k, | h |, au lieu d’étre constant dans chaque courant, 
croit avec | / |, donc il en est de méme de ¢,: cela fixe 4 la fois le minimum 
de vitesse des ondes de cyclone et le maximum de leur longueur d’onde, 
soient respectivement 15 m/sec et 4000 km environ pour un jet-stream moyen. 


Ces considérations, bien que tres grossi¢res, montrent cependant que 
les ondes de résonance du jet-stream sont bien quantitativement compa- 
tables aux ondes de cyclone réelles, ce dont il était évidemment indispensable 
de s’assurer: Ajoutons que le tassement des lignes fluides zonales, donné 
par la théorie de l’onde de résonance, dans la zone initialement perturbée, 
pourrait peut-étre aussi fournir une explication du phénomene de fronto- 
génése qui est toujours associé 4 la cyclogénése, et qui joue un réle si impor- 
tant dans la formation des systemes nuageux et pluvieux accompagnant les 
ondes de cyclone. 


Enfin il nous reste 4 voir quel peut étre le processus physique qui conduit 
a la déformation du tourbillon vertical déterminant le développement de 
onde de résonance. A ce sujet il est tout naturel de faire appel aux phéno- 
ménes ayant leur si¢ge dans la troposphére inférieure, ot la baroclinité de 
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lair permet un transport vertical du tourbillon: or c’est bien en effet au- 
dessus des régions 4 forte baroclinité (régions de forts contrastes horizontaux 
de température) que les cyclones extra-tropicaux apparaissent, conformément 
au pricipe méme de la théorie norvégienne. D/ailleurs, une fois que l’onde 
de cyclone est devenue prépondérante, elle est nécessairement accompagnée 
d’une déformation ondulatoire de méme longueur d’onde dans la troposphére 
inférieure, et celle-ci 4 son tour doit logiquement produire dans le jet-stream 
une déformation du tourbillon ayant aussi cette méme longueur d’onde; 
or cela ne peut que renforcer l’onde de résonance: en d’autres termes cette 
derniére une fois amorcée peut s’entretenir d’elle-méme, de la méme facon 
que les oscillations d’un poste d’émission radioélectrique par exemple. 


Ainsi se trouveraient conciliées les théories thermiques et dynamiques 
de la cyclogénése: il y aurait couplage, par l’intermédiaire de la troposphére 
moyenne, des effets thermiques prédominant au-dessous avec les effets 
dynamiques prédominant au niveau de la tropopause, et c’est de ce couplage 
que résulteraient en définitive les ondes de cyclone. Evidemment une 
théorie compléte devrait tenir compte du fait que le couplage existe en 
réalité 4 tous les niveaux, mais il ne semble pas que cela puisse diminuer 
sensiblement le rdle de l’effet de résonance que nous venons de mettre en 


lumiére. 


A NEW RADIO THEODOLITE 
Preliminary Communication 
By VILHO VAISALA 
Received January 2, 1953 


Tue radio theodolite here described has been developed especially for 
tracking the ordinary signals of the Finnish (Vaisala-7) radiosonde. Because 
of the rather long wave (about 12-5 m.) of this radiosonde the usual methods 
employing movable antenna systems are incapable of giving a satisfactory 
accuracy of direction. The dimensions of the antenna system must be 
enlarged to several times the tracked wave-length and therefore only fixed 
antenna systems can give better accuracy. 


The principle of the method here described is the following: when the 
radiosonde (S, Fig. 1) is relatively far away from two fixed antennas (A and 
B) the wave surface of the signal sent by the radiosonde can be considered 


Fie. 1 


a plane. After reaching the one antenna (B) the wave surface (BC) must 
still proceed the way CA in order to reach the other antenna (A). This 
distance CA = 4 has to be measured. When that has been done we obtain 
the direction angle CAB = a between the incoming radio ray (SA || SB) and 
the direction AB. Let the antenna distance AB=a. Then 


cosa = 4 (1) 
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To get the direction of the ray we must measure two direction angles in 
respect of two fixed directions. For this purpose a second antenna pair 
A’B’ = a’ is provided and we obtain the direction angle 8 in respect of said 
second antenna pair from 


cos = (2) 


In practice a’ =a and A’B’ is perpendicular to AB, both antenna pairs 
being situated in one and the same horizontal plane. 


To measure 4, the wave-length A of the radio ray is used as unit. N, 
may be the phase difference of the wave surface in the positions BC and DA, 
i.e., between B and A, and N, the phase difference between B’ and 4’, 
Then 

4 = N,A, 4’ =N,A (3) 
and we get 


A A 
cosa = cosB = (4) 


The principal feature of our radio theodolite is that the phase differences 
N, and N, are directly recorded on the recording sheet as functions of time. 
By means of cos a and cos f the elevation angle A of the direction to the 
radiosonde as well as the azimuth angle A relative of the x-axis can be 
reckoned : 


cos a = cosh cos A, 


cos 8 = cosh sin A, (5) 
cosh = costa + cost = 


cos8 WN 
(6) 
These expressions show that A is independent of wave-length (or fre- 
quency) but A dependent upon it. To obtain the direction of the radio ray 
we must alongside N, and N, record the wave-length A (or frequency v) as 
the wave-length varies in some degree during the sounding. Thus the 
record consists of curves for N,, N, and A (or ¥). ‘ 


For practical wind computation we develope here a numerical method 
suitable especially for calculating machines. First we calculate the co-ordi- 
nates of the radiosonde supposing the height H as known from the sounding 
of PTU. Then the horizontal co-ordinates of the radiosonde are 
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sin A 


Yon 


sin A 
We get with (6): 


sinth = ($) —cost) = (5) es 


sinh = (3) — (8) 


For practical reasons we have decided to use 


a= 120m. 
Further we write 
A, => 12m. (10) 
Then 
a_anr v 
i. = 10% = 10 Ye (11) 


Wherein v and vg are frequences of A and Ay. vy= 25 mc/s. Because the 
wave-length A sent by the greater constant condenser of our radiosonde is 
about 12-5m. and always > 12m., 


A> A,and v< yor 


(10 < 100. (12) 


The-expression (8) can be written as 


«sin h = 100 — + AN) (13) 
wherein 
= 100 [1—(2 4) 
Denoting 
C= N,2+N,?+ 4N (15) 
100 
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we get from (7) 


Y = ) (17) 


Addendum: 
Here H, X and Y are in metres. 


Now the wind components are obtained as usually from — 


(18) 


wherein 4t is the time interval (one minute of perhaps two or three minutes). 


If we compute X and Y at all minutes it is suitable, especially when the 
radiosonde goes farther away, to smooth the results, e.g., in following 
manner: By means of the wind components (v,, v,) the hodograph of the 
wind vector is designed. A mean vector always for three successive points 
is determined by means of gravity centre method. These gravity centres 


give the smoothed hodograph wherefrom all requisite wind vectors can be 
read. 


The practical steps in calculating the co-ordinates (17) are as follows: 


1. The PTU-sounding gives H (in kilometres) or it is computed by 
means of the approximately known ascending rate of the balloon. 


2. N,, N, and 4N are read from the radio theodolite record. 


3. The quantity C (15) is calculated by means of a calculating machine. 
Here it is important the machine calculates the quadrates directly. 


4. S§(16) is taken from a table with C as argument. 
5. The product HS is calculated. 


6. HS is multiplied successively with N, and N, (17). 


X and Y are then obtained in dekametres whole dekametres give a suitable 
accuracy. 


In this way the computation sheet contains the following columns: 
|Minutss| H | N, | N, | 4N | X¥ | Y | 4X¥ | AY | 


Vy At > |_| 
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My opinion is that this computation method employing a calculating 
machine is better than any graphical method. 
The following table is an abstract of the original S-table. 


S= Vi00 —¢ 


| § 


10-000 | 18-257 
10-541 20-000 


11-180 | 22°36] 
11-952 
12-910 
14-142 


23-570 
25-000 
26-726 


15-81] 


31-623 


65 33-333 


16-903 


| | 
| 
| | 
| | 
14-907 28-868 
| | 
| | 
| | 


In our new (third) radio theodolite prototype* 100 mm. on the record 
corresponds to one wave-length A. As N, and N, can.vary between — 10 
and + 10 the width of the record would be 2000 mm. As the place on the 
recording sheet is only about 200mm. the record must jump one period 
back when the record is about to overpass the paper. 


Because of the periodicity of the phase angle the question arises as to 
which period is being recorded. For this purpose an other antenna system 
(search antennas) with a = A, is used and here the phase differences n, and 
n, vary only between about — 1 and + 1. The period can in this way be 
easily determined. For instance when the search antennas are switched in 
we may have n, = + 0-675; then N, = 10n, = + 6°75, i.e., the record 
is going on between + 6 and +7. Because of the very great changings 
of balloon direction the search antennas are also used to record the wind 
at the outset, after release of balloon. 


* This will be ready in the beginning of 1953, 


| | s 
92 35-355 
10 93 | 37-796 
) | 20 | 94 40-825 
’ 30 95 | 44-721 
e 40 96 | 50-000 
g | 
2 50 97 | 57-74 
8 55 98 70-71 
60 99 100-00 
100 
| 
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As to the accuracy of the radio theodolite I can not as yet give any 
exact numbers. The comparisons of the two first (yet very incomplete) 
prototypes with an optical theodolite have given in azimuth a relative accu- 


racy of about + 0°-1 and in elevation angle an accuracy + 0°-1 at 45° but 
+ 0°-5 at 12° elevation. 


It has not been the purpose to treat in this connection the radiotechnical 
side of the equipment which is designed by Y. Wilska, M. Eng. 


University of Helsinki, 
Institute of Meteorology in December 1952. 


— 


ON THE DYNAMICS OF THE INDIAN MONSOON 


By SVERRE PETTERSSEN 
(The University of Chicago) 


Received January 9, 1953 


1. ALTHOUGH it is now more than 250 years since it was suggested (Halley, 
1686) that the monsoon is caused by differential heating, very little progress 
has been made to explain the Indian monsoon, and its variations. on thermo- 
dynamica! principles. 


The difficulties in the way of a thermodynamical treatment of the details 
of the monsoon are, of course, enormous. Even if the precise distribution 
of the sources and sinks of heat over India were known, the resulting circula- 
tion would depend upon many other factors, such as the existing global 
circulation, the mechanical effects of topography, the physical properties of 
the earth’s surface, etc. 


A further difficulty arises from the circumstance that it is difficult to 
obtain a suitable scientific definition of what is meant by the term 
“monsoon”. The circulation over India, even at the height of the season, 
is not steady, and the question of the role of the perturbations in the mainte- 
nance of the mean state naturally arises. 


2. In the following we shall consider the Indian monsoon as a steady 
state which in all essentials resembles the mean state in midsummer. In 
particular, we shall assume that the potential temperature (@) and the absolute 
vorticity (Q) are steady. Thus, our basic assumptions are that 


20 _ 
0 and 0. (1) 


Since these assumptions cannot be justified by reference to dynamic 
or thermodynamic principles, the validity of these assumptions must be deter- 
mined by comparing derived results with observation. Essentially, the 
assumptions (if broadly true) mean that the mean state is largely self-sus- 
taining, and that the perturbations are not a prominent feature in the mainte- 
nance of the monsoon. 


3. The essential feature of the monsoon that we desire to consider is 
the vorticity of the motion. The mean state of the motion at low levels 
‘over India is characterised by a large eddy of cyclonic relative vorticity. The 
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question is whether it is possible on dynamic and thermodynamic principles 
to explain the maintenance of this eddy. 


In endeavouring to answer this question it is necessary to bear in mind 
that the Indian subcontinent is situated to the south of the so-called Jet 
Stream, with the result that the relative (and absolute) vorticity decreases with 
elevation. This feature is largely dependent upon the global circulation of 
the atmosphere, and cannot be regarded as part of the monsoon itself. 


Furthermore, consideration must be given to the mechanical influences 
due to topography. Even a cursory inspection of the circulation charts 
produced by Dr. Ramanathan and others reveal a pronounced influence 
due to the Western Ghats, the Burma and Assam Mountains and the 
Himalayas. 


4. It has been shown elsewhere* that the vorticity equation may be 
written in the form 


Here, Q is the component of the absolute vorticity normal to a given surface, 
V the wind vector in that surface, S the number of isosteric-isobaric solenoids 
per unit area, and F the vorticity of the frictional force in the same surface, 
while V,, is the wind component normal to the surface considered. 


Adding and subtracting the term Qodv,/dn, the vorticity equation may 
be written 


where 7-(QV) denotes the three-dimensional divergence. 


5. The effect of heat and cold sources may be expressed through the 
variations of potential temperature experienced by a moving unit of air, viz., 


dé 20 


where the second term on the right is the advection along, and the third term 
the advection through the surface considered. 


6. We shall now choose to represent the conditions in an isentropic 
surface, rather than in a level surface or an isobaric surface. On account 
of this choice 

S=0 and V-7,8 = 0. (5) 


* Petterssen, S., ‘Some Aspects of the General Circulation of the Atmosphere,” Centenary 
Proceedings, Roy. Med. Soc., 1950, pp. 120-55, 
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Introducing now assumptions (1), Eqs. (2) and (4) may be combined to give 


6 
which specifies the combined dynamic and thermodynamic control of any 
stationary circulation. 


The term on the left is a measure of the intensity of the vorticity source. 
Eq. (6) states that this intensity is equal to the vorticity of the frictional force, 
plus a component which depends upon the intensity of the heat or cold source. 


We may therefore speak about a frictional and a thermal component of the 
vorticity source. 


7. The sign of the frictional component near the earth’s surface is 
readily determined. If the frictional force is to act mainly against the 
motion, the vorticity of the frictional force must be positive or negative 
according as the vorticity of the motion is anticyclonic or cyclonic. There- 
fore, F < 0 in the lower layers of the monsoon. 


8. The sign of the thermal component depends not only upon the sign 
of the heat source but also upon whether the absolute vorticity increases or 


decreases in the direction of increasing entropy. As indicated in para. 3, 
2Q/34 < 0 over India. Since any air that arrives in India from the adjoin- 
ing seas will be heated, d@/dt > 0, and, hence, the thermal source of vorticity 
is positive. 


Part of the vorticity which steadily is being created by thermal processes 
is consumed by friction and given back to the earth; some of it may be 
exported to other regions, particularly downwind to the Bay of Bengal. The 
steady state is realized when the thermal source is balanced by the frictional 
dissipation and the export. 


If the pressure distribution is to adapt itself to the resulting vortex 
motion over India, such that some measure of geostrophic balance tends to 
be established, the result must be the creation and maintenance of an area 
of low pressure over India. 


As a result of the diurnal variation of the intensity of the heat source, 
one would expect to find a corresponding variation in pressure, and this 
variation should not be limited to the coastal strip effected by the land and 
sea breeze system. 


9. There remains now to be considered the mechanical influences due 
to the topography. If India were level with the sea and the surface physically 
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uniform, the isobars would tend to be parallel to the coast line. However, 
the main import of air into India is across the Western Ghats, with ascending 
motion on the western slope and descending motion on the eastern slope. 
The effect of this vertical motion may be seen from Eq. (3). Making use of 
Eqs. (1) and (5), Eq. (3) reduces to 


P-(QV) = F+Q 


Now, on the windward side of the range dV,,/dn < 0, and (apart from the 
frictional effect) absolute vorticity is consumed. On the lee side 0V,,/dn > 0, 
and absolute vorticity is created. Since friction must be regarded as a modi- 
fying, but not as an over-riding effect, the net result will be a ridge of anti- 
cyclonic vorticity over the Western Ghats, with a gradual change to a maximum 


of cyclonic relative vorticity further to the east, and also north of the northern 
shoulder of the Ghats. 


The effect of the Eastern Ghats in regard to Burma is similar to that 
of the Western Ghats in regard to India. The Assam Mountains and the 
Himalayas are not effectively crossed by air currents, and their effects reside 
mainly in a frictional consumption of vorticity. 


10. The foregoing discussion has dealt exclusively with the problem 
of the maintenance of the mean state at low levels over India, without appeal 
to perturbations or dynamic coupling with neighbouring circulation systems. 
Since the deductions are in close agreement with the observed mean state, 
it seems permissible to conclude that the Indian Monsoon, to a very large 
extent, is a self-sustaining system which is driven by the local heat source 
and balanced by frictional dissipation and, that perturbations from neigh- 
bouring regions are not essential for the maintenance of the mean state. 


The variability of the monsoon may, however, be due mainly to such per- 
turbations. 


11. The Indian Monsoon is not essentially different from any other 
geographically bound circulation system; they are all driven by differentia! 
heating and controlled by friction. 


If the earth’s surface were completely uniform there would be no reason 
for the formation of geographically bound climatological patterns, except 
zonal bands. The existence of climatological patterns, which differ from 
zonal bands, indicates the role that topography plays in the production of 
weather and climate. The importance of topography is not limited to fric- 
tional effects; of even greater importance are the topographical influences 
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on the location and intensity of the sources and sinks of heat, and the mecha- 
nical influences of topographical obstacles. 


In view of the complexity of the topography of the Indian region it is 
gratifying to observe that the Indian circulation system has been more 
intensely studied and more adequately described than any other geographically 
bound circulation system. Much of our knowledge of the Indian Monsoon is 
due to the researches of Dr. K. R. Ramanathan. On the occasion of the 
publishing of this volume, the writer has great pleasure in paying homage 
to Dr. K. R. Ramanathan for his numerous contributions to the science of 
the atmosphere. 
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ON A CLASSIFICATION WITH REGARD TO DAILY 
RADIATION INCOME FROM SUN AND SKY 


By Anpers ANGSTROM 
Received January 19, 1953 


Tue records of radiation income now obtained by means of various kinds 
of pyranometers—Callendar, Angstrém, Moll-Gorczynski, Robitzch a.o— 
give us the possibility to evaluate in general the total radiation income 
referred to periods of day, pentade, month and vear and in the case of the 
more accurate of the mentioned instruments also corresponding to hours 
or even shorter time intervals. The evaluation includes often a rather 
cumbersome and painstaking work, which can often not be peiformed im- 
mediately, but first after a certain lapse of time. On the other hand, the 
records seem to give many information concerning the character of the 
weather of the day, which are immediately apparent at a short inspection 
of the records, and it seems evident that a classification with due regard to 
these characteristics would add with a minimum of labour to the climato- 
logicai knowledge concerning the place in question. 


At the treatment of the records from Stockholm I have arrived at the 
following classification :—- 


I. Perfectly clear or almost clear days with depreciations in the 
record chiefly below 5 per cent. 
II. Partly clear days with frequent variations from clouds: 
(a) depreciations generally between 5 and 25 per cent., 
(b) depreciations generally larger than 25 per cent. 
III. Sun persistently shadowed by clouds in the first part of the day, 
later sun unshadowed or frequently shining. 
IV. Sun persistently shadowed by clouds in the later part of the day, 
earlier sun unshadowed or frequently shining. 
V. Practically no direct sun radiation: 
(a) income highly variable, 
(b) income low and steady. 


From the photographic records at Stockholm instances of these five 
cases are given in the Figs. 1-5. The following remarks may be added. 


1. It seems appropriate to include in this group days also where small 
depreciations have occurred, such as may easily be produced by accidental 
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Daily Radiation Income from Sun and Sky 


TABLE I 
Stockholm 


June 


1947 1949 1950 1951 1952 


December 


1947 1948 1949 1950 1951 1952 


1 3 


3 


13 (9 


increases in the turbidity of the air through smoke or dust. Such pheno- 
mena seem in general to be accounted for if we allow a limit of about 5 per 
cent. for the corresponding depreciations. 


Il. The group corresponds to cases which often occur chiefly in the 
summer. Convective clouds are then causing almost periodic deprecia- 
tions in the middle of the day, especially at the time when the temperature 
reaches its highest value. 


Most frequently the case occurs that cumulus clouds develop due to 
convection from the ground. But convection may also arise within higher 
atmospheric layers when a lower layer is heated through absorption relative 
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to an upper one. Even convective cirrus clouds (cirrocumulus) may in such 
cases arise. It scems appropriate to distinguish in the classification between 
cases with high convective clouds (il a), and low convective clouds (II 5). 


IlIl-IV. The groups correspond generally to cases where uniform cloud 
sheets are passing the station, as often occurs in connection with the passage 
of fronts. 


V. Corresponds generally to cases where the sky has been totally over- 
cast during the whole day. 


During such days the radiation income may have very different character. 
Figs. 5-8 show some characteristic cases, all corresponding to totally over- 
cast skies. The variations, which at times are ,almost periodical, are as 
regards their period naturally dependent as well on the size of the cloud 
elements as on their velocity in passing the sun. Density variations evidently 
always occur even in apparently homogeneous cloud sheets. The actino- 
metric records often suggest a regular wave formation giving rise to periodic 
variations in density or thickness of the cloud sheet. 


A classification like that suggested above may form a desirable comple- 
ment to the usual indications of cloudiness. If a more detailed classifica- 
tion should prove of value it may easily be worked out on the foundation 
here outlined. 


Meteorological and Hydrological Institute of Sweden, January 1953. 
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A RATIONAL APPROACH TO RAIN MAKING 


By HORACE R. BYERS 
The University of Chicago, Chicago, Illinois, U.S.A. 


Received January 26, 1953 


THE artificial nucleation of clouds to induce precipitation so obviously gives 
promise of great benefits to agriculture and related interests that it is already 
being tried commercially in a number of places before the scientific world 
has even approached agreement as to its efficacy. This means that care- 
fully designed experiments and tests still are needed in order to establish 
the facts. The simplest type of test, in which a cloudy area is seeded with 
artificial nucleating agents and an attempt is made to determine statistically 
whether or not the seeded clouds released more rain than they would have 
if left alone, has failed to produce answers that are universally accepted. 
Visual observations of what happens to individual cloud masses shortly 
after seeding reveal some spectacular modifications, particularly in the upper 
parts of the clouds, but tell little in a quantitative sense about rain produced 


in comparison with what might have resulted from undisturbed natural 
causes. 


A rational approach to the problem should consist of a physical as well 
as a Statistical study. To work intelligently with artificial methods one 
should know how Nature operates to produce precipitation and should have 
at hand statistical and physical information about the natural cloud popula- 
tion. 

PHYSICAL-CHEMICAL EFFECTS 

As a Starting point in examining the precipitation mechanism in natural 
clouds, one asks why rain or snow falls out of some clouds and not out of 
others. In the first place, it is recognized that precipitation is a different 
process from condensation and it occurs only under special conditions. In 
an ordinary cloud, as more and more water vapour changes to liquid through 
condensation, more and more droplets of the same size as existed before 
in the cloud are formed. Condensation does not occur on already existing 
droplets, so they do not grow beyond a certain initial cloud droplet size. 
For most ordinary clouds, this prevailing size is of the order of 5 to 20 
microns diameter and the size range is very narrow. According to observa- 
tions at the Mount Washington Observatory (U.S.A.F., 1948), a broad 
spectrum of size distribution is rare. Diem (1948) collected cloud droplets 
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on oil slides from an airplane and prepared composite curves of size distri- 
butions. For non-precipitating types of clouds he obtained the results 


Shown in Table I. H. J. aufm Kampe and H. Weickmann (1952) have 
reported similar results. 


TABLE | 
Diem’s Size-Distributions for Droplets in Non-Precipitating Clouds 


Mean diameter Range of diameters 
Gant: (microns) (microns) 


Fair weather cumulus 


8-5 
Cumulus congestus 14+5 3—40 


Stratocumulus 7-9 


Stratus 12-9 


Altostratus 10-6 


Droplets of these sizes are too small to have an appreciable velocity of fall, 
as shown in Fig. | a. Even if permitted to settle slowly in perfectly still 
air they would evaporate within a few centimetres below the cloud base, 
according to calculations of Findeisen (1939). Drops must have a dia- 
meter of the order of hundreds of microns to have an appreciable fall rate 
(Fig. | 6) and to be able to maintain themselves without evaporating as soon 
as they leave the cloud. 
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Fig. 1. Velocity of fall of droplets of various diameters. In (a) the droplets fall a 
accordance with Stokes’ law. The values in (6) and (c) were determined experimentally in 
still air near normal temperature and pressure by Gunn and Kinzer (1949). The sizes were 
computed from the mass to give an equivalent diameter assuming a spherical mass. At sizes 


in the upper ranges of curve (c) it is known that the drops are markedly nonspherical. 
The question to ask ourselves before we can be confident of the efficacy 


of seeding agents is: What is the manner of growth from suspended drop- 
lets (c1. 104) to falling drops (> 100) in natural clouds? Two processes 
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have gained acceptance. One is the Bergeron hypothesis (Bergeron, 1933) 
requiring the presence of ice crystals in the supercooled portion of the cloud. 
The other is the hypothesis of gravitational separation requiring a broad 
size distribution in the cloud so that there are large droplets, say with 60 
diameter, which fall with respect to the smaller ones, collide and coalesce 
with them to attain raindrop size. This process can account for the forma- 
tion of precipitation at temperatures above freezing, while the Bergeron 
process cannot. 


The best way to determine whether or not the development of rain takes 
place at temperatures above or below freezing and thus to see which process 
is likely to occur, is to use radar. A range-height indicating radar making 
vertical cross-sections through clouds all around the horizon will give the 
height of the rain echoes, and through height one can infer the temperature 
from radiosonde observations in the vicinity. If the radarscope is photo- 
graphed as the antenna makes its sweeps, one can go over the film and detect 
the first appearance of each echo and tabulate the temperatures of that echo. 
If the echo appears first at heights where the temperature is above freezing, 
one has some evidence that the precipitation-forming process occurred at 
above-freezing temperature. For, at ordinary ranges, the range-height 
indicating radar sets now available, such as the 3 cm. wavelength AN/TPS-10, 
first detect an echo when the drops have attained average diameters of 100 
to 200 microns, the precise value depending on the range and on the content 
of liquid water. The theoretical perception thresholds for this set are 
represented graphically in Fig. 2. 


In Fig. 3 is shown an example of a case studied by Battan (1952) in which 
a first echo of a growing cumulus cloud appeared with its top at a tempera- 
ture of 2°C. One might argue that there may have been ice crystals falling 
from an upper part of the cloud which could not be detected by the radar 
because of the lower dielectric constant of ice. This argument does not 
seem plausible for two reasons: First, experience has shown that such snow 
particles melt very near the freezing line and, as a matter of fact, their melt- 
ing is thought to be the cause of the well-known radar “ bright band ”’ 
which does not appear in these pictures. In the second place, as the series 
of photographs shows, the cloud top is growing upward at about 1,000 feet 
per minute, probably with an updraft of this magnitude, and even the heaviest 
snow particles have a velocity of fall of no more than 200 feet per 
minute. 


To verify the conclusions, one should have and could have in the case 
of an isolated cumulus cloud such as this one, visual observations of the 
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Fic. 2. Theoretical thresholds of perception by AN/TPS-10 radar set depending on 
droplet diameter and cloud water content. Curves from Battan (1952), computed on the assump- 


tion of uniform droplet size and complete cross-sectional interception of the radar beam 
by the cloud. 


cloud top or, better still, carcfully calibrated panoramic cameras at two points 
from. which the cloud heights could be measured by triangulation. 


If a large number of natural clouds examined in this way produce echoes 
before they have reached subfreezing temperatures, then rainmaking attempts 
based on the Bergeron hypothesis and requiring subfreezing cloud tops in 
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many cases will have come into the picture too late as far as the initiation 
of precipitation is concerned. ‘‘ Warm” cloud seeding, such as practised 
by Bowen (1952) in Australia might be more to the point. Bowen intro- 
duces a water spray which provides out-sized drops to start coalescence 
through gravitational separation. The case presented in Fig. 3, representing 
natural precipitation, could be a case where nature provides the out-sized 
drops. Langmuir (1948) has described how a “chain reaction” can be 
set off by introducing only a few large drops which grow so big they finally 
break, creating other drops which go through the same process. He has 
computed the collection efficiencies of drops in collision with smaller ones 
inthe cloud. The rates of growth in typical cloud situations through gravita- 
tional separation have been computed by him, by Houghton (1950, 1951) 
and by Bowen (1950). Once raindrops are initiated, by whatever cause, 
it is clear that this form of coalescence is the major factor in their further 
growth. 


If gravitational separation is a common cause of growth to 100-200 
micron sizes, then one must ask the question: How does the unusual cir- 
cumstance of a broad spectrum of size distribution come about? Some- 
thing must produce out-sized droplets in order for gravitational separation 
to begin. Woodcock (1952) collected some unusually large sea-salt particles 
from the atmosphere which could serve as nuclei for out-sized droplets. 
Some of them were so large that shortly after reaching a cloud at ordinary 
water saturation they would have diameters of 56 to 88 microns. They 
were few in number among a great many salt nuclei that under the same 
conditions would produce droplets of 14 to 40 microns diameter. Gravita- 
tional separation of these two size ranges could initiate and carry out precipita- 
tion growth by Langmuir’s “ chain reaction’. It would be unwise to accept 
the hypothesis of the large sea-salt nucleus without further experimental 
studies. There may be other substances of importance in this connection 
or other processes in warm clouds that contribute to droplet growth. Certain 
types of small-scale turbulence might increase the rate of collision and 
coalescence, or there may be some electrical effects. 


A study of atmospheric particulates and their effects requires both labo- 
ratory and field experiments and measurements. For fieldwork, techniques 
have been devised for examining atmospheric aerosal content both on the 
ground and from airplanes. For sea-salt particles there are the methods 
used by Woodcock, a special method of studying chloride particles by their 
chemical action on specially coated impactor slides or drums such as used 
by Crozier, Seely and Wheeler (1952) and the flame counter for sodium, 
described by Soudain (1951). The chemical industry has developed some 
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new filters which are ideally suited for collection and chemical analysis of 
aerosol particles. The variety of laboratory experiments suggested by field 
measurements and theoretical studies is almost without limit. 


If the ice-crystal mechanism is of predominant importance in initiating 
droplet growth, it is necessary to know at what temperatures ice crystals 
begin to form in natural clouds. The variability of these temperatures also 
should be determined. Then the effective temperature of crystal formation 
with various artificial nucleating agents should be found. Experimental 
determinations of these effective temperatures have been made, notably 
those by Schaefer (1951). If natural ice nucleation occurs at temperatures 
considerably lower than those for artificially-induced nucleation, then cloud 
seeding has a considerable range of usefulness in growing clouds. These 
effects also can be checked with radar. If a suitable technique for collecting 
ice particles from airplanes is used, direct measurements in the clouds can 
be made. Furthermore, the same experiments made on the temperature 
of spontaneous formation of ice crystals in laboratory air at saturation can 
be made in the free air by means of an airplane carrying a freezing box into 
which outside air is drawn. 


Another physical effect which has been noted in clouds as they begin 
to form rain is a marked change in the electrical field under such clouds, 
particularly if they are of a convective type. At the present state of know- 
ledge one suspects that this is an effect rather than a cause. Further 
measurements are necessary to see to what extent electrical effects play a 
role not only in the growth of droplets but in some of the ice crystal pheno- 
mena. High electric fields are known to exist in some snow conditions. 
The effect noticed by Workman and Reynolds (1950) when water freezes 
with certain contaminants in it should be pursued further. A combined 
set of field and laboratory experiments would be profitable in this 
connection. 


METEOROLOGICAL FACTORS 


Cloud seeding tests conducted up to now have shown varying degrees 
of success which appear dependent upon meteorological conditions—not 
just the temperature conditions or appearance and dimensions of the clouds, 
but larger-scale factors. If the clouds are in an area of widespread conver- 
gence and ascending motion the results are vastly different from those in a 
sinking, laterally diverging air mass [see Coons, et al. (1948, 1949, 1951) and 
Anderson (1950)]. Most clouds suitable for seeding tests probably are some- 
where midway between these two extremes. 
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The question to be asked here is: What determines the activity and 
liquid water content of clouds and thus their potentiality to produce rain ? 
For convective types, it is found that the amount of water vapour in the upper 
levels is important, or, in the language of synoptic meteorology, the presence 
or absence of moist tongues above about the 850-mb. level is critical. Also 
in developing convective activity one finds the occurrence of local regions 
of horizontal convergence, often too small in area to be detected by the 
existing network of upper-wind stations. An experiment designed to solve 
problems of rain inducement from convective clouds should include the 
establishment of upper-wind stations on a suitable grid for making the 
necessary computations of convergence. It might be necessary to establish 
through trial and error what the optimum station grid might be for detecting 
these effects. 


Another meteorological problem comes up in connection with attempts 
to seed clouds from chemical generators located on the ground. Most 
commercial operators use silver iodide smoke, depending upon the vertical 
stirring of the atmosphere to accomplish the transport of the silver iodide 
particles to the supercooled regions of the clouds. Whether or not the 
chemical reagent is carried to the proper levels by the air currents has never 
been carefully determined. Before going too far with claims of rain induce- 
ment, workers using this method should find out whether they are really 
treating the clouds in the manner claimed. Tests of smoke plumes, possibly 
using fluorescent particles for tracing, such as done by Braham, Seely and 
Crozier (1952) should be carried out. 


STATISTICAL DESIGN AND EVALUATION 


Most of the criticisms of artificial cloud nucleating work done so far 
have been in the nature of questions of statistical significance of the results. 


In many cases the tests are not so designed as to give statistically significant 
results. 


Precipitation conditions suitable for artificial modification can be 
classed in two types: (1) widespread general rains from solid cloud layers 
and (2) rain of the shower type, with more or less isolated convective clouds. 
The statistical analysis of the results would be different for the two types. 
A much-used method of statistical analysis has been to compare the season’s 
rainfall in a seeded area with that in a neighbouring unseeded region. A 
variant of this is to deal only with departures of the year’s rain from normal, 
making comparisons between a seeded area and a neighbouring “ control ”’ 
area. The significance of the results for a given number of years will depend 
on the homogeneity of the rainfall pattern in the region which includes the 
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seeded and the control areas. It turns out that in regions where rainfall 
is most needed and where rain making is most often tried, such as in the 
arid sections of the United States, the precipitation pattern is not homo- 
geneous, so that comparisons of this kind over a few seasons are not conclu- 
sive enough to satisfy the critics. An examination of rainfall records in the 
United States shows that the most homogeneous precipitation patterns both 
in time and space are found in the middle-latitude rain-forest area of western 
Washington and Oregon. There the comparison of rainfall totals may have 
some validity even though tested over only two or three years. 


For rain of the convective shower type, significant results can be obtained 
best by examining and treating the clouds as individuals. This requires 


much more knowledge of the physical processes in natural clouds than does 
the steady-rain problem. 


The United States Weather Bureau has announced a project under the 
direction of Mr. Ferguson Hall to test the effects of artificial cloud nuclea- 
tion in the adjoining areas of southwestern Washington and northwestern 
Oregon, the region of homogeneous rainfall. The seeding operations which, 
for greater certainty of treatment and in order to make in-cloud measurements, 
will be carried out by airplane, will be statistically designed to provide the 
greatest feasible sensitivity in obtaining quantitative results. Of the total 
number of seedable occasions, approximately half will be selected by random 
means for seeding operations. The remainder will be used as control cases, 
and results will be based on difference between the seeded and unseeded 
occasions. In an attempt to be sure of the portion seeded, the trajectories 
of the seeded cloud masses as well as the precipitation falling from them, 
will be followed by means of a tracer powder such as zinc sulphide, which 
fluoresces under ultra-violet light. 


The usual method of determining seeding effects is to start with the 
normal relationship, based on past records, between the rainfall of the 
target area and an adjacent control area. A regression is then formed for 
the two and evaluation of a seeding period takes the form of determining 
whether the normal historical relationship has been disturbed. Brier and 
Enger (1952) have shown that very different estimates in seeding effects 
result from different lengths of records used in establishing the regressions, 
as well as with different sets of adjacent control stations. To get away from 
this difficulty, the Weather Bureau will vary the procedure between seeded 
and non-seeded rainy days, using non-seeded rain days as controls. 


The type of evaluation just described is not suitable for the convective 
type of rainfall because of the extreme variability both in time and in space. 
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A more fruitful approach appears to be one in which the convective clouds 
are studied as individuals. The rain-gauge networks available are insuffi- 
ciently dense to give valid comparisons of rainfall produced from each cloud 
during its life-history, but radar can be used effectively. By means of radar 
it is possible to record whether or not a cumulus cloud produced an echo 
and therefore whether or not precipitation formation was initiated. In order 
to gain a sufficiently large sample of clouds for study, it is advisable to go 


hunting for them in airplanes, and to study them with radar mounted in 
an airplane. 


Let us say that, for seeding, a cumulus cloud must be at least 6,000 feet 
tall from base to top and its top must extend to levels where the temperature 
is between — 2 and — 10°C. These may not be the most logical figures 
to choose, but for the purposes of this discussion we will assume that they 
are. If one seeds some of these and rain results, the question may be asked 
as to whether or not they might have produced rain if left alone. To answer 
this, one may set up a programme where only half of the suitable clouds are 
seeded. The question of whether to seed or not to sced after a cloud meet- 
ing the dimensional and temperature requirements is chosen, may depend 
on the flip of a coin. Each of the clouds of the two categories, seeded and 
nonseeded, could be watched with the radar to see if an echo develops. One 
could take the results and place them conveniently in a two-by-two table 
of seeded and non-seeded versus echo and no echo. Since there might be 
a day-to-day variation in susceptibility to seeding, it might be well to work 
on at least one such pair of clouds each day. 


It is desirable to know how many clouds will have to be studied in order 
to get significant results. For this purpose it is necessary to know how 
many suitable non-seeded clouds develop echoes compared with the seeded. 
Suitable clouds are defined as those meeting the prescribed dimensional 
and temperature requirements. One determines the probability of echo 
formation from a suitable non-seeded cloud versus the same for a seeded 
cloud and statistical methods are available for determining the power and 
significance of any given number of trials. Since these probabilities are 
not known, considerable work lies ahead. Combined radar and camera 
installations at a ground centre can furnish information on natural clouds. 


A more difficult problem in connection with convective clouds is that 
of determining whether or not the rainfall can be increased, especially from 
clouds that are already producing rain. Perhaps the same type of approach 
as that used by the United States Weather Bureau on layer clouds will, from 
a prolonged series of tests, lead to valid conclusions, At the present stage 
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of investigation, the problem of initiation of the precipitation mechanism 
appears more reasonable of solution. By making careful measurements 
of time required to develop an echo after reaching suitable status in both 
seeded and non-seeded clouds, one can determine whether or not the seeding 


hastened the process. 


In conclusion, it should be pointed out that there is still a vast amount 
of fundamental work and rigid testing to be done before a scientific evalua- 
tion of rain making can be accomplished. 
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Fic. 3. Series of photographs of radarscope for vertical sweeps through a developing 
cloud at a range of 47 miles, from 12h 50m 48s to 12h 53m 53s on 14 Aug. 1947 
The approximately horizontal white lines mark altitude at intervals of 5,000 feet and the 
faint, dark vertical lines at the sides of each picture are the range markers of 40 and 50 miles 


The first echo at 12: 50: 48 has its top at 13,500 feet where the temperature is 2° C. 
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1. 


INTRODUCTION 


IN infra-red spectroscopy one is concerned with the detection and measure- 
ment of extremely minute quantities of radiant energies, which are available 
at the exit slit of the spectroscope after undergoing dispersion by the optical 
system. This minute energy is detected by its heating effect by means of 
sensitive thermocouples and measured by a high sensitivity galvanometer. 
This detecting arrangement which is basic to almost all types of infra-red 
spectroscopes, is subject to certain limitations in respect of sensitivity 
as well as of the speed with which the measurement can be carried out. While 
the Brownian motion sets the ultimate limit to the sensitivity which can be 
attained in the system, the response time of the thermocouple and the period 
of the galvanometer are the limiting factors for the speed of measurements, 


The development of the Pfund resonance radiometer! enabled a great 
reduction to be made in the effects of drift due to changes in ambient tempera- 
ture and air currents and also due to the Brownian motion. In Pfund’s 
resonance radiometer this is achieved by periodically interrupting the beam 
of radiation to be measured by means of an electro-magnetically maintained 
pendulum operated shutter, at a frequency at which the detecting galvano- 
meter resonaies; thus the system becomes more sensitive to impulses of the 
period to which the system resonates than to the random fluctuations of 
the Brownian motion and drifts due to changes in ambient temperature and 
air currents. Firestone? has further developed this idea in his periodic 
amplifier, and his recorder, which has a period of 3-6 seconds as compared 
to 140 seconds used by Pfund in his resonance radiomeier, has been very 
successfully used by various workers (viz., Randall and Weber and others), 
and is widely used even to-day. Judging from the beautiful records obtained 
with this recorder, while there can be no question about the sensitivity 
attained by these methods, the limitations in the speed of recording imposed 
by the galvanometer system are serious disadvantages, especially as it excludes 
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the use of this recorder for studying rapid phenomena, viz., in solar spectro- 
scopy and in the study of kinetics of chemical reactions, etc. 


Much progress has been made recently in respect of the speed of record- 
ing by the development of self-baiancing potentiometer type or ‘servo’ 
type of recorders. Here again the greater advantage is in obtaining a pen 
and ink record than in the speed of the recorder which is again limited 
by the electromechanical balancing system of the recorder. 


It is proposed to describe here a new technique of high speed recording 
for use with infra-red spectroscopes, which has been developed by the writer 
in this laboratory. In this new system the speed of recording is limited only 
by the response time of the radiation detector and not by the galvanometer 
system as in all other systems. As modern infra-red radiation detectors, 
yiz., the fast thermocouple, the thermister bolometer, the Golay pneumatic 
detector, are being made with response time which are only 1/1000th of 
those used by Pfund and Firestone, and the advances made in photo- 
conductive cells which can respond to extremely high frequencies and can 
be used in the infra-red region up to about 9-0,, it will be obvious that 
the new technique achieves an improvement in the speed of recording by 
a factor of 1000 or even higher if used with photo-conductive cells. 


2. THE NEW RECORDING SYSTEM 


In modern infra-red spectroscopes the infra-red radiation falling on the 
entrance slit is usually modulated mechanically at a frequency which can be 
accommodated by the response time of the detector used, which may be a 
thermocouple or a bolometer. This modulation is carried out in order that 
the detector may yield an alternating e.m.f. as output which greatly simpli- 
fies the problem of electronic amplification and also enables discrimination 
to be made against the background radiation from the walls of the spectro- 
scope which becomes almost comparable to the signal at longer wave-lengths, 


The frequency at which the beam is modulated depends on the response 
time of the detector, and in most present-day equipment this varies from 
Sc/s to 25 c/s. Although with some of the new detectors the use of modu- 
lating frequencies of several hundred cycles per second is possible, the above 
low frequencies are chosen as the use of higher frequencies will be of little 
value so long as deflection systems are slow. As in the new recording system 
the mechanical deflection system has been replaced by a purely electronic 
one, i.e., a cathode-ray tube, the highest modulating frequency which can 
be accommodated by the speed of the detector can be used. In this labo- 
ratory, the new technique is at present being used in conjunction with a 
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Beckman Model IR-2 infra-red spectrophotometer which uses a vacuum 
thermocouple and has a modulating frequency of 10c/s. Even with this 


low frequency and slow thermocouple, an improvement in speed by a factor 
of 50 to 100 has been achieved. 


Fig. | shows a schematic diagram of the new recorder. In this recorder, 
the output of the detector of the light modulated infra-red spectroscope, 
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Fic. 1. Schematic Diagram of a high-speed recorder using cathode-ray tube 
for use with infra-red Spectrescope 
which consists of an approximately square wave-shaped alternating e.mf. 
of the modulating frequency (10 c/s in this instrument), is first fed to a filter 
network with a narrow band pass characteristic for the modulating frequency. 
The signal is now passed on to a high gain tunned amplifier, which is a phase- 
shift oscillator working below its oscillation point. The details of this 
narrow band amplifier (Fig. 6) will be given later on in this paper. As this 
amplifier can be made to have a band width as small as 1 c/s, it gives a very 
high signal to noise ratio. This is a very great advantage in electronic ampli- 
fying systems where the limit of ultimate usable amplification is set by the 
thermal noise in the input circuit. This advantage may be considered 
equivalent to that obtained in the Pfund Resonance Radiometer in respect 
of Brownian motion. The signal is now fed to a straight resistance-capacity 
coupled amplifier which works into a phase inverter stage fof providing 
push-pull deflection to the horizontal plates of a cathode-ray tube. 


The trace obtained on the fluorescent screen of the cathode-ray tube 
consists of a horizontal trace whose amplitude is proportional to the amount 
of radiant energy falling on the detector of the spectroscope. A camera 
with a continuously moving 35 mm. film records this trace as the wave-length 
drive of the spectroscope is rotated by a synchronous motor which is linked 
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to the motor drive of the camera. The record obtained on the film is thus 
a spectrogram showing the amplitude of the trace which is proportional to 
the energy falling on the detector, as a function of the wave-length. It is how- 
ever not essential to use a camera for this purpose as a roll of film or photo- 
graphic paper moved in contact with the face of the fluorescent screen would 
give the same result. 


A few typical records of the infra-red spectrum and the emission of a 
Nernst glower are given in Plates VIII and IX, Figs. 2, 3, 4 and 5. These 
records were each made in about 12 to 15 seconds. It may be mentioned 
that it will take from 10 to 15 minutes to record the above spectrum with 
‘servo’ type of strip chart recorders and a few hours with the Pfund reso- 
nance radiometer and periodic galvanometer. 


ORWE 


3. THe TUNED AMPLIFIER 


The narrow band tuned amplifier which gives a great advantage by 
improving the signal to noise ratio by giving a discrimination against the 
unmodulated background radiation as well as the random thermal noise 
originating in the detector and the input circuit of the amplifier, is a very 
important part of this equipment.. This amplifier actually consists of a' 
phase shift oscillator working with its feed back adjusted to a level just below 
the point where it goes into oscillation. 


ook 


The basic schematic circuit diagram of this type of amplifier is given in 
Fig. 6. The frequency to which the amplifier is tuned is given by the relation: 


1 
= 
where f, is the resonant frequency and R and C are the resistance and capa- 
city used in the phase-shifting network. The phase shift produced by the 
RC network is such that the feed back is positive only at f, the resonant 
frequency and negative at all other frequencies higher and lower than fy. 
The amount of feed back can be controlled by changing the screen grid 
_ Voltage by adjusting the value of Rg and R;. When the feed back is adjusted 
to a point just below that at which oscillations occur the response of the 
amplifier is of the type shown in Fig. 7. Even a single stage of this type 
of amplification in the low level stages of the amplifier gives a tremendous 
increase in the signal to noise ratio. Additional stages will of course make 
the response sharper and can be advantageously used. 


4. CONCLUSIONS 


The technique of recording the infra-red spectrum at high speed extends 
the application of infra-red spectroscopy to rapid phenomena which were 
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Fic. 6. Phase Shift Tuned Amplifier 
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Fic. 7. Response of Tuned Amplifier 


hitherto beyond the range of this useful technique, which enables the study 
of molecular structure and reactions in a manner which is simple, elegant 
and of very great accuracy. As an example may be quoted the application 
of this technique to solar spectroscopy. Solar spectroscopic measurements 
were possible only when the sky was clear enough to give bright sunshine 
for the whole period required for making one complete spectrogram which 
could be of the order of hours with the older technique. While with the 
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technique described here, the number of days on which such records can 
be obtained is very greatly increased as it is possible to make a record even — 
through openings in clouds, if the sun shines even for half a minute. Unless 
the sky is completely overcast, it is always possible to obtain some records. 


The use of this technique in the study of the kinetics of chemical reac- 


tions is obvious. Many other problems will no doubt be amenable to study 
by this technique. 


Another advantage of this technique is that it greatly increases the 
output of work of an instrument and this is of importance since infra-red 


spectroscopic equipment is very expensive, and a single apparatus becomes 
available for more work than is possible otherwise. 


In conclusion, the writer wishes to thank Dr. L. A. Ramdas, the Investi- 
gator-in-chief, Radiation Scheme, and the Director of Agricultural Meteoro- 
logy, for his kind interest and encouragement during the course of this work. 
The writer’s thanks are also due to the Council of Scientific and Industrial 
Research who are financing the Scheme in which the writer is at present 
working. The writer is also grateful to Shri S. P. Venkiteshwaran, Met- 
in-charge, Instruments and Supplies, Poona, and the Staff of the Meteoro- 
logical Office Workshop for facilities and help. 
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1. INTRODUCTION 


Tue infra-red spectrum of magnesium oxide in the form of thin cleavage 
. plates as well as of thin films fumed on to nitrocellulose films and NaCl 
plates, has been the subject of study by several workers. While there appears 
to be some disagreement between the earlier workers in the results obtained, 
recent work has brought out a more consistent picture of the spectrum of 


magnesium oxide. The results of previous workers may be summarised as 
follows :— 


_ Tolksdorf (1928)! studied the transmission of thin layers of powdered 
MgO held between thin plates of NaCl and found absorptions at 3-85, 7-65 
and 14:2 microns. The residual ray technique of Strong (1931)? has revealed 
the existence of a strong absorption in the region of 23 microns in fumed 
films of MgO and a strong reflection of almost 80% in the same region in the 
case of a freshly cleaved plate of MgO. Barnes, Bratain and Seitz (1935) 
who made a detailed study both of thin cleavage plates as well as of films 
of MgO evaporated on to nitrocellulose films in vacuum found a considerable 
amount of fine structure consisting of a number of weak lines between 6 
and 12 microns. Fock (1934)* studying films fumed on to nitrocellulose 
found a strong absorption maximum at 17-3 microns and two weaker 
absorptions at 14-9 and 25 microns. Of the recent studies, Elias Burstein, 
J. J. Oberly and Plyler (1948)° have studied the absorption spectra of thin 
cleavage plates between 5 and 14 microns and the reflection spectra between 
15 and 37 microns. They report absorption bands at 8-3, 10-18 and 11-82 
microns and a reflection which starts at about 13 microns and rises rapidly 
to a peak of almost 100% at about 21 microns except for a small kink between 
15 and 16 microns and decreases rapidly beyond 24 microns and is about 
40% at 37 microns. Wilmott (1950)® has studied both plates and films. 


The work of Strong clearly establishes the selective reflection of MgO 
about 23 microns, but his reststrahlen technique is not capable of giving any 
details regarding the structure of these bands and only indicates the rough 
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position of the band. The measurements of Barnes, Bratain and Seitz 
appear to suffer from errors which may be attributed to instrumental factors, 
viz, an unregulated infra-red source or their method of making point by 
point measurements. Although their general results have been confirmed by 
later work, the considerable amount of secondary structure reported by 


them has not found confirmation in recent work carried out with more refined 
techniques. 


In view of the lack of agreement on some points between the various 
previous workers, and the fact that a new model of infra-red spectrophoto- 
meter which the writer felt was a more advanced instrument than those used 
by other workers, became available to him, it was felt that it may be worth- 
while to repeat these measurements of MgO films. Another important 
reason for undertaking this investigation was that thin layers of MgO are 
in use in this laboratory for the purpose of filtering out the visible portion 
of the spectrum in an instrument for the measurement of infra-red radiation 
from the atmosphere developed by Kale® and it was important to know the 
exact transmission characteristics of this material in the form of thin films 
with various thicknesses in the infra-red region between | and 20 microns. 


Before presenting the results obtained it may be useful briefly to describe 
the experimental arrangement used in these measurements as it will help 
in making a comparative study of the measurements made with other types 
of instruments by previous workers. 


2. EXPERIMENTAL 


The Beckman model IR-2 Infra-red spectrophotometer used in this 
investigation uses a potassium bromide prism and optics and covers the 
range from | to 25 microns. The lay-out of the instrument is indicated in 
Fig. 1. The optical path of the infra-red radiation through the spectro- 
photometer is as follows and is indicated by an arrow line in Fig. 1: 
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. Schematic optical path of Beckman Model IR2 infra-red spectrophotometer 
(Courtesy of National Technical Laboratories) 


Light from the Nernst glower A is modulated at 10 cycles per second 
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by the condensing mirror D, forming an image of the source in the focal 
plane of the concave KBr lens E. The collimated beam passes through the 
filter slide and shutter F and the liquid cell compartment G. The lens H 
converges the beam through the gas cell compartment I and the lens K, form- 
ing an image of the source on the entrance slit L. Light from the slit is colli- 
mated by the spherical mirror M (5-5° off axis) and dispersed by the KBr 
prism N (6cm. base and 4cm. height). It is reflected by the rotatable 
Littrow mirror O, dispersed a second time by the prism and focussed by the 
collimator M on the plane of the exit slit, which lies just below the entrance 
slit. The beam is directed by the plane mirror J against the condensing 
mirror P which focuses an image of the exit slit on the thermocouple Q. The 
phototube R monitors the source to maintain constant output; its received 
light is adjustable by a shutter so that the glower current may be regulated 
to the proper value. 


The radiation source is a Nernst glower operating on 60 Watts of power 
at about 0-8 ampere. The glower is self-starting and its output is auto- 
matically held constant within 0-1% by the special photoelectronic monitor. 
The watercooled housing maintains a constant temperature. 


Two interchangeable shutters, one of metal for use up to 9 microns and 
the other of glass for use beyond 9 microns, are provided for modulating 
the beam. The glass shutter eliminates the effects of contamination of the 
longer wave-lengths by the scattered short-wave radiations which are con- 
centrated in the energy peak of the glower near about 2 microns. Since 
this scattered radiation is not modulated by the glass shutter it is ignored 
by the tuned amplifier which only responds to radiation modulated at 10 c/s. 


3. PREPARATION OF THE MgO FILMs 


The MgO films were prepared by fuming on to a KCI plate from a burn- 
ing magnesium ribbon. The measurements were first started with an 
extremely thin layer which was just visible and then the thickness increased by 
adding more layers of MgO on the same plate by fuming. As the KCI plate 


is almost completely transparent up to about 21 microns these measurements 
cover the region | to 21 microns. 


4. DISCUSSION OF THE RESULTS 


Measurements have been made on 14 different thicknesses starting 
from a layer which was just visible to a thickness of about 2-4 mm. when 
the film begins to break away from the KCI plate. Fig. 2 shows four curves 
giving the infra-red spectrum of MgO films for selected thicknesses, for the 
region | to 21 microns. The other curves have been omitted for the purposes 
of clarity as they lie in between the curves given in this figure. 
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The most obvious conclusion that can be drawn from a study of these 
curves is that infra-red absorption spectrum of MgO is not a continuous one 
but consists of a series of well defined bands. The strongest of these absorp- 
tion bands lies between 18 and 19 micron. There is another strong absorp- 
tion band between 14 and 16 and it is clear that this is not a single band but 
consists two or more bands very close together. As we go towards the shorter 
wave-lengths a number of bands can be seen developing as the thickness of 
the film is increased, until we reach a region near about 7 microns where the 
transmission begins to decrease rapidly. This is obviously due to the reflec- 
tion and scattering of the shorter wave-lengths by the MgO film. More 
broadly speaking a film of MgO has a broad transmission from about 1 to 2 
microns to beyond 21 (except for a band at about 18 microns) with its peak 
at about 10 microns and its short and long wave-length cut off points are 
functions of the thickness of the film. This transmission characteristic of 
MgO films is of great use in measurements of atmospheric radiation in the 
infra-red region, during daytime where we are dealing with emission sources 
with effective temperatures of about 300° A which have their energy peaks in 
the 10 micron region. The use of MgO films by Kale® for his instrument 
for the measurement of atmospheric infra-red radiation appears to be fully 
justified, if a film of correct thickness is chosen so that the short-wave cut off 


257 


358 A. U. Momin 


is at about 2 to 3 microns and the 18 microns band is not very strong. It 
may be pointed out that MgO films reflect and scatter away almost 
100% of the energy in the visible region while remaining almost com- 
pletely transparent in the 10u region and therefore if a layer of this material 
is coated on a blackened thermopile or other receiver, it becomes sensitive 
to the atmospheric infra-red radiations while remaining completely un- 
affected even by bright sunlight. 


5. APPLICATIONS TO THEORY OF VIBRATION SPECTRA OF CRYSTALS 


The infra-red spectrum of magnesium oxide is of great interest to the 
crystal physicist. According to the theory of vibration spectra of crystals 
put forward by Sir C. V. Raman,’ a magnesium oxide crystal has nine 
fundamental modes of vibrations whose frequencies are listed below:— 


| 
Frequency in cm.-? 704 | 680 | 652| 584| 498| 184 


Wavelength in microns -| 14-2 | 14-7 | 15-35] 17-1 | 19-0 | 21-1 | 23-4 | 38-8 | 54-4 


| 
| 
Assignment | v1 | V4 Vs % Vr Vs V9 


Of these all frequencies except the last two, siz., 38-8 and 54-4 microns 
are expected to be infra-red active. 


The results of the present investigation permit an experimental verifica- 
tion of these discrete vibration frequencies predicted by the theory. As 
the region covered in this investigation is from | to 21 microns only the 
23-4 band cannot be verified by the present measurements, the last two fre- 
quencies being inactive in the infra-red. 


Fig. 3 gives a chart showing the predicted active fundamentals for 
magnesium oxide on the basis of Raman theory and the clear absorption 
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bands found in these measurements. As there are a number of well- 
developed bands in the overtone region these are compared with the second 
harmonics of the fundamentals. A fairly good agreement between theory 
and experimental result is obtained. There are indications that some of 
the bands have additional structure which is either beyond the resolving 
power of the instrument used or is incapable of being resolved on account 
of the overlapping or broadening of the bands due to coupling between the 
modes of vibrations. 


The author intends to further investigate this region with a NaCl prism 
which has a higher dispersion in this region than the KBr used in the present 
measurements and with a new recording technique which it is hoped will 
resolve the additional structure of which there are strong indications, if it is 
resolvable at all. 
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Ir has been considered for some time that measurements of electrical poten- 
tial gradient may provide an exact method of identification of air masses 
in the upper air, but it could not be measured so far due to the want of a 
simple and reliable device for the same. A number of methods are available 
for the measurement of the electrical potential gradient near the ground, 
but, in most of them the equipment used is not easily portable and cannot 
therefore be adapted to measure the distribution of this element in the upper 
air with the help of balloons. A simple method of using a tetrode valve 
as an electrometer was devised by R. P. Lejay! and used by Idrac? for 
measuring the electrical potential gradient in the atmosphere. This was 
not found satisfactory as the losses were considerable. In this system the 
controlling electrode of the valve, namely the grid, was in the electron stream 
from the filament to the plate and can therefore account for the losses. 
Recently R. E. Belin® has described another system employing discharge 
points in association with the audio frequency modulated type of radiosonde 
used in the U.S.A., but objections can be raised on these methods also, as 
they cannot function satisfactorily when the potential gradient varies 
suddenly ; the temperature of the air will also have an effect on the discharge 
points. 


INTRODUCTION 


L. Koenigsfeld and Ph. Piraux* have described a simple portable electro- 
meter using the HL 23 valve to measure the atmospheric potential, and have 
shown how it can be adapted for use with the type of radiosonde employed 
in Belgium. The transmitter in the radiosonde described in this paper was 
working on a wave-length of about 10 metres and the frequency of modula- 
tion of the carrier wave gave successively the measurements of temperature, 
pressure and humidity. The frequency modulation due to the changes in 
the meteorological elements is caused by varying the inductance due to the 
changes in the air gap. In adapting this instrument to measure the electrical 
potential gradient in the atmosphere with their valve electrometer, they kept 
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the air gap constant, and the coefficient of self-induction of the iron bobbin 
was controlled by the current from the valve electrometer. The current 
from the valve electrometer depended upon the potential gradient. 


A radiosonde equipment of the type used by the U.S. Army during 
World War II was available at the Meteorological Office at Poona. The 
meteorological elements are measured with this type of instrument from the 
audio frequency modulation of a carrier wave on about 4 metres (72-2 mcs.). 
The signal is modulated by varying the resistance in the meteorological 
control circuit by means of resistors. An attempt was made to see whether 
this type of radiosonde can be used with the valve electrometer devised by 
Dr. Koenigsfeld to measure the electrical potential gradient in the atmosphere. 


The method adopted and the observations made at Poona on a few days 
are described in this paper. 


2. THE VALVE ELECTROMETER OF DR. KOENIGSFELD AND PIRAUX 


The electrical potential gradient in the atmosphere can vary near the 
ground in the clear weather between 10 and 200 volts per metre and the electro- 
meter employed for this purpose should consume only an extremely feeble 
current. Generally, lead nitrate fuses or radioactive substances like polo- 
nium or ionium are used as “collectors”? for measuring the electrical 
potential, and the air-collector resistance is of the order of about 5 x 10!° ohms. 
Therefore the smallest leakage current in the measuring apparatus would 
cause voltage changes that would completely falsify the measurements. 
Koenigsfeld and Piraux have shown a method of measuring these potentials 
with a Mazda HL 23 valve, by keeping the leakage losses at a minimum value. 


Firstly, they removed the base of the valve whose insulation was found 
to be unsatisfactory. 


Secondly, there can be losses due to the grid current in the valve depend- 


ing on the mode of use. The chief causes when used as Type A amplifier 
are: 


(1) Electronic grid current due to the initial speed of electrons which 
teach the grid insufficiently polarised. 


(2) lonic grid current due to the attraction of a negative grid on positive 
ions produced by the collision of electrons emitted by the cathode against 
the molecules of residual gases. 


(3) Electronic grid current. 


(4) Photo-electric current generated by external light or by X-rays 
emitted by the anode at the time of impact of electrons, 
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It is clear from the above that the valve cannot be expected to function 
as an electrometer in the normal way. Koenigsfeld and Piraux have modi- 
fied the use and have employed the plate instead of the grid as the governing 
electrode of the valve. They have also shown that by this innovation, the 
leakages referred to above can be reduced to a negligible value. The only 
extra precautions to be taken in designing the apparatus are in enclosing 
the valve in a light tight box to reduce photo-electric effects and by employ- 
ing polysterene wires to obtain suitable insulation. They have also recom- 
mended the use of the acorn triode 958 A in which the leakages by insulation 
defects were less than in the case of HL 23. 


3. THe METHOD ADAPTED FOR USING THE VALVE ELECTROMETER 
WITH THE AMERICAN RADIOSONDE TO MEASURE THE 
ELECTRICAL POTENTIAL GRADIENT 


The American radiosonde consists of two oscillators operating on a 
frequency of about 1 mc. and 70 mcs. respectively (Fig. 1). The 1 me. 
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oscillator is the modulating oscillator (M) and is of the squegging type. 
The other is the carrier oscillator (C), also of the squegging type and is 
“on” when M is “ off’ and vice versa, thus in effect producing a frequency 
modulation of the carrier oscillator. The modulating frequency of the 
oscillator M is controlled by the RC combinution in the grid circuit, Hence 
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by varying the value of the resistance R in the grid circuit, the modulating 
frequency is altered. In the American radiosonde the signal is modulated 
by varying the resistance in the meteorological control circuit by means of 
resistors sensitive to temperature and relative humidity. It also contains a 
baro-switch with a contacting arm moving on a commutator due to the chang- 
ing pressures and it indicates the pressure values during the sounding and 
switches into the control circuit in a definite order, the temperature, humidity 
and standard reference resistors. 


In adapting this radiosonde for use with the valve electrometer, the 
resistances for measuring the temperature and humidity were removed and 
in their place a negative bias developed in the valve electrometer circuit was 
applied in the grid circuit of the modulating oscillator M. The modulating 
frequency depended on the amount of bias applied. 


A positive voltage of the order of 6 volts is applied through a resistance 
(R) of about 350 ohms in the grid circuit of the HL 23 valve (Fig. 1). The 
voltage and resistance in the grid circuit being fixed, the current through 
R will be controlled by the negative potential applied on the traditional plate. 
If this bias voltage is applied in the grid circuit of the modulating oscillator 
M of the American radiosonde, the modulating frequency of the carrier 
oscillator will be a measure of the negative potential on the plate circuit 
of the HL 23 valve. 


To measure the electrical potential gradient in the atmosphere, the two 
collectors (either lead nitrate fuses or polonium collectors) were fixed on 
polysterene insulators and hung at a known distance of one or two metres. 
The lower collector was connected to the plate of the HL 23 valve and the 
upper collector to the negative side of the filament. As generally in the 
atmosphere, the electrical potential increases with height, the effective 
potential on the plate of the HL 23 valve will be negative, and the modu- 
lating frequency of the carrier oscillator will be a measure of the potential 
gradient in the atmosphere. 


In an actual sounding, the HL 23 valve and its leads were well insulated 
with polysterene sheets and alkathene sleevings and placed in a cardboard 
box with holes at proper places for the leads and fixed inside the American 
radiosonde in the space for the temperature and humidity elements, and 
covered with the flaps to reduce external light to a minimum. Before 
releasing the instrument it was calibrated by applying different known 
voltages from a H.T. battery to the HL 23 valve at the collector terminals 
and the audio frequency modulation measured on the radiosonde recorder. 
A calibration curve is drawn with the data thus obtained. The baro-switch 
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on the radiosonde is also adjusted suitably after reading the atmospheric 
pressure on the ground. 


As the modulated frequency which is a measure of the electrical potential 
gradient, depends on the voltage developed across the resistance R (350 ohms) 
in the grid circuit of the HL 23 valve, it is essential to see that the value of 
this resistance does not change with the large fall in temperature experienced 
during the ascent. This resistance was therefore specially wound with 
constantan wire. 


It is also necessary that the batteries used in the sounding are of suitable 
capacity. For this purpose light acid batteries, after charging suitably, 
were used with the radiosonde transmitter. 


Both lead nitrate fuses and polonium were used during the few ascents 
so far made to measure the potential gradient. Dr. Koenigsfeld was kind 
enough to spare a few polonium collectors for our use. These collectors 
were of pure silver plates, about one square centimetre in area and about 
1/10 mm. thick, dipped in a solution of Radium D, and the emission from 
these collectors were estimated to be more than 40,000 per minute. They 
were fixed in polysterene sheets and separated by a distance of two metres, 
one vertically below the other. The leads from these were also covered with 


alkathene or polysterene sleevings to keep the leakage losses at a negligible 
value. 


When polonium collectors were not available, burning lead nitrate fuses 
were employed. These fuses were made by dipping narrow strips (about }" 
wide, }” thick cardboard pieces, about 2’ in length) in a solution of 5% 
lead nitrate and drying them. Two of these strips were tied together from 
opposite sides of a steel wire, supported on polysterene on a wooden frame. 
Two such units, at a distance of two metres, formed the collectors. They 
were ignited just before release, and laboratory trials showed that these fuses 
could burn for nearly half an hour. 


But it may be mentioned in this connection that polonium collectors 
should be preferred to lead nitrate fuses, as they are more handy and can be 
expected to function in any weather. Moreover, polonium collectors are 
more active and will attain the surrounding potential much faster than lead 
nitrate fuses—2 to 3 seconds as against 8 to 10 seconds. The radiosonde 
and collectors were suspended from the balloon by a chord of 50 metres 
with polysterene insulators in the middle in order to isolate the radiosonde 
from the balloon which can acquire electrical charges by friction in 
the air, 
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During the ascent, the baro-switch on the radiosonde keeps the valve 
electrometer connected to the transmitter except occasionally when the con- 
tact arm connected to the aneroid moves to the commutators corresponding 
to previously calibrated pressure values, when the modulated signal at a 
fixed frequency of about 180 c/s. is received. Thus a record of the modulated 
frequency due to electrical potential gradient is obtained almost continuously 
with occasional breaks of short duration when the pressure signals are 
obtained. From these one can easily evaluate the variation of potential 
gradient with height. 
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To record positive potential gradient, the plate of the HL 23 valve is 
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gradients are to be m:asured, the connections are reversed. But to measure 
positive and negative potentials during the same sounding, a battery, say of 
100 volts, is introduced between the upper collector and the plate of the 
HL 23 valve (positive of the battery to the collector). In that case the 
potential difference across the plate and filament of the HL 23 valve will be 
the difference between the battery potential and the atmospheric potential, 


Fig. 2 shows the records obtained at Poona during soundings on the 
16th and 22nd December 1952. The ascent on the 16th was with lead 
nitrate fuses two metres apart. At the time of ascent the sky was covered 
4/8 with altostratus and 2/8 with cirrus. On the 22nd December, polonium 
collectors were used at a distance of two metres and at the time of ascent 
there was a trace of fair weather cumulus. 


—— LEAD NITRATE FUSE (WO. 
COLLECTOR (0.2) 
POONA 
650/- AT 1646 HRS 
(2) AT 4649 HRS 1.5.7, 
” 
z 
i 
650 | a 
| 
20 40 60 680 100 V 
POTENTIAL GRAQIENT in 2) 
Fic. 3 


Fig. 3 shows the distribution of potential gradient obtained from two 
soundings made consecutively with polonium and lead nitrate fuses on the 
25th November 1952, when the skies were clear, The variation in the 


the 
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potential gradients are similar and the differences may be due to rapid fluctua- 
tions which occur towards the evening. 


The present paper represents the results obtained from preliminary trials 
and further experiments are in progress to improve the technique and to 
obtain regular data. The authors are grateful to Dr. L. Koenigsfeld for 
his kind interest and help in developing the technique. 


REFERENCES 
1. Lejay, R. P. .. Comptes Rendus, 1924, 178, 2871. 
2. Idrac .. Mem. de l’office National Met. de France, No. 18. 
3. Belin, R. E. .. “A radiosonde method for atmospheric potential gradient 
measurements,” Proceedings of the Physical Society, 
London, 1948, 60, 381-87. 
4. L. Koenigsfeld and “Un nouvel electrometre portatif pour la mesure des 


Ph. Piraux charges electrostatiques par system clectronique,” Inst. 
Roy. Met. de Belg. Mem., 45. 


4 

d 

n 

t 

two 
the 


SPECTROPHOTOMETRIC STUDY OF SUNSPOTS 


By R. ANANTHAKRISHNAN, F.A.Sc. 


(From the Solar Physics Observatory, Kodaikanal) 
Received December 23, 1952 


INTRODUCTION 


THE well-known paper by Prof. H. N. Russell (1921) entitled “* Cooling by Ex- 
pansion in Sunspots ”’ is the first attempt towards a quantitative explanation 
of the observed fact that sunspots are some 1500° to 2000° K. cooler than 
the surrounding photosphere. Russell considers the sunspot as a region 
whose equilibrium is governed by the equations of adiabatic change of state 
while the surrounding photosphere is assumed to be in radiative equilibrium. 
The observed difference of temperature between the spot and the photo- 
sphere is then made use of to calculate the pressure, temperature and density 
at the depth from which the spot is supposed to originate. Retaining the 
basic idea of Russell but modifying his simple theory in some important 
respects Petrie (1930) and Milne (1930) made calculations of the tempe- 
rature and pressure at the base of sunspots and the depth below the photo- 
sphere of the base of the umbral column. From the horizontal pressure 
gradients between the umbra and the adjacent photosphere, Milne also made 
estimates of the velocity of horizontal flow of gases in sunspots first observed 
by Evershed (1909) at Kodaikanal. 


In two important papers entitled ** Konvektion in der Sternatmosphare,” 
Unsdld (1930, 1931) showed that in a stellar atmosphere composed mainly of 
hydrogen, there should be a zone just below the surface where convection 
currents can be sustained because of the super-adiabatic lapse rate prevailing 
in this zone consequent on the rapid increase with depth in the ionisation of 
hydrogen. Under plausible assumptions Unséld found that in the solar 
atmosphere convective instability would set in at a depth of about 100 km. 
below the surface and continue downwards for the next 100 km. He 
identified this zone as the source of sunspots and granulations and also calcu- 
lated the minimum possible cooling in sunspots. 


In 1930 Pettit and Nicholson of the Mt. Wilson Observatory published 
the results of measurements of the ratio of the energy in the umbre of spots 
to that in the adjacent photosphere over the wavelength interval 0-4, to 
2:2. From this they derived the spectral energy curve of sunspots and 
deduced a mean black body temperature of 4750°K. for the sunspot as 
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against 5955° K. for the centre of the solar disc. They also measured the 
total energy in the umbre of spots in terms of that for a corresponding area 
of the adjacent photosphere and found a value of 0-471 for this ratio. 


The experimental results of Pettit and Nicholson formed the starting 
point of an important theoretical investigation by Minnaert and Wanders 
(1932) on the theory of sunspots. They investigated the wavelength de- 
pendance of the intensity ratio umbra/photosphere and also the variation of 
this ratio with the distance of the spot from the centre of the solar disc for 
the two cases in which the thermal equilibrium in the optically accessible 


regions of the sunspots is: (a) convective, and (b) radiative. The main 
results brought out were:— 


(a) Sunspot in Convective Equilibrium: 


(i) The contrast between the spot and the photosphere should increase 


(or the intensity ratio umbra/photosphere should decrease) rapidly towards 
the solar limb; 


(ii) The intensity ratio umbra/photosphere should increase gradually 
with increasing wavelength. 


(b) Sunspot in Radiative Equilibrium: 


(i) The contrast between the spot and the photosphere (or the intensity 


ratio umbra/photosphere) should remain practically constant over the entire 
solar disc. 


(ii) The intensity ratio umbra/photosphere should increase more rapidly 
with wavelength than in the case of convective equilibrium. 


Minnaert and Wanders found that the experimental results were incompatible 
with adiabatic equilibrium, but were in good agreement with the assump- 
tion that sunspots are in radiative equilibrium at an effective temperature 
of 4300° K. They concluded that the levels at which cooling by ascending 
motion takes place in sunspots are below the depth accessible for optical 
observation, a result in conformity with the fact that displacements of the 
spectral lines towards the violet have not been observed in the umbre of 
spots near the centre of the solar disc. 


Following this investigation, Wanders (1935) made an experimental study 
of the intensity ratio umbra/photosphere for a large number of spots in the 
44000 A region. He found that the ratio remains practically constant over 
the solar disc, and that larger spots show smaller central intensities. The 
latter conclusion was confirmed by Wormell (1936) who made measurements 
of the intensity of total radiation from sunspots. During the sunspot 
maximum epoch of 1938, Richardson (1929) made measurements of the 
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centre to limb variation of the intensities for a number of large stable spots 
at four different wavelengths (AA 4100, 5100, 5800 and 6600 A). His results 
appeared to support the conclusion that sunspots behave as though they 


are in radiative rather than in adiabatic equilibrium. He deduced an effective 
temperature of 4300° K. for the umbra. 


Abetti and Colacevich (1939) made a spectrophotometric study of a 
large sunspot in the blue (A 4227 A) and red (A 6568 A) regions of the spec- 
trum. They found the values 3-7°% and 12% respectively for the intensity 
ratio umbra/photosphere corresponding to an umbral temperature of 
3690° K.—the lowest recorded temperature of a spot.* Ten Bruggencate 
and Von Kliiber (1939) studied the total absorption in the Fraunhofer lines 


in the spectra of some large sunspots from which they deduced an excitation 
temperature of 3800° K. for the spots. 


Thackeray (1940) showed that if the absorption coefficient of the solar 
atmosphere is assumed to increase linearly with pressure, (Minnaert and 
Wanders considered the two cases of k constant and k « pT-) then the 
observed ‘“‘ angle variation” and “ wavelength variation”’ of the intensity 
ratio spot/photosphere can be explained equally well either as the assumption 
that both disc and spot are in convective equilibrium or that both are in 
radiative equilibrium. 

On the assumption that sunspots are in radiative equilibrium, Waldmeier 
(1941, 1942) showed that it is possible to account for the Evershed Effect 
based on a model solar atmosphere. 


2. OBJECT OF THE PRESENT WORK 


Experimental data on the “angle dependance” and “ wavelength 


dependance ” of the intensity ratio umbra/photosphere which are obviously 
of great importance for the theory of sunspots are comparatively scanty. 
The conclusions of Minnaert and Wanders as well as of Thackeray are mainly 
based on the experimental results of Pettit and Nicholson. Although the 
subsequent work of Wanders and Richardson has given support to the idea 
of radiative equilibrium for sunspots there appears to be no other mecha- 
nism excepting convective ascent of gases from below which can account for 
the cooling observed in sunspots. The collection of more quantitative 
experimental data on the spectral energy distribution for a larger number 
of sunspots appears desirable and hence it was decided to undertake a syste- 
matic programme of spectrophotometry of sunspots extending over the 
entire photographic region of the spectrum whenever comparatively large 


* Transactions of the International Astronomical Union, Voi. VU, p. 381 (1950). 


Spectrophotometric Study of Sunspots 271 


stable spots could bé studied under good sky and seeing conditions. This 
work was taken up in the beginning of 1951. Unfortunately the number of 
favourable occasions has not been many. The results of measurements 
made on some spots during 1951 are reported in this paper. 


3. EXPERIMENTAL ARRANGEMENTS 


The optical system employed in the present work for forming the solar 
image was the same as that used for taking the daily spectroheliograms at 
the Kodaikanal Observatory. It consists of an 18-inch siderostat and 
a 12-inch photovisual lens of 21 ft. focal length giving a solar image 
of about 60 mm. diameter. The spectra were photographed with an auto- 
collimating spectrograph of 20 ft. focal length employing a 6} inch Michelson 
grating. All photographs were taken in the first order, the dispersion being 
practically uniform and very nearly 3 A/mm. between A 6600 A and A 5000 A. 
In place of the grating a set of 3} prisms can be employed as the dispersing 
system in this spectrograph. Some spectra in the H-a region were taken in 
this manner, the dispersion being about 4 A/mm. 


By means of the slow motion controls of the siderostat the solar image 


. was so maintained that the vertical slit of the spectrograph passed through 


the centre of the umbra of the spot which in turn bisected the length of the 
slit. A slit width of -04 mm. was employed for photographing the spectra. 
Since the spectrograph was not capable of rotation and since the optical field 
of the siderostat gradually rotates, the orientation of the spectrograph slit 
with reference to the spot was computed in each case from the times at which 
the spectra were taken. 


Ilford Special Rapid Panchromatic plates specially sensitised for H-a 
were used for the spot spectra in that region; ordinary Ilford Special Rapid 
Panchromatic plates were used for the Na-D and Mg-b regions. Five 
spectra were taken on each plate in rapid succession, the exposure times 
varying from | to 5 seconds for the H-a and Na-D regions and from 2 to 10 
seconds for the Mg-b region. The width of each spectrum was about 6 mm., 
the central spot band covering a width of 1 mm. in the case of the larger spots. 
A set of intensity marks was also impressed on each plate using a photo- 
graphic step wedge which was prepared in the manner described in the next 
paragraph. The centre of the solar disc was used as the source of light in 
this case, the slit width of the spectrograph being so adjusted that the time 
of exposure required to get intensity marks with suitable density with the 
wedge in front of the slit being about the mean of the exposure times for the 
set of spot spectra. All plates were developed for 5 minutes in M-Q deve- 


loper at room temperature, taking the usual precautions. 
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Since the solar image was only about 60 mm. diameter, the step wedge 
had to be correspondingly small to avoid the error due to limb darkening. 
The wedges employed were prepared in the following manner. A master 
wedge containing 18 steps with suitable density gradation was prepared on an 
Ilford Process plate by exposing the plate in the blue region of the second 
order spectrum in a high dispersion grating spectrograph the slit of which 
was illumined by a Philips tungsten ribbon lamp. The different density 
steps were obtained by appropriately varying the exposure times from 10 
seconds to 150 seconds. The 18 steps and the intervening gaps together 
occupied 28mm. on the plate. Because of the narrow spectral range 
involved, there was little variation of intensity along the length of any one 
of the steps. This wedge was next placed in an enlarger-cwm-reducer and 
a large number of reduced negatives of it were made on fine-grained Ilford 
Lantern and Thin Film Half Tone plates. In this process, special care was 
taken to see that the master wedge was uniformly illumined. By suitably 
controlling the exposure time and development, negatives having various 
gradations of density could be easily prepared. The total width occupied 
by the 18 steps and the intervening spaces in the reduced wedges was 9 mm. 
While the gap between successive steps is transparent in the master wedge, 
it is opaque in the reduced wedges, the density gradation of the steps being 
in the opposite direction to that on the master wedge. Since the steps in 
the reduced wedges employed in the work extended only 4-5 mm. on either 
side of the centre of the solar disc, the error due to limb darkening was 
negligible compared with the other errors of photographic photometry 
(For sin 8 = 0-2 which corresponds to a distance of 6 mm. from the centre 
of the disc on a 60 mm. image, the intensity of the solar disc varies from 98% 
of the central intensity at A 3700 A to 99% at A 6700 A). 


The step wedges employed in the photometric work were calibrated in 


the manner described by Shane (1932). The relative transmissions of two 
wedges employed in the present work are given in Table I. 


TABLE I 

Relative Transmission Factors of the Step Wedges 

Step No:! | | 
Wedge No.| | | 
| 


| 


12 | 13] 14 | 15 | 16 | 17/18) 


(100,75+1 4:2 


| | 
19-013-8 11 8 


10-0 8-5 7-0) 5-4| 4-8 3°7 3-4 


1l 10-1) 
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4. RESULTS 


A good series of spectra were obtained on a single large spot in tho 
great sunspot group of 1951 May 9-22. This group. whose mean helio- 
graphic latitude was 14° N, crossed the central solar meridian on May 16 
and attained a maximum area of 4850 millionths of the sun’s visible hemi- 
sphere. According to H. W. Newton (1951) it is the third if not the second 
largest sunspot group on record. The successive appearances of this group 
during its passage across the disc are illustrated in Pl. X, Fig. 1. 


At Kodaikanal atmospheric conditions were good on May 12 and 13 
and a large number of spectrograms were taken on these two days in the 
H,, Na-D and Mg-b regions of the spectrum. All spectra showed excellent 
definition and fine details. A typical set of spectra in the Mg-b region is 
reproduced in Pl. XI, Fig. 2 (a). 


Spectrograms in the H, region were also taken on three other spots 
during May 1951, while the spectra of a single large spot were photographed 
on three days in October 1951. 


Particulars regarding the spots whose spectra were studied are given in 
Table Il. 


The appearances of the spots on the respective dates are illustrated 
in Pl. XII, Fig. 2 (0). 


It is well known that apart from the umbra and the penumbra which 
are the most conspicuous features of a sunspot, there is a bright ring sur- 
rounding the umbra which is slightly more intense than the photosphere. 
This is quite distinct from the facule which are generally noticed around 
spots especially when they are near the limb. The bright ring around the 
penumbra can be seen on all good photoheliograms particularly for spots 
near the central meridian when the facule are quite inconspicuous. Accord- 
ing to Waldmeier (1939) this bright ring is some 2 to 3% brighter than the 
photosphere on photoheliograms taken in violet light (~4000A). Just 
as the transition from the penumbra to the photosphere is discontinuous, 
the transition from the umbra to the penumbra is also discontinuous. 
It was noticed by Secchi (1875) long ago that the inner edge of 
the penumbra—the edge towards the periphery of the umbra—is 
brighter than its outer edge. This observation was doubted by Pettit and 
Nicholson but has been confirmed by later workers (Strebel, 1932; 
Ananthakrishnan, 1951). Visual and microphotometric examination of 
regular spots on good photoheliograms show that the inner edge of the 
penumbra is brighter in many cases, the intensity maximum at the outer edge 
of the umbra being quite conspicuous in some cases. The microphotometer 
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records of two sunspots illustrating the discontinuous nature of the transi- 
tion from the umbra to the penumbra are reproduced in Pl. XIII, Fig. 3. 
Thus, starting from the undisturbed photosphere we have first the outer 
bright ring, then the penumbra, then the inner bright ring and finally the 
umbra of the spot. In conformity with this structure of sunspots it is found 
that on many of the good spectrograms of spots the umbral and penumbral 
spectra are bordered by bright margins along their entire length. 


Each of the spectrum plates covered a range of about 300A. They 
were microphotometered perpendicular to the direction of dispersion at 
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selected wavelengths in the continuous spectrum using the Cambridge photo- 
electric microphotometer of the observatory. All records were taken under 
a magnification of x50, care being taken to avoid the absorption lines. A 
typical record is reproduced in Pl. XIV, Fig. 4 (a). The density variations 
were converted into intensity variations using the characteristic curves 
worked out from the intensity marks Fig. 4(b). Figs. 5 (a), (b), (c) and (d) 
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illustrate typical curves of the intensity variations across one of the sunspots 
studied. The abscissa for these curves is seconds of are (geocentric) and 
the ordinate relative intensities. (The intensity of the photosphere suffici- 
ently far away from, the spot is taken as 


an 

I 

100 . 

90 

80 

70 

50 } 

40 

30 

10 

I 

100 

- 

g0 

10 

GO 

50 

40 

30 

20 

10 

20 30 40° 


278 R. ANANTHAKRISHNAN 


Table III gives the relative intensities of the umbra and the penumbra 


for the various spots studied. The penumbral intensities are the average of 
the values on either side of the umbra. 


TABLE III 


‘ Rel. Int. (%) | Eff. Temp. (° K) 
Spot No. : 


KKL 9620 


5460 
(13-5-1951) 6540 


5440 
5450 


who 
iste 
= 


5450 
5871 23-8 80- 4300 5460 
5430 


5399 
5170 14-4 70-2 4100 | 5350 
q 7 5390 


KKL 9687 27: 76-5 
(19-35-1951) 68 | 5 


KKL 9688 28- 
(19-5-1951) 


KKL 9692 
(24-5-1951) 


KKL 9768 
(14-10-1951) 


Do +4 
(16-10-1951) : 6540 26-2 


Do --| 6375 | 22 
(18-10-1951) -+| 6620 29- 


U = umbra; P = penumbra 


Because of the rapid fall in the sensitivity of Ilford special rapid pan- 
chromatic H-a plates on the longer wavelength side of H, it was found 
that the same spectrum could not be utilised for photometry for wavelengths 
on the red and violet sides the H, line. Spectra with longer exposures 
(about 3 to 4 seconds) had to be used for the red side of Hy. This fact would, 
to some extent, account for the sudden change in the intensity ratio for a 
comparatively small change of wavelength shown by some spots listed in 
Table III. Nevertheless, it is clear that the spot of May 24 had a smaller 
umbral intensity as compared with the other spots, 
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6620 30-3 85-8 4380 5510 
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The major sources of error in the spectrophotometry of sunspots arise 
from the scattering by the earth’s atmosphere and by the instrumental optics 
as well as from “ poor seeing’’. The result of all these is to throw photo- 
spheric light into the spot and thereby enhance the measured intensity ratio 
umbra/photosphere. The manner in which the measured values may be 
corrected has been theoretically discussed at some length by Wanders (1934). 
The method is based on the study of the observed curve of intensity fall just 
outside the solar limb photographed simultaneously with the sunspot spectrum 
under identical sky and instrumental conditions. However, the practical 
application of the formula worked out by Wanders is far from simple and 
he himself has made use of a procedure originally employed by Pettit and 
Nicholson to correct his measurements. These authors measured the ratio 
of the radiation from the sky and instrument 15” of arc outside the solar 
limb in terms of the radiation 15” of arc within the limb. If this ratio is 
x then they assumed that the measured energy ratio umbra/photosphere 
has to be diminished by 2x. Wormell corrected his measurements by 
making some simplifications in the formula given by Wanders to render 
its practical applicability more easy. 


In the present work the spectra of the limb and adjacent sky were photo- 
graphed on some days almost simultaneously with the spot spectra to find 
out the magnitude of the correction for the measured umbral intensities. 
An attempt was made to evaluate this correction on the lines of Wormell. 
However, the method was found to be rather uncertain and in some respects 
arbitrary. Hence the attempt was given up. All the umbral and penumbral 
intensities given in Table III are only the values as obtained by direct 
measurements on the plates. Almost all the spot spectra were photographed 
between 8 and 10 A.M. when the definition of the solar image is comparatively 
good at Kodaikanal; sky conditions were generally good on all the days 
on which the spectra were photographed. Hence the corrections required 
for the measured umbral intensities are expected to be small. 


5. EFFECTIVE TEMPERATURES OF SPOTS 


According to the measurements of Wanders and Wormell, the intensity 
ratio umbra/photosphere varies within wide limits depending upon the size 
of the spot. Wanders found values of umbral intensities ranging from 5% 
to 40% of the photospheric intensity in the A4000A region, while in 
integrated light Wormell found values ranging from 28% to over 70%. The 
values in Table III also show some variations, but the range of variation is 
smaller. This is probably because all the spots were comparatively large 
and the sky and seeing conditions generally good on all the days. 
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Assuming T = 5740° K. for the effective temperature of the photo- 
sphere, the values of the ratio I (A, T*)/I (A, T) for various values of \ and T* 
computed according to Planck’s radiation formula are given in Table IV 
and represented graphically in Fig. 6, for the temperature range 3000° K. 
to 8000° K. and for the wavelength range A 3000 A to A 10000 A. 


The range of values for the intensity ratio umbra/photosphere found by 
Wanders corresponds to effective spot temperatures from a little under 
4000° K. to a little over 5000° K. For total radiation the values of Wormell 
give effective spot temperatures ranging from 4175° K. to 5250°K. The 
effective temperatures for the spots studied in the present work are given in 
Table III. 


6. ORIGIN OF THE “‘ BRIGHT RINGS” IN SUNSPOTS 


Odgers (1946) has given a theory of sunspots which indicates the prob- 
able origin of the bright rings seen around the umbra and the penumbra. 
He regards the sunspot as a region where the product of the mass absorption 
coefficient and the density is higher than that in the photosphere. It is well 
known that this product is a measure of the “ obstructive power” of the 
material of the stellar atmosphere to the passage of radiation through it. 
Hence the assumption of Odgers implies that radiation encounters greater 
obstruction to its passage through the sunspot region than through the sur- 
rounding photosphere. The mathematical problem of evaluating the 
energy density and the flux of radiant energy inside and outside the sunspot 
region assuming radiative equilibrium in both cases becomes analogous to 
the classical electrostatic problem of the effect of introducing a dielectric 
in a uniform electric field.* Three cases have been considered by Odgers, 
viz., (i) when the sunspot can be idealised by a sphere of radius a and 
absorption coefficient k,/p, embedded in the homogeneous radiation field of 
the surrounding photosphere whose absorption coefficient is k/p; (ii) when 
the sunspot has the shape of a prolate ellipsoid; (iii) when the sunspot has 
the shape of an oblate ellipsoid. In all cases the problem reduces to the 
solution of Laplace’s equation with the appropriate boundary conditions. 


For the first case referred to above, the solutions are:— 


= Fx[ (2) 


* Sir James Jeans, The Mathematical Theory of Electricity and Magnetism, Fifth Ecn., 
p. 228 (1951). 
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where E, and EB, are the energy densities outside and inside the spot, 7 is the 
distance from the centre of the sunspot sphere which is taken as the origin 
of co-ordinates, x is the distance measured in the direction of the flux, and A 
is a constant. 


Sufficiently away from spot (r = oo) the energy density in the photo- 
sphere is given by: 


E,, = Fx +A 


If F; is the energy flux inside the spot we have: 


F; k, = F | = constant. 


ox 


In the yz plane the energy flux outside the spot is given by: 


(3) 


From (3) and (4) we see that the ratio of the flux just outside the spot (r=a) 
to that inside is given by k,/k. From the observed intensity ratio 
umbra/photosphere we can calculate k,/k according to equation (3). The 
manner in which the energy flux outside the spot falls off with distance from 
the centre of the spot for various values of the intensity ratio umbra/photo- 
sphere is shown in Table V. 


The flux is a maximum just outside the spot and falls off at first rapidly 
and then slowly till it merges with the value for the undisturbed photosphere. 
Physically this means that the flux held back by the increased obstructive 
power of the spot region reappears as an excess of flux immediately surround- 
ing the spot analogous to the flow round a sphere immersed in a fluid moving 
with uniform velocity. The increased energy flux around the spot appears 
as the “ bright ring ”’. 


Odgers has also discussed the case of a spherical spot with umbra and 
penumbra. In this case three different values of k corresponding to the 
umbra, penumbra and the photosphere have to be considered. The ana- 
lysis shows that there should be two bright rings one around the umbra and 
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the other around the penumbra. The consideration of spots which resemble 
prolate or oblate ellipsoids shows that the numerical value of the sudden 
increase in flux just outside the spot depends on the shape of the spot. For 
a spot whose depth is great compared with the aperture the increase of flux 
just outside the spot is smaller than for a shallow spot. 


According to Odgers’ theory the variation in the intensity ratio 
umbra/photosphere as well as the variation in the depth/aperture ratio of 
spots can apparently account, at least qualitatively, for the observed variation 
in the intensity of the bright rings. There is no other theory which attempts 
to explain this almost universal phenomenon associated with sunspots. 
There is, however, an important point which has to be remembered in this 
connection. It is generally believed that the depth of a sunspot is much 
smaller compared with its aperture. In the case of medium and large spots 
the aperture is of the order of 10*km., while according to Unsdld’s theory, 
for instance, the depth will only be of the order of 10?km. In such a case, 
the intensity of the bright ring should be quite appreciable according to 
Odgers’ theory. However, the observed intensities are generally small 
judged from visual and microphotometric examination of sunspots. For 
the spots studied in the present work the inner bright ring was about 2% 
more intense than the penumbra in one case and less than this in the other 
cases. The brightness of the outer ring was also estimated to be the same 
order. Interpreted in terms of Odgers’ theory this would imply that spots 
extend to appreciable depths below the photosphere. 


SUMMARY 


After a brief review of some of the earlier work on sunspots the experi- 
mental technique employed for the spectrophotometry of sunspots and the 
results obtained in the case of 5 spots are described. In the case of one of 
the spots studied on 1951 May 13, the intensity ratio umbra/photosphere 
decreased from about 32% at A6620A to 15% at 15078 A. The lowest 
value of 22% at 6620 A was recorded for the umbral intensity of a spot 
on May 24. Microphotometer records of two spots are reproduced showing 
the discontinuous nature of the transition from the umbra to the penumbra 
first observed visually by Secchi. Tables and curves of I (A, T*)/I(A, T) 
based on Planck’s radiation law are given and the effective temperatures of 
the spots calculated. Odgers’ theory of the “ bright rings” in sunspots is 
considered and it is pointed out that the observed small values of excess 


intensity of the bright rings probably imply that sunspots extend to appre- 
ciable depths below the photosphere. 
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METEOROLOGICAL AND EXTRA-TERRESTRIAL 
CAUSES OF THE DAILY VARIATION OF 
COSMIC RAY INTENSITY 


By V. SARABHAI, U. D. DesAlI AND R. P. KANE 
(Physical Research Laboratory, Ahmedabad, India) 


’ Received January 27, 1953 


1. INTRODUCTION 


It has long been realised that a study of the solar and sidereal time daily 
variations of the cosmic ray intensity could give some clue to the location 
of regions where cosmic rays originate. However, a successful attempt to 
follow up this line of thought has not so far been possible due to the un- 
certainty in separating the variations caused by terrestrial influences from 
those due to an anisotropy of the primary radiation. 


In investigating the time variations of cosmic rays, the ionisation chamber 
offers the great advantage of constancy of operation, but being an omni- 
directional detector of radiation, it is hardly a satisfactory instrument for 
the study of an anisotropy of the primary radiation. Nevertheless, very 
valuable data have been collected with it. Apart from the Carnegie Insti- 
tution studies made at widely separated places on the earth and reported 
by Lange and Forbush,! observations have been made, amongst others, by 
Hess and Grazeadei? at the Hafelekar, by Schonland, er al.,> at Capetown 
and by Hogg* at Canberra. Unidirectional measurements of the diurnal 
variation of the vertical meson intensity, performed with narrow angle 
geiger-counter telescopes could be more revealing than omni-directional 
measurements. But the only extensive data with vertically pointing tele- 
scopes comes from Duperier® in whose experiment the angles of the tele- 
scopes were fairly wide. 


Interpretation of the ionisation chamber and Duperier’s experiments 
has been much confused by various differing corrections for meteorological 
factors that have been applied, and it has not been possible finally to deter- 
mine how much of the diurnal variation is due to an anisotropy of the 
primaries. To overcome this difficulty Alfven and Malmfors* and Elliot 
and Dolbear?»* have studied the daily variation of cosmic ray intensity with 
telescopes pointing in the North and South directions. While the North- 
South daily variation difference curve is substantially independent of atmo- 
spheric effects and constitutes evidence for an anisotropy of theprimary radia- 
tion, it is difficult to interpret it further. For, an anisotropy of primaries, 
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such as may be caused by solar emission of cosmic rays, can produce a 
daily variation in both North and South pointing telescopes. The differ- 
ence curve in consequence reflects an arithmetic difference between the daily 
variations in the two directions due to anisotropic cosmic ray primaries, 
but does not reveal the true nature of the daily variation due to primary 
anisotropy in either direction. 


It is felt that a satisfactory solution to the problem must begin with an 
understanding of the nature of terrestrial effects on the solar daily variation 
of cosmic ray intensity. These have then to be corrected for, leaving a 
residual daily variation essentially of extra-terrestrial origin. The daily 
variation of meteorological elements is more pronounced and regular at 
places in low latitudes than at high latitudes. It is particularly appro- 
priate therefore to study the daily variation near the equator. With this in 
view, apparatus has been designed to carry out comparable studies of the 
daily variation of the total intensity as well as the meson intensity at 
Ahmedabad (Mag. Lat. 13° N., Alt. 50 metres) and at Kodaikanal (Mag, 
Lat. 1° N., Alt. 2,340 metres). As it is important to study the intensity of 
particles incident in a narrow cone in a fixed direction, compromise has to be 
made in the design of the apparatus to make the angles of the telescope 
narrow and still to retain an adequate counting rate for good statistics. We 
shall describe in this paper details of the apparatus and the results obtained 
during the past two years at Ahmedabad. A summary of the results obtained 
at both stations has already been communicated elsewhere.® Details of 
the Kodaikanal results will be presented later as soon as more significant 
data are available. We discuss here reasons which lead us to believe that 
the solar time daily variation of meson intensity corrected for barometric 
pressure is caused by an anisotropy of the primary radiation, probably con- 
nected with the emission of charged particles from the sun. 


2. THe APPARATUS 


A schematic diagram of the apparatus is given in Fig. 1. 


Five trays, each with four ‘self-quenched geiger counters connected in 
parallel, form three triple coincidence vertically pointing telescopes of identi- 
cal dimensions. The counters have copper cathodes 30cm. long with 
diameter of 4cm. The counters are placed in the N-S direction, and each 
telescope subtends a semi-angle of 22° in the E-W plane and a semi-angle 
of 37° in the N-S plane. Since the purpose of the experiment is to measure 
the daily variation of cosmic ray intensity connected with the rotation of 
the earth, the apparatus is oriented so that the telescopes present the smaller 
angle in the E-W plane. 


Causes of the Daily Variation of Cosmic Ray Intensity 


COOO OOOO ¢ | 


Fic. 1. Schematic diagram of apparatus showing counter trays A, B,C, D, E, A’,B’, C’, D’ 
and E’, quenching units Q, triple coincidence units N and scale of four units S feeding the 
electro-mechanical recorders R. 


Lead absorbers are placed between the 3rd and 4th trays and the 4th 
and Sth trays. While therefore, the uppermost telescope ABC measures 
the total cosmic ray intensity T, the lower two telescopes BCD and CDE 
measure intensity that can penetrate through 7cm. and 17cm. of lead 
respectively. The soft component E is almost completely eliminated by 
7cem. of lead and is given by the difference between counting rates of tele- 
scopes ABC and BCD. The penetrating component ‘ m’ consisting mostly 
of w-mesons is measured by the telescope CDE. The small difference 
between the counting rates of telescopes BCD and CDE represents the 
intensity of an intermediate component ‘I’ consisting partly of the very 
energetic soft component and partly of the slow meson component. Since 
‘I’ is negligible compared to ‘m’, the counting rates of telescopes BCD 
and CDE may be considered together to represent the meson component M, 
when distinction is being made only with the electronic component E. 


The complete apparatus has three duplicate telescopes A’B’C’, B’C’D’ 
and C’D’E’ operating alongside the ones described above. These are also 
shown in Fig. 1. The object of providing these additional telescopes, 
measuring identical components of the cosmic ray intensity, is not only to 
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improve the statistics but also to provide for continuance of data during 


periods when faults develop in counters or circuits connected with one or 
other of the telescopes. 


All geiger counters are placed in a heat insulated box whose temperature 
is thermostatically regulated to 105° + 2° F. Each counter tray is connected 
to an external electronic quenching unit which, for every discharge of one 
of the counters, feeds to their central wire a square negative voltage pulse 
of about 300 volts and 800 micro-seconds duration with a very sharp leading 
edge. These quenching units improve the flatness of plateau and prolong 
the life of self-quenched counters considerably by suppression of multiple 
discharges. This is important in time variation experiments where reliable 
operation over long periods of time is essential. The low impedence cathode- 
follower outputs of the quenching units are fed to fast triple coincidence 
units. Finally, the coincidences are scaled by a factor of 4 or 8 and recorded 
on telephone call registers which are automatically photographed hourly 


on standard 35 mm. film. All power supplies are electronically regulated 
to ensure stability of operation. 


Hourly values of atmospheric pressure and temperature are obtained 
from daily charts of an accurate micro-barograph and a thermograph. Upper 
air meteorological data are obtained from radiosonde ascents, with I.M.D. 
F-type or the Vaisala type instruments, conducted by the atmospheric 
physics division of the laboratory under Prof. K. R. Ramanathan. Details 
of these experiments will be published elsewhere. 


3. ANALYSIS OF DAILY VARIATION DATA 


Even though the primary data for the intensity of cosmic ray compo- 
nents and the surface atmospheric pressure and temperature are available 
for hourly intervals, the analysis has been done for bihourly intervals com- 
mencing from midnight Indian Standard Time, which is 40 minutes in 
advance of the local time at Ahmedabad. The criterion used for elimina- 
tion or inclusion of data for any particular day is the range of bihourly devia- 
tions. The data is discarded for days on which any individual bihourly 
' value is more than 5% different from the mean for the particular 

day. This corresponds to a deviation exceeding three times the expected 
standard deviation for a bihourly value. Such cases are generally attri- 
butable to some faults either in the electronic circuits or in the counters. 
Cascs of abnormally large daily variation in cosmic rays which may rarely 
occur are however also rejected on this criterion. The useful data, as 
presented here, extend from May 1950 to September 1952 and include about 
(OU days with a fairly cven distribution over the four seasons. 
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- The annual mean daily variation given by the bihourly percentage 
deviations from mean of the total intensity T, the meson intensity M and 
the electron intensity E are sown in Fig. 2. The bihourly deviations from 


08 
Hours (1.S.T.) 

Fic. 2. The smoothened daily variation of total cosmic ray intensity T, mesons M, 
electrons E, barometric pressure P and surface atmospheric temperature @. No correction 
has been applied to the cosmic ray values and the solid lines showing the daily variations 
are formed by the superposition of the diurnal and the semidiurnal components for each 
intensity. 


mean of the atmospheric pressure P and the surface temperature @ are also 
shown. For the purpose of smoothening the data, moving averages over 
three consecutive bihourly intervals have been taken for all variates shown 
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in the figure. The standard deviations for each bihourly value of cosmic 
ray intensity are also indicated. 


The amplitudes and hours of maxima of the first four harmonic compo- 
nents of the unsmoothened daily variation of T, M, E as well as P and @ are 
indicated in Table I. The hour of maximum is expressed in terms of the 
angl: in the harmonic dial representation between midnight and the vector 
for the particular harmonic component of the daily variation. 


TABLE I 


Amplitudes and hours of maxima of harmonic components of 
the daily variation of cosmic ray intensities and atmospheric 
pressure and temperature 


lst Harmonic 2nd Harmonic 3rd Harmonic 
24 hourly 12 hourly 8 hourly 


Variate 


Ampl. Max, Ampl. Max, Ampl. Max, Ampl. 


38° 


M 144° 


E 


P mm. of Hg 135° 


6-60 54° 1-60 114° 


It will be observed that for all variates the predominant harmonic 
component is either the first or the second one. In mesons and in atmospheric 
pressure, the two components are about equally important. For total in- 
tensity and for electrons, the second harmonic is larger than the first; but 
for surface temperature, the converse is true. In all cases, the third and 
higher harmonics are small and may be neglected. In what follows there- 
fore, only the 24 hourly diurnal and the 12 hourly semidiurnal components 
are taken into consideration. In discussing these components, it is con- 
venient to adopt a notation explained in Table II. 


4. INFLUENCE OF METEOROLOGICAL FACTORS ON THE DAILY 
VARIATION OF MESON AND ELECTRON INTENSITIES 


Extensive studies have been made to relate the day-to-day variations 
of the cosmic ray intensity with meteorological changes in the atmosphere. 
Duperier!® has shown that changes of meson intensity are connected with 
a mass absorption effect, an effect due to alteration of the probability of 


4th Harmonic 
6 hourly 
| | 
T -09% 29° 
108° 
—61° 1-40% 112° = +30% 76° 
-00 
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TABLE II 


Notation used for describing the amplitude and the hour of maximum 
of harmonic components 


% amplitude of meson (M), diurnal variation (D) observed at Ahmedabad (A) 


angle corresponding to hour of maximum of the meson (M) diurnal varia- 
tion (D) observed at Ahmedabad 


ms = % amplitude of meson (M) semi-diurnal variation (S) observed at Ahmeda- 
bad (A) 


M@, = angle corresponding to hour of maximum of the meson (M) semidiurnal varia- 
tion (S) observed at Ahmedabad 


Similarly T?, a, pp and ed represent the amplitudes of the diurnal 
variations at Ahmedabad of T, E, P and @ respectively 


meson decay accompanying changes of heights of isobaric levels and an effect 
of the temperature or density of the atmosphere near the 100 mb. level. The 
physical processes responsible for the positive upper air temperature effect 
are not clearly understood but it should not be expected that in the daily 
variation the influence of these three factors on meson intensity would be 
identical to what is found for day-to-day variations. This is because in 
barometric pressure as well as in atmospheric temperature, the day-to-day 
changes are brought about under very different circumstances from those 
that produce the daily variations. 


Processes responsible for day-to-day changes of barometric pressure are 
entirely different from those causing the dynamical periodic oscillations of 
the barometric pressure. The use of a barometric coefficient obtained from 
studies of day-to-day variations for correcting cosmic ray daily variation 
data in respect of the daily variation of pressure is therefore questionable. 
But this is exactly what has been done by most authors in the past. A better 
method appears to be to derive a barometric coefficient from daily variation 
studies for subsequent application to the same data. In doing this, we have 
to keep in mind available knowledge on the physical processes responsible 
for the daily variation of the meteorological elements and the special features 
of the atmospheric oscillation. 


It is difficult to draw conclusions about the effect of meteorological 
factors on cosmic ray intensities by a comparison of the daily variation curves 
of Fig. 2. Solar radiation and gravitational forces are the most important 
causes of the daily variations observed in geophysical elements. These 
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variations, as well as one that could be caused in meson intensity by an 
anisotropy of the primaries due to solar emission of cosmic rays, would have 
a predominant 24 hourly diurnal component. Therefore, it is not clear 
how much of the M® variation is connected with PP and 6° or a hypothetical 


upper air diurnal temperature variation, and how much is due to a solar 
cosmic component. 


In the 12 hourly semidiurnal components however, the position is differ- 

The atmospheric pressure, unlike temperature, has a very appreciable 
PS component. At Ahmedabad, P§ is as important as PR. At Kodaikanal, 
nearer the equator and at a higher level, P§ is 6 times PR. As was originally 
pointed out by Kelvin,?* the semidiurnal variation of pressure is due to reso- 
nance in the atmosphere which has a free period of oscillation of nearly 
12 hours. Thus, even though the exciting solar force is diurnal, the semi- 
diurnal component in pressure becomes important and is predominant at 
low latitudes. If attention is therefore confined only to the semidiurnal 
components, we have a means of studying the influence of pressure, un- 
contaminated with effects due to temperature variations in the atmosphere 
or due to an anisotropy of cosmic ray primaries. 


ent. 


The first two harmonic components of T, M and E along with those 
of P and @ are shown in Fig. 3 on 24 and 12 hourly harmonic dials. It will 


~< 


\ oolHr, / 7 7 


P| 
mm. 06 


4°C. 


FIRST | HARMONIC SECOND | HARMONIC 


Fic. 3. Diurnal and semidiurnal harmonic dials showing amplitudes 
maxima of uncorrected variations, 


and hours of 


b 
al 


4°C. 
0-5mm. 
9 
Hr. | 
b mum 
M 
E 
: 


Causes of the Daily Variation of Cosmic Ray Intensity 295 


be observed that while on the 24 hour dial, the vectors lie all around the clock, 
there is on the 12 hour dial a very striking grouping of the cosmic ray vectors 
almost completely opposite in phase with the pressure vector. 


Correlation analysis of the semidiurnal vectors of cosmic ray compo- 


nents with atmospheric pressure give correlation coefficients and barometric 
coefficients as shown in Table III. 


TABLE III 


Correlations with barometric pressure and the barometric coefficients 
of the semidiurnal components of cosmic ray intensity 


| 


Cosmic ray Correlation coefficient | eee 

companana with ps | Barometric coefficient 
= — -96 Py = — 4°2% percm. Hg 
| = — -98 Bu = — 2°4% per cm. Hg 
E§ | ree = — -95 Bs = — 14°0% per cm. Hg 


The barometric coefficient By, = — 2°4% per cm. Hg for mesons may 
be compared with the value found by other workers from day-to-day varia- 
tions of cosmic ray meson intensity. Coefficients of — 3:0% per cm. Hg 
for Huancayo and —1-8% per cm. Hg for Cheltenham, Christchurch 
and Godhavn have been determined and used for barometric pressure correc- 
tion by Lange and Forbush! for the Carnegie Institution ionisation chamber 
data. It is not clear why at Huancayo the coefficient should be so much 
larger than at the other stations in spite of the shielding being the same for 
all instruments. Duperier’s!° barometric coefficient is —1-50% per cm. Hg and 
Dolbear and Elliot!! have reported a value of — 1-88% per cm. Hg obtained 
from the seasonal variation of intensity. These authors, by partial correla- 
tion analysis, give estimates for the three meteorological coefficients which 
affect meson intensity. These are shown in Table IV. Our value of the 
barometric coefficient for mesons is larger than the true mass absorption 


coefficient of Duperier but agrees well with the coefficient of Dolbear and 
Elliot. 


It is important to examine if there is a substantial decay contribution 
in the semidiurnal variation. Pekeris,:? and more lately Wilkes and Weekes!* 
have examined the details of the modes of oscillation of the atmosphere. 
Nicholson and Sarabhai! have estimated the effect on meson intensity of 
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TABLE IV 


Estimates of barometric coefficient 8, true absorption coefficient , 
true decay coefficient »' and positive temperature coefficient a for mesons 


Coefficient Duperier?? Dolbear and Elliot!? 


— 1-50% /cm. Hg* 


— 1-88%/cm. Hg 


— 1:05%/cm. Hg 
— 3-90%/km. 


— 2:07%/cm. Hg 
— 422% /km. 


+ 0°12%/°C. + 0:14%/°C, 


* Weighted mean for 5 periods of observations. 


the semidiurnal change of height of isobaric levels due to atmospheric oscilla- 
tions. For meson production near 16 km., there should be a semidiurnal 
oscillation of the isobaric level which would not exceed 4 metres and thus 
would not change significantly the contribution of the pure mass absorption 
effect to the barometric coefficient. There is reason to believe therefore 
that the barometric coefficient derived by us from the semidiurnal variation 


corresponds mainly to the true absorption coefficient for the meson compo- 
nent. 


The appropriate barometric coefficients experimentally determined from 
semidiurnal components can be used to correct the unresolved daily varia- 
tions of T, M and E for the barometric daily variation. Smoothened 
bihourly values of the barometric pressure corrected variations designated 
as T’, M’ and E’ are shown in Fig. 4. The harmonic components of these 
are indicated in Fig. 5. 


It will be noticed that T’ and M’ are left with a residual diurnal variation 
of amplitude -33 + -05% and -35 + -04% respectively and hour of maxi- 
mum near 0900 hours I.S.T. E’ on the other hand has no significant varia- 
tion exceeding the standard deviation of the individual bihourly points. 
This indicates that the daily variation of the electron component at sea level 
can be explained almost completely by a mass absorption effect connected 
with the barometric variation. 


Duperier has corrected the daily variation of mesons for a decay effect 
due to an estimated diurnal change of height of about 50 metres in the iso- 
baric levels near 16 km. in consequence of a diurnal heating of the atmo- 
sphere. The process has been considered to be analogous to the seasonal 
variation of meson intensity where, during summer, the general expansion 
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Fic. 4. The daily variation of total cosmic ray intensity, mesons and electrons, each 
corrected for barometric pressure. 
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of the atmosphere produces a decrease of meson intensity. Dolbear and 
Elliot have suggested a further correction for a diurnal positive tempe- 
rature effect similar to the one demonstrated by Duperier!! for day-to-day 
changes. The contribution of both these factors to the daily variation of 
mesons would depend on the magnitude of the diurnal temperature change 
at different levels in the atmosphere and particularly near the tropopause. 


All available meteorological evidence goes to show that the daily varia- 
tion of air temperature gets rapidly attenuated going upwards from ground, 
and becomes negligible beyond 2 km. above the surface of the earth. It 
may again become important in the ozone ultraviolet absorption region, 
well above the tropical tropopause. Duperier,!° and Elliot and Dolbear® 
have based their arguments on radiosonde data obtained by the Meteoro- 
logical Office, London, at Larkhill (100 km. S-W of London) and Downham 
Market (75 km. N-E of London). The interpretation of the data is however 
very suspect and Kay'® has critically examined it to come to the conclusion 
that much of the apparent diurnal variation of temperature in the lower 
atmosphere is spurious and caused by inadequate radiation shielding of the 
thermal element. Flights have been made in India to test whether there is 
a diurnal variation in low latitudes. It has been tentatively found in ascents 
made with Vaisala instruments at Ahmedabad that there is no significant 
diurnal increase of temperature near the tropopause. There is ozone above 
100 mb. but the heating effect due to absorption of solar radiation is consi- 
dered to be of significance only at much higher levels. 


An examination of Fig. 5 reveals that the barometric pressure corrected 
vectors for total intensity and mesons on the 24 hourly dial are both signi- 
ficant and are negligibly correlated with the surface atmospheric temperature. 
If the heating in the upper atmosphere were to take place from lower levels, 
the maximum temperature would occur at a later hour than the surface 
temperature and the correlation with T’ and M’ would be almost zero. For 
heating of the layers of air near 16 km. from above, the maximum tempe- 
rature may occur nearer noon, but even so the correlation between the diurnal 
vectors for T’ and M’ and a temperature vector at noon would be quite low. 


One is finally led to conclude that the atmospheric temperature has little 


or no part to play in producing the daily variation of meson intensity corrected 
suitably for barometric pressure. 


5. THE BAROMETRIC PRESSURE CORRECTED MESON DrURNAL VARIATION: 


Having eliminated meteorological effects, it is necessary to consider 
possible geomagnetic and helio-magnetic influences on the residual daily 
variation of mesons. Janossy'’ has suggested the possibility of a diurnal 
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variation of cosmic rays at latitudes beyond 40° due to the helio-magnetic 
field. Dwight’® has worked out detailed implications, but this theory can 
be safely excluded in view, amongst other things, of the evidence from several 
quarters concerning the non-existence of an appreciable general helio- 
magnetic field at the present time. Vallarta and Godart'® have discussed 
the influence in low latitudes of ionospheric current systems responsible for 
the geomagnetic diurnal variation. While the latter alters fundamentally 
in character with latitude, Thompson?® has shown that the meson variation 
has similar features over a wide range of latitudes. It appears therefore that 
the barometric pressure corrected variation of mesons is of extra-terrestrial 
origin and is connected with an anisotropy of the primary radiation. 
Taking into consideration the occurrence of the maximum of the diurnal 
variation at about the same period of the day according to local time at 
widely separated places on the earth, it is reasonable to conclude that the 
anisotropy is caused by the solar emission of cosmic rays. 


Duperier® has indeed made a similar suggestion by a consideration of 
the seasonal change of amplitude of the meson diurnal variation corrected 
for barometric pressure and decay coefficient. As however, the application 
of the decay effect is questionable for reasons mentioned earlier, the close 
agreement between the ratio of summer and winter diurnal amplitudes with 


what would be expected due to change of the solar zenith distance at the two 
periods may be fortuitous. 


6. THe EFFECTS OF CosMIC RAYS FROM THE SUN 


The sun is known to emit corpuscular streams which take about 23 
hours to travel to the earth and produce geomagnetic and auroral activity. 
It is also known to emit during some intense solar flares, moderate and low 
energy cosmic ray particles which travel with a velocity not very different 
from that of light and produce measurable effects at sea level on the cosmic 
ray neutron and charged particle intensity. The magnitude of the effect 
has a marked longitude dependence, and no effects have been observed at 
equatorial stations. Increases in neutron intensity reported by Simpson, 
et al.,”" and charged particle intensity reported by Neher and Forbush2* have 
been associated with the central meridian passage of active regions on the 
solar disc. These demonstrate the emission of more energetic particles from 
the sun which make their effects felt even at Huancayo on the geomagnetic 
equator. The present association of the meson diurnal variation with an 
anisotropy due to solar cosmic rays shows that the sun is a continuous 
emitter of cosmic radiation. Unlike the bursts of radiation following the 
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observation of flares, this continuous emission is energetic enough to cause 
measurable effects in the charged particle intensity at sea level at all latitudes. 


There is some evidence to show that the energy distribution of cosmic 
ray particles from the sun is displaced towards low energies as compared 
to the general energy distribution of cosmic rays from all other sources, 
For, the percentage amplitude of the diurnal variation increases with eleva- 
tion as revealed by our comparative studies at Ahmedabad and Kodaikanal, 
by the Carnegie studies at Huancayo compared to the low level stations, 
and the studies made at the Hafelekar. Neher and Forbush have reached 
a similar view from the increase of the worldwide fluctuations of ionisation 
with altitude, and the fluctuations being less pronounced at the equator. 


An important question arises about the observed hour of maximum 
M@? of the meson diurnal variation. For high energy particles from the 
sun which are not appreciably deflected by the geomagnetic field, one should 
expect the maximum to occur at noon local time. For less energetic posi- 
tive primaries the maximum will be shifted to earlier hours, and for less 
energetic negative primaries to later hours. There is a divergence amongst 
the reported results of various workers concerning the precise hour of 
maximum. It varies in extreme cases from 0800 hours to 1600 hours. A 
great deal of this divergence is perhaps due to differences in methods of 
correcting for meteorological effects. 


Hogg** has compared on a harmonic dial the diurnal vectors for meson 
dialy variation observed by various workers at different places. For 
Canberra data, a vector has been given for barometric pressure corrected 
meson variation as well as for one which has, in addition, been corrected for 
a temperature effect. There is a considerable difference between the ampli- 
tude and the hour of maximum of the resultant variation in the two cases. 
Forbush,™ for Cheltenham data, has shown how the uncorrected meson 
diurnal vector at 1400 hours shifts to 1100 or 1000 hours according to the 
magnitude of the barometric coefficient which is applied for correction. 


Some of the differences in the hours of maxima and amplitudes observed 
by various workers are probably connected with the nature of the measuring 
apparatus and the angle within which it allows incident radiation. Gene- 
rally an omni-directional instrument would reveal a smaller amplitude of 
variation than a unidirectional one. In latitudes where there is an E-W 
asymmetry of the cosmic ray intensity, an ionisation chamber would effectively 
function like a West pointing telescope having a later maximum than a vertical 
telescope. In view of all these factors, comparisons between the diurnal 
variation at different latitudes and elevations can only be made where similar 
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experimental technique is followed at the various stations, and appropriate 
similar corrections are applied to the basic experimental data. Carnegie 
Institution studies are therefore very valuable for this purpose, and when 
further significant data is available from our unidirectional studies at 
Kodaikanal and Trivandrum (Mag. Lat. 1° N., Alt. 0 metres) it might be 
possible to get a better insight in this subject. From our own studies, 
there is every indication that the maximum occurs before noon, and the hour 
becomes earlier when the diurnal amplitude increases in going from 
Ahmedabad to Kodaikanal. Though all sea-level stations run by the 
Carnegie Institution have maxima in the early afternoon, the mountain 
station of Huancayo has an earlier maximum before noon and of increased 
amplitude. 


An important point that arises now concerns the relationship that can 
exist between the amplitude and the hour of maximum of the diurnal varia- 
tion. The amplitude should be controlled, among other things, by the cut 
off in the solar cosmic ray energy spectrum either by geomagnetic blocking 
or atmospheric absorption. The mean energy of the allowed radiation 
determines the bending in the geomagnetic field and hence the hour of 
maximum of the diurnal amplitude. When changes in amplitude of the 
diurnal variation are due to alteration of the cut off energy, a change in the 
hour of maximum may also be expected. This consideration requires a 
revision of the past analysis made for detecting a sidereal time variation in 
cosmic rays from the seasonal change of the diurnal variation on lines sug- 
gested by Thompson.” 


Due to bending of the trajectories of cosmic rays from the sun, there 
is every reason to expect a diurnal variation of intensity in both the North 
and the South pointing telescopes of Alfven and Malmfors,® and Elliot and 
Dolbear.8 On account of the asymmetry of the geomagnetic field with 
respect to the earth’s axis of rotation, the two telescopes may exhibit 
with respect to each other, a small shift of phase of the measured diurnal 
variation. Just as an East-pointing telescope would have an earlier maximum 
than a West-pointing telescope, a North pointing telescope in England 
should have an earlier maximum than one pointing to the South. This is what 
Dolbear and Elliot have found. While qualitatively the explanation is 
attractive, it remains to be seen whether it would hold quantitatively. 


Elliot and Dolbear® have made the very significant observation that 
during days of increased geomagnetic activity, the diurnal variations in both 
N and S directions are enhanced, and the N-S difference changes from being 
a semidiurnal to a diurnal curve. Since it is known that geomagnetic dis- 
turbances are associated with solar corpuscular streams, it would appear 
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that when there is increased activity of the solar M-regions, there is also 
increased cosmic ray emission from the sun. The radical change in the 
nature of the difference curve during magnetically disturbed days supports 
the view that it has no special physical significance apart from being the 
arithmetic difference of the diurnal variations in the two directions. 


Recently, two of us (V. Sarabhai and R. P. Kane) have demonstrated 
large and long-term world-wide changes in the amplitude and the hour of 
maximum of the meson diurnal variation. These changes, at least during 
years of low solar activity, are well correlated with relative sunspot number 
and the American magnetic character figure. There is therefore good reason 
to believe that continuous emission of cosmic rays from the sun is an im- 
portant cause for the diurnal variation of meson intensity. However, there 
are a number of points about the diurnal variation which are still to be 
cleared up. Perhaps the most important, concern the dependence of ampli- 
tude and hour of maximum of the diurnal variation on latitude, elevation 
and seasons. Our understanding of the problem is very much confused by 
the complex trajectories of charged particles from the sun. The very consi- 
derable deflection of cosmic rays from the sun is demonstrated by the 
occurrence both on the sunlit and the dark hemispheres of the abnormal 
increases of cosmic ray intensity associated with flares. It would be very 
valuable if the rather complicated problem of the effect of the geomagnetic 
field on non-isotropic cosmic ray primaries from the sun is tackled not only 
for a static case, but when the field wobbles with respect to the sun on 
account of the rotation of the earth. 


We are indebted to Prof. K. R. Ramanathan for many valuable dis- 
cussions. We wish to thank the Atomic Energy Commission of India for 
financial support given for our investigations and the India Meteorological 
Department for the supply of radiosonde instruments. We also owe our 
thanks to Mr. K. A. Gidwani for computational assistance. 


SUMMARY 


The daily variations of total cosmic ray intensity and the intensities 
of meson and electron components have been studied ‘at Ahmedabad with 
vertical geiger counter telescopes. The influence of meteorological factors 
on these variations has been examined, and it has been found that appro- 
priate barometric coefficients for correcting the cosmic ray intensities can 
be obtained from a consideration of the semidiurnal components of the varia- 
tions. The barometric coefficients for the three intensities are 


By = — 4:2% per cm. Hg, By = — 2°4% per cm. Hg, 
By = — 14-0% per cm. Hg. 
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The cosmic ray intensity variations are corrected with the appropriate co- 
efficients for the daily variation of barometric pressure. No signifi- 
cant variation is then left in the electron intensity, implying that variations 
of this component are mostly caused by the mass absorption effect with a 
variation of barometric pressure. In total intensity and in meson intensity, 
on the other hand, there is a significant residual variation of about -3% in 


amplitude. This is mainly diurnal in character with a maximum at 0900 
hours I.S.T. 


Reasons are given for concluding that the meson residual variation is 
not primarily caused by either the diurnal variation of temperature in the 
atmosphere or of geomagnetic elements. It is finally concluded that the 
bulk of the meson residual diurnal variation is extra-terrestrial in origin and 
is caused by continuous solar emission of cosmic ray particles. This con- 
clusion is discussed in terms of the interpretation of omnidirectional and 
unidirectional measurements of the diurnal variation by other workers. 
A connection between changes in the amplitude and the hour of maximum 
of the diurnal variation has been suggested. 
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INTRODUCTION 


A CRITICAL study of the phenomena near a heated surface is of fundamental 
interest to the physicist as well as the atmospheric physicist or the meteoro- 
logist. The convection pattern above a hot plate tells us much about the 
nature and origin of the instability of the air layers above the insolated 
ground in relation to the vertical gradients of air temperature, density, vapour 
pressure, etc. On the other hand, below the hot plate, conditions are obvi- 
ously stable and, with suitable experimental techniques, it is possible to 
investigate the phenomenon of thermal repulsion which is molecular in 
origin. 


The bibliography (Nos. | to 36) appended to this paper gives a list of 
papers on these topics published by the present writer and his collaborators 
during more than two decades after 1930. Some of these investigations 
deal with the meteorology of the air layers near the ground or micro- 
meteorology corresponding to what may be called phenomena “ above the 
hot plate”. We have also carried out numerous studies on the phenomena 
which may be grouped under the “ below the hot plate” type. 


Professor K. R. Ramanathan, while he was working at the Meteoro- 
logical Office, Poona, was much attracted by these investigations of ours. 
The present writer, therefore, takes this opportunity to dedicate this brief 
review to him on the occasion of his 60th Birthday. 


In what follows we shall present a very condensed review of the results 
achieved so far and then describe further investigations now in progress. 


2. PHENOMENA VERY CLOSE To A HoT SURFACE 


If a small heated! ? plate, say about 2 inches square, is suspended hori- 
zontally at the centre of a box with glass windows the air inside which is 
rendered smoky and the space above and below the surface examined under 
suitable illumination, one observes that immediately near the hot surface 
there is a thin film of dust-free air which looks dark. The thickness of the 
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dark layer is uniform and comparatively small, of the order of 1/10 mm. or 
so, below the plate. Above the hot plate the thickness of the dark layer 
increases near the centre where it ascends as a pillar or tongue or rising 
column with rapidly decreasing cross-section. The general appearance 
of the dark layer and the movement of the smoke particles is as shown in 
Fig. 1 (a). 

This pattern remains fairly steady both above and below the hot plate. 
If the hot plate is made quite extensive, say one foot or more square, the 
appearance below the plate remains steady as before, but conditions above 
the plate become very much more lively and fascinating. The dark layer, 
which really represents the hottest air close to the plate, now ascends, not in 
one steady dark pillar as in Fig. 1 (a), but at a large number of places as ascen- 
ding tongues of hot air rushing up into the colder air above, extending several 
centimetres upwards. While the base of the ascending tongues near the 
hot plate is 14 to 2} cm. broad, their width diminishes rapidly upwards. 
These ascending columns develop and move about in a random manner, 
sloping with the wind, if any, across the plate. The appearance at any 


tat 


Hot Plate (2cm. square) 
Fic. 1 (a) 
instant is shown. in vertical section in Fig. 1 (b) where H, H,.... represent 


the rising hot air filaments and C, C, C, etc., the compensating downward 
movements of the colder air from above. 


Fig. 1 (b) shows how the dark layer is in unstable equilibrium with the 
colder layers above it so that it breaks through the cold air, producing also 
the compensating downward currents of cold air. 
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Even when there is wind (horizontal), this pattern{is not wiped off, but is 
only inclined in the direction of the general air movement. 


It will be obvious that the picture shown in Fig. 1 (5) really represents 
what happens above the heated ground during the day-time. This special 
type of convection is responsible for the well-known “ shimmering” of 
distant objects when they are viewed through the air layers near the ground 


during a clear day with strong insolation. Indeed, we may call this layer 
the shimmering layer ”’. 


We may now continue the examination of the “ shimmering layer”, 


reserving the discussion of the phenomena “ below the hot-plate ” for a later 
section of this paper. 


3. THE SHIMMERING LAYER AND SOME OF ITs SPECIAL FEATURES 


(a) Enormous Temperature Lapse Rates—From the hot surface the 
temperature falls very rapidly with height at first and less and less rapidly 
as one moves away further and further.* Using an interferometric technique 
it has been shown"! that within the layer of a millimetre or two nearest 
to the surface, the variation of temperature is most rapid, but sensibly linear. 
This is the layer too which develops inferior mirages.* * 


As an example of the variation of temperature and lapse-rate with height 
above bare ground we have Table I which gives the mean daily values recorded 
during April 1950 at the maximum temperature epoch at the Central Agri- 
cultural Meteorological Observatory at Poona. 
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TABLE I 


Variation of Temperature °C. and Lapse-Rate in Degrees per Kilometre at 
the Maximum Temperature Epoch at the Central Agricultural 
Meteorological Observatory at Poona, April 1950 


Height above- 
ground in — Lapse-rate in °C, 


centimetres : per kilometre 


2-02 x 10° 
7+33 x 104 
2-60 x 104 
1-33 x 104 
5+33x 102 
2-00 x 102 
1-67 x 
1-33 x 108 
5-00 x 10? 
5-00 x 10? 
1-67 x 10? 
1-67 x 102 
1-11 x10? 
2-00 x 10? 
| 1-33 x 102 
| 


NW 


6-06 x 102 


Note that the lapse-rate can be about 2,00,000 times the adiabatic lapse- 
rate just above the ground and that it decreases very rapidly with height 
attaining a value of the order of only 6 times the adiabatic lapse-rate about 
8 metres above ground. 


(b) Fluctuation of Temperature in the Shimmering Layer.—A sufficiently 
sensitive temperature recording element will record a fluctuating temperature 
if placed in the shimmering layer. Using extremely small thermocouples 
of copper and constantan (45S.W.G.), one of them kept in a bath ata 
standard temperature and the other inserted in the position of the thermo- 
meter in the smaller-sized Assmann Psychometer, connected to a sensitive 
Mill-Galvanometer with a period of 4 second, it has been possible to record 
practically instantaneous values of air temperature at intervals of a few 
seconds. After aspirating the psychrometer, a series of 50 readings was 
taken at 10 second intervals at various levels above ground up to 35 feet, 
at the Central Agricultural Meteorological Observatory at Poona. After 
smoothing, so as to eliminate long-period or slow changes, values of (i) the 
mean temperature, (ii) the standard deviations of short-period fluctuations 
and (iii) the highest and lowest temperatures recorded at each of the levels 
were calculated. The values obtained at the maximum temperature epoch 
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on the 6th January 1942, a typical clear day, are plotted in Fig. 1 (c) the 
height scale being logarithmic. XX and NN are the highest and lowest 
temperatures, while MM shows the mean temperature. AA and BB show 
x 


VANE 
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AM 8 


i i 
32° 34° 
Temperature 


Fic. 1 (c). Fluctuation of air temperature at different levels above ground during the 
interval 1430 to 1530 hrs. L.S.T. on 6th January 1942, Heights above ground in inches’ 


30° 


d feet’. 


the mean temperature “‘ minus” and “‘ plus” the standard deviation res- 
pectively. The rapid variation with height in the amplitude of the tempe- 
rature fluctuations within the “‘ shimmering layer” is shown clearly in 
Fig. 1 (c). On plotting the standard deviation against the logarithm of the 
height it is found that the points all practically lie on a straight line which 
cuts the height axis at a point corresponding to about 170 feet above the 
ground. This shows that at the maximum temperature epoch the shimmering 
layer extends up to 150 or 200 feet above ground. The above method pro- 


vides an easy method of estimating this height from observations near the 
ground.'® 1° 


Krishna Rao, at present working in our laboratories, has succeeded in 
photographically recording the fluctuations with the aid of a quick-run Moll 
Recorder. The temperature is recorded with a 40 S.W.G. Copper-Cons- 
tantan thermocouple coated with magnesium oxide and exposed at various 
levels above ground (Central Agricultural Meteorological Observatory) and 
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connected to the same sensitive Moll Galvanometer referred to above. 
Using this technique; Rao is studying the characteristics of these fluctua- 
tions and their variations with height above ground as well as with time of 
the day. These results will be discussed in detail in a later paper. 


Fig. 2 shows some typical two-minute records at various heights above 
ground and at 0630 hrs. (minimum temperature epoch), 1000 hrs., 1400 hrs. 
(maximum temperature epoch) and 1700hrs. on the 16th of November 
1952. Near each record the mean temperature in degrees centigrade is 
recorded. The fluctuations of temperature are very slight or feeble at all 
levels at the minimum temperature epoch (0630 hrs.), but there is a slight 
increase with height up to 4 or 5mm. above ground and a rapid deérease 
thereafter; the fluctuations being negligible above 5 cm. level. 


Some time after sunrise the temperature is rising rapidly and the records 
obtained at 1,000 hrs. clearly show that the shimmering pattern is develop- 
ing rapidly and has already grown to a level beyond 9 metres. The next 
series of records relate to the maximum temperature epoch (1400 hrs.). 
Here the amplitude of the fluctuations is relatively greater than before at all 
levels. The variation with height both at 1000hrs. and at 1400hrs., 
however, bring out some very interesting features. The amplitude of the 
fluctuations is least at the ground surface and increases rapidly with height 
up to a few mm. (say 3 mm. at 1000 hrs. and 1 cm. at 1400 hrs.), attaining 
a maximum value. The amplitude of the fluctuations decrease thereafter 
with height as indicated in Fig. 1 (c). 

A reference to Fig. 1 (5) will show that up to a few millimetres above 
ground the temperature measuring element will be more often than not 
within the dark layer and will begin to be influenced by both the uprising 
hot and the down-coming cold currents of air only at some distance above 
the ground. The amplitude of the fluctuations will, of course, depend also 
on the contrasts of temperature in these air currents. 


From the last set of records at 1700hrs., the rapidity. with which the 
upper limit of the shimmering layer sinks downwards towards the ground, 
and the decrease in the amplitude of the thermal fluctuations will be quite 
clear. 


A fuller discussion of these results will be presented elsewhere. 


Y. D. Altekar*® working in our laboratory is engaged in building an 
electronic variance-meter designed by him for recording the sums of squares 
of fluctuations of air temperature, humidity or wind velocity. When 
completed this equipment will provide another useful tool for these 
investigations, 
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(c) The Shimmering Laver in Relation to the Growth and Decay of Thermal 
Inversion during Clear Nights.—The upper limit of the shimmering layer 
extends upwards beyond 200 feet (about 60 metres) at the maximum tempe- 
rature epoch and thereafter comes down rapidly towards sunset. In tropical 
and sub-tropical and lower temperate regions. during clear nights with light 
or no air movement, the shimmering layer about a foot or less in thickness 
is found to persist throughout the evening and night. The coldest air is 
therefore at some distance above ground at the top of the shimmering layer 


which, of course, is a iayer of “lapse” (or temperature decreasing with 
height).*-* 


The investigations at Poona show that in an open space with bare ground 
the nocturnal inversion layer begins forming some time before actual sunset. 
Its lower boundary merges into the upper boundary of the shimmering layer 
which also begins to contract. The upper boundary of the shimmering 
layer indeed approaches the ground rapidly about sunset and within a short 
time after sunset attains its minimum height, about a foot or so above the 
ground; thereafter, it remains at this level, more or less, during the rest 
of the night and until insolation sets in next morning when it grows rapidly 
again. The transition between the developing inversion and the upper 
boundary of the shimmering layer is a region of zero temperature gradient 
and may be called the “‘ lower isothermal layer”’. This layer moves down- 
wards in the evening about sunset and upwards in the morning after sunrise. 


It may be observed that while the base of the “ inversion” layer is 
associated with the top of the ‘* shimmering layer ” in the manner described 
above, the inversion layer itself builds up during the night, growing in thick- 
ness. Its upper boundary merges into the “ upper isothermal layer ” above 
which again there exists a layer with a lapse or decrease of temperature with 
height. Thus the inversion layer near the ground is cushioned, as it were. 
between two layers of lapse, one of which, the lower one, is the remnant of 
the shimmering layer and the other is the free atmosphere above in which 


also temperature decreases with height. The rise of temperature with height 
within the inversion layer is quite pronounced. 


The growth of the inversion layer during the night and its destruction 
after sunrise have been investigated at Poona in some detail. While the 
psychrometric readings of air temperature up to 35 feet (with the aid of the 
observatory tower) taken at short intervals of time provide adequate informa- 
tion about the base of the inversion, the behaviour of its upper boundary 
has been investigated recently with tethered radiosonde instruments of the 
F. Type constructed by the Instruments Division at the Meteorological Office, 
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Poona.’ *8. 99 The actual instrument used in our experiments had a more 
open temperature scale than the standard instruments used for upper air 
investigations. From these observations recorded at the Central Agri- 
cultural Meteorological Observatory on the 25-26th November and 
98-29th December 1949, it is observed that by 2200hrs. the thickness of 
the inversion layer is of the order of 100 metres, increasing to 200 metres by 
0300 hrs. and to 250 metres just before sunrise, on the next day. These 
results will be discussed elsewhere, but it is clear that at Poona the layer can 
grow up to 250 or 300 metres above ground. 


Figs. 3 (a) and (b) are idealised diagrams constructed to show the deve- 
lopment and destruction respectively of the inversion layer. In both these 
figures the ordinates represent the height above ground (logarithmic scale) 
and the abcisse temperature. The temperature height curves are intended 
to show the nature of the variations which set in before, during and after 
sunset [Fig. 3 (a)] and just before, during and after sunrise [Fig. 3 (d)]. 


Height above ground in metres 


Temperature ——> Sun-set Time 


Fic. 3(a). Development of Nocturnal Inversion. The above is an idealised diagram in 
which the curves 1 to 10 would correspond roughly to the following hours for a clear winter 
day at Poona. 

Approximate time in hr. 
No. New LS.T. 


10 0700 


N.B.—The temperature range from B to C is of the order of 10°C. The height scale is 
logarithmic, ground surface to level of BC is of the order of 4 metre and O is about 
6 to 7 metres above grouud. 
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Fic. 3(6). Destruction of Inversion after Sunrise. The above is an idealised diagram in 


which the curves | to 4 would correspond roughly to the following hours for a clear winter 
day at Poona. 


Approximate time in hr. 
Curve No. New LS.T. 


1 0700 
2 0800 
3 0900 
4 1000 


N.B.—The temperature scale is as in Fig. (3)a. BC is about 4 metre above ground and 
AO' about 300 metres above ground. 


Curves 1, 2, 3....10 of Fig. 3 (a) represent the vertical thermal struc- 
ture at different times. As the temperature decreases with time at all levels, 
the time sequence is from right to left. Curves 1 and 2 show the thermal 
structure before the formation of the inversion layer. Sometime before 
sunset, when the inversion is about to commence, the upper limit of the 
shimmering layer is at the point O. As soon as nocturnal cooling starts, 
the inversion starts developing. Its lower boundary, ‘.e., the ‘* lower iso- 
thermal layer” approaches the ground rapidly from O to B within a short 
while after sunset and thereafter remains more or less steady at about a foot 
or so above the ground as BC. The upper boundary of the inversion layer, 
starting at O, begins to rise rapidly at first and more gradually thereafter 
as OA in the figure. AOBC is the inversion region in the diagram. OA 
and OBC are in the upper and lower isothermal layers respectively. Below 
OBC we have the lapse of the shimmering layer! and above OA the lapse 
of the free atmosphere. 


Curves 1, 2, 3 and 4 of Fig. 3 (6) represent the thermal structure at 
different times before and after sunrise. As the temperature increases with 
time, the time sequence is now from left to right. The inversion region is 
AO’BC. AO?’ is the upper boundary of the inversion layer in the morning 
which remains more or less steadily at the same level during the process of 
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destruction of the layer. CB represents the base of the inversion layer 
early in the morning.’ As soon as insolation begins at the point indicated 
by the vertical arrow (sunrise), the shimmering layer thickens very rapidly 
from B to C. As soon as it reaches O’ the inversion is completely 
obliterated. 


(d) Other Phenomena.—We have also carried out numerous investiga- 
tions on the invisible condensation of water vapour during the night on the 
soil,®: 22-8 on the microclimates of plant communities,” * 32 on the fate 
of solar radiation,**-%* on the correlation between the vertical gradients of 
wind velocity and air temperature,”° etc. 


4 THERMAL REPULSION 


Below the heated plate, owing to the stable arrangement of the air layers, 
the convective phenomena are not violent as in the previous case but very 
much simpler. In fact, when matters are so arranged that an air cell (filled 
with smoke and illuminated suitably) is formed with a hot surface above 
and a cold surface below and the sides are suitably enclosed, the dark layer 
is found to be restricted in its scope by a pair of vortices. As the cold lower 
surface is brought nearer and nearer to the hot surface it is found that the 
vortices separate towards the two sides leaving a calm layer at the centre 
where the dark layer attains a maximum thickness; ultimately, when the 
cold surface approaches the hot one within 2 mm. or less the vortices die 
away completely, all convective phenomena cease, and thermal repulsion 
has full play. When this “ convection-free ” state is reached, dust particles 
are repelled with a uniform velocity as defined by Stokes’ Law and are 
deposited on the cold plate. The velocity is found to be proportional to 
the thermal gradient. These results have been discussed fully by Ramdas,}2 
and Paranjape.!* Later, Ramdas and Joglekar® studied the movements of oil 
droplets in a vertical convection-free cell (between a vertical hot surface 
and a vertical cold surface) and found that the falling particles being acted 
upon by gravity vertically and by the thermal repulsive force horizontally, 
move in straight lines inclined to the vertical from the hot to the cold surface. 


In the same paper they have described experiments on the steady deflec- 
tion of a mica vane suspended parallel to the two surfaces by means of 
quartz fibre. The mica vane is repelled from the hot towards the cold sur- 
face, the deflections being proportional to the temperature gradient. The 
apparatus used in the above experiment is shown in Fig. 4. Here the mica 
vane M is suspended by a quartz fibre between the two surfaces G H (hot) 
and KL (cold). These surfaces are maintained at the desired temperatures 
by circulating hot and cold water respectively through the pipes C,-C, and 
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Fic. 4. The Thermal Repulsion Apparatus 


C,-C,. The distance between GH and KL is adjustable. The joints 
C, D, E and F are made air-tight and the interspace can be evacuated to any 
desired pressure. When there is a difference of temperature between GH 
and KL the mica vane is deflected to the right by an amount proportional 
to the temperature gradient. 


Fig. 5 represents a typical curve showing the dependence of thermal 
repulsion on the gas pressure. Here the distance between the hot and cold 
surfaces was 5mm. 64 is the deflection in divisions of the microscope eye- 
piece scale and AT the difference of temperature. One division on the 
microscope scale equals 1/20mm. and a unit value of 5/AT corresponds 
to a thermal pressure of 5-659 x 10-3 dyne per cm.? 


It will be seen that the deflection is very small at high pressures but 
increases to a maximum value as the air pressure is reduced to 6-3 x 10-3 cm. 
of Hg. The deflection decreases thereafter with further decrease and aif 
pressure. Obviously the left-hand portion of curve where the deflection is 
proportional to the gas pressure, enables us to design a simple gauge for 
recording low pressures. Paranjape’® !? working under the present writer 
has discussed in great detail the results he obtained for various gases, kinds 
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of vanes suspended, and with the vane parallel and perpendicular respec- 
tively to the surfaces GH and KL. He has also shown how the accommo- 
dation coefficient can be deduced from the last mentioned experiment. It 
may also be mentioned in conclusion, that in our experiments, thermal 
repulsion is shown in all its simplicity as a molecular phenomenon without 
the complications introduced in the classical Crooke’s Radiometer experi- 
ments. 


The adaptation of the apparatus shown in Fig. 4 for use in radiation 
measurements is also under investigation. 


5. CONCLUSION 


In the present paper a brief survey of some of the important experi- 
mental results obtained above and below a hot plate has been made. The 
subject is still full of new problems. 
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The present writer wishes to take this opportunity to thank all his colla- 


borators who have helped to make progress in these studies. 
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EAST-WEST ASYMMETRY OF COSMIC RAYS 


By D. S. R. Murty 
(From the Physical Laboratories, Osmania University, Hyderabad) 


Received August 19, 1952 
1. INTRODUCTION 


Tur theory of motion of electrically charged particles in the earth’s magnetic 
field has been developed by Stormer! and by Lemaitre and Vallarta.” 
According to their theory, at a given latitude and in a given direction only 
particles having energy above a critical energy are able to break through 
the earth’s magnetic field. The cones into which the allowed particles sepa- 
rate from the forbidden ones, start opening from the western sky for the 
positively charged particles and from the eastern sky for negatively charged 
particles. If there is preponderance of particles of one charge over those 
of the opposite charge there must exist the east-west asymmetry. Rossi® first 
tried to verify this experimentally with a two-fold coincidence telescope but 
failed to establish it. The first experimental evidence of this asymmetry 
came with the work of Johnson and Street.’ It was later confirmed by many 
others. From the most extensive survey of the east-west asymmetry by 
Johnson® it is found that the asymmetry is maximum (13%) near the geo- 
magnetic equator and decreases towards the higher latitudes being 2 to 3% 
at A= 50° at sea level. For a given latitude, asymmetry increases with 
altitude. 


The present work is undertaken at Hyderabad, India; A= 7-7° the 
height above mean sea level being 1,800 feet. 


2. EXPERIMENTAL ARRANGEMENT 


The cosmic ray telescope employed for the present study of east-west 
asymmetry consists of three Geiger Muller counters mounted in parallel posi- 
tions on a light wooden framework which is capable of rotation about 
horizontal as well as vertical axis, the angles being read accurate to } of a 
degree by a pointer moving on a graduated circular scale. The counters 
are self-quenched with copper cathodes and tungsten anodes and filled with 
a mixture of argon and alcohol to a pressure of 11 cm. of mercury in a glass 
envelope. The threshold voltage of the counters starts at about 1,050 volts 
with a plateau of about 200 volts. The cathodes are 3-5cm. in diameter 
and 40cm. in length. With a distance of 25-5cm. between the two 
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extreme counters the triple train subtends an angle of 15-6° in the plane in 
which the telescope is rotated. 


To keep the high tension applied to the counters independent of fluctua- 
tions of the A.C. mains, a circuit due to Neher and Pickering® is used. The 
voltage could be controlled and kept constant at the required value, the 
regulation with the circuit being highly satisfactory. 


In the circuit used for recording the triple coincidences the counters 
are connected in the usual Rossi parallel connection preceded by one stage 
of amplification. A typical hard tube circuit is used for recording (Janossy’). 
Pulses are fed to the grids of the amplifying valves through shielded wire. 
To test the accuracy of the mechanical recorder, that is, to decide whether 
it missed any counts or not, the coincidence pulse from the plate circuit of 
the Rossi stage is applied to the vertical plates of a cathode-ray oscillograph. 
The triple coincidences were visually counted on the oscillograph screen and 


also recorded by the mechanical recorder and the two counts were found 
to be equal. 


The assembly was mounted on the roof of the Physics Laboratory, 
Osmania University. The telescope was practically open to the sky. The 
axle of the framework on which the counters were mounted was oriented 
magnetic north and south so that the rotation would be in the local east-west 
plane. The following method was adopted in taking observations. The 
telescope was oriented at a particular zenith angle to the east for a certain 
time and the counts recorded. The telescope was then rotated and 
orientated at the same zenith angle to the west for the same time (30 minutes 
to one hour) and the counts recorded. This process was repeated a number 
of times and the total number of counts from the east and west obtained. 
Readings were taken after every interval of 15°. This method eliminates 
any error due to the changes in instrumental selectivity, short period changes 
of cosmic ray intensity or any other changes due to variation in baro. 
metric pressure or any magnetic disturbances. 


Johnson® has shown that the only systematic error which merits discus- 
sions is due to accidental counts. When the horizontal rate was less than 10% 
of the vertical rate no correction for accidentals was needed. In the present 
case it is about 8% and, therefore, no correction was applied for accidentals. 


The probable error (P.E.) was calculated from the total number of counts, 
where P.E. = 0-67 1/1/C, where C is the total number of counts in time T. 
The asymmetry was calculated from the formula a = 2 (1,, — 1,)/(1, + 1), 
where I,,, and I, are respectively the counting rates from west and east. The 
results are given in Table I and shown graphically in Fig. I. 


East-West Asymmetry of Cosmic Rays 


% Asymmetry 


15 30 45 60 75 
Zenith Angle in Degrees 
Fic. 1. E.W. Asymmetry Curve (Lines Probable Errors) 
TABLE I. East-West Asymmetry 
(Hyderabad, India ; \ = 7-7° N. and altitude 1,800 Ft.) 


Time in | Total | Counts per minute % Asymmetry 
Zenith angle minutes | counts _ 2 (Iw — Ig) 
6 


T Cc pt PE. x 100 


— We) 


7-160 
7-599 5-920 


5-610 - 
6-148 9-139 


3-790 
4-279 


2-295 
2-499 


1-165 
1-211 
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3. RESULTS AND DISCUSSION 


From the graph in Fig. 1 it is evident that the asymmetry first rises to 
a maximum of 12-3% at a zenith angle of 43° and then falls off continuously 
with the zenith angle. The results are in general agreement with those of 
Johnson® !° and Korff,!! though some of the workers, Bhattacharya,!* Sarma 


and Sarna!* in India reported slightly higher values of asymmetry at higher 
latitudes. 


Purely from theoretical considerations, asymmetry should increase 
continuously with the zenith angle. That will be true only for primary 
particles. The field sensitive particles are absorbed by the thickness of air 
that they have to traverse before they reach the surface on the earth. The 
same explanation is offered for higher values of asymmetry at high altitudes 
where more of field sensitive particles can reach due to lesser absorption. 


The author expresses his grateful thanks to Prof. S. Bhagavantam for 
his helpful guidance and constant encouragement in the progress of the 


present work. He also thanks Mr. M. Krishnamurthi for his help in the 
course of the work. 


4. SUMMARY 


A triple coincidence counter telescope has been employed for the study 
of the east-west asymmetry of cosmic rays at Hyderabad, India (A = 7-7°; 
h = 1,800 ft.). The telescope, subtending an angle of 15-6° in the plane 
of rotation, is mounted on a light wooden frame capable of rotation about 
horizontal and vertical axis. Asymmetry measurements have been made 
at a number of zenith angles with an interval of 15°. The asymmetry rises 
to a maximum of 12-3% at a zenith angle of 43° and falls off at higher angles. 
The results are in general agreement with some of the other workers. 
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HEIGHT DISTRIBUTION OF ATMOSPHERIC OZONE 


By K. R. RAMANATHAN AND R. N. KULKARNI 
(Physical Research Laboratory, Ahmedabad) 


Received January 27, 1953 


In extra-tropical latitudes, fluctuations in ozone amount occur from day 
today! These are most marked in late winter and early spring and least 
marked in late summer. In low latitudes, the inter-diurnal variability is 
generally small; in N. India, during the months December to April, when 
active western disturbances pass across the country, fluctuations of a similar 
character occur, but are less frequent and of smaller amplitude. This can 
be seen from Fig. 1 which shows the daily values of ozone over Mount Abu 
(Lat. 24° 26’N., Long 72° 43’ E., 3,900 ft. above sea level) in the period 
agaenal 1951 to June 1952. 


DAILY OZONE AMOUNTS 
Mount Asu — Nov. 1958 To JuNE 1952 


VM 


Vio Jan. 20 ‘ 10 Mar. 20 


i0 Apr. 20 0 May 20 9 June i9 3 1952 
Fic. 1. a ozone amounts at Mount Abu 1951 to June 1952 
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The inter-diurnal variation of ozone at a station in Europe is known 
to be closely related to meteorological conditions, the closest correlation 
perhaps being between the ozone amount and the height of the tropopause, 
the ozone amount increasing with a lowering of the tropopause.” * In India, 
with the limited high-level upper air data available, it appears that even 
during winter, the correlation between the height of the tropopause and the 
pressure field at ground level is much less close than in temperate latitudes; 
so also is the association between ozone and surface pressure.’ Still, it was 
felt that some insight would be obtained into the connection between ozone 
and weather if in the period of western disturbances, not only the changes 
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disappeared off Kathiawar on 17th. 
Feb. 5—Following W. disturbance which was active in Punjab, U.P. 


Feb. 7—Disturbance from south affecting Madhya Pradesh. 
Mar. 14—W. disturbance centred between Bikaner and Gwalior, 


Mar. 22—W. disturbance located N.-E. of Abu. 


Height Distribution of Ozone at Mount Abu 
Noy. 13-17—Arabian Sea storm appeared near 12° N., 65° E. on 13th, moved north and 


Apr. 2—W. disturbance over N. Rajasthan, U.P, 


Weather Remarks :— 
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of ozone amount but also the changes in vertical distribution could be 
studied. This is a difficult problem, because, close to a centre of low pressure, 
clouding prevents the collection of a connected series of clear zenith sky 
observations which are necessary for determining the height distribution of 
ozone. Mount Abu is however well situated for attempting such observa- 
tions, and it has been possible to collect umkehr observations of the zenith 
sky on 112 days during the period October 1951 to November 1952. Some 
of the observations were on consecutive days. The present paper contains 
a summary of the results so far analysed, followed by a discussion of the 
variations in the height distribution of ozone corresponding to different 
ozone amounts, first over Mount Abu, and then over other places on earth 
for which data are available. The consideration of the Indian data in rela- 
tion to weather will be taken up on a later occasion. 


2. The zenith sky intensity observations were taken on two pairs of 
wavelengths, the usual pair A 3112/A 3323 and another pair A 3075/A 3278.* 
The method of calculating the height distribution is the same as that explained 
in a previous paper.’ It was found that practically identical results were 
obtained from observations on either pair of wavelengths made on the same 
day (Sample umkehr curves on the two pairs of wavelengths are given in 
Fig. 2 (5-2-1952 and 7-2-1952). 


3. Table I gives the height distribution of ozone on 34 days in the 
period October 1951 to November 1952 for which the data have been calcu- 
lated. Of these, 29 days fall in the months November to April. The ozone 


amounts and the evening pressures at Mount Abu at 1730 hours I.S.T. have 
also been given. 


The following points are noteworthy. 


(1) Most of the changes in ozone take place in the layer 18 to 27 km., the 
amount in that layer increasing with increasing total amount. The largest 
increases have occurred shortly before, or shortly after, the passage of centres 
of western depressions across the longitude of Abu. 


(2) There is no detectable ozone in the layer 0-9 km., and only a very 


small quantity in 9-18km. The latter shows a slight increase when the 
ozone amount is high. 


* The usual method of computing the height distribution of ozone from umkehr obser- 
vations with the assumption that the vertically scattered skylight is entirely made up of 
primarily scattered light is open to objection and probably leads to too low a value of ozone 
in the first, second and third layers. This question is under investigation (Q.J. Roy. Met. 
Soc., 1952, 78, 625). It is however believed that the general nature of the discussion followed 
in this paper will not be invalidated.—K. R. R. 
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Fic. 3. Ozone amount in 18-27 km. and above 27 km. for different total amounts at 
Mount Abu. 
the (3) When the ozone amount is abnormally high as on 5-2-1952 and 
Best 7-2-1952, there is a tendency for the amount above 27 km. to be slightly 
tres smaller than usual. 
Fig. 3 is a dot diagram showing the ozone amounts in 18-27 km. and 
/ery above 27 km. plotted against the total ozone amount on individual days. The 
the mean data relating to Delhi and Kodaikanal are also plotted in the same 
figure. The high correlation of the amount in 18-27km. with the total 
aii amount is obvious. This suggests that air-movement and convergence 
> a associated with low pressure centres of western depressions sometimes 
ozone extend well into the stratosphere. 
Met. 
lowed In Table II, the data have been grouped according to the total ozone 


amount, and the corresponding mean distributions are given, 
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TABLE II 


Mean Height Distributions for Different Ozone Amounts (Grouped Means) 
at Abu, Delhi, Poona and Kodaikanal 


(Height distribution 10-8 cm.) 


| 
Station | 3 
| 10-* cm. 0-9 9-18 1827 27-36 36-54km. 


ooooceo 


Poona 


Kodaikanal 


It may be remarked that the ozone distributions in September, October 
and November do not show any appreciable difference from those in the 
months December to April, when they relate to the same ozone amount. 


4. Comparing the ozone distributions over Mount Abu with those 
over Delhi and Kodaikanal (Table III), the main features are similar. A 
small amount of ozone seems to get down to the lowest layer 0-9 km. at 
Delhi during winter. While the amount above 27 km. is practically the 
same as at Delhi and Abu, it is nearly independent of the total amount. Over 
Kodaikanal in South India for which umkehr curves are available for only 
5 days in March,® a larger amount of ozone is situated between 27 and 
36km. The abnormally high figure for the ozone amount above 36 km. 
for Poona is not understood. 


5. The highest ozone amount for which we have height distribution 
data in India is 0-217 cm. It is interesting to extend the range by consi- 
dering the data of Arosa and Troms¢ and see how they compare. For this 
purpose, it is convenient to estimate the amounts of ozone over these places 
for the same height intervals 0-9 km., 9-18 km., etc. This has been done 
from the data given in the following two papers. 


(i) Ozone height distribution data for Arosa. By Gétz, Meethan 
and Dobson’ (Method B). 


(ii) The same for Troms¢. By T¢nsberg and Olsen® (Method A). 


|| 
| 
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Abu 153 (4) 9 60 68 15 
159 (11) 9 68 68 15 
3 172 (8) re 76 70 18 
180 (6) 10 86 60 19 
: 195 (4) 11 106 57 18 1 
209 (1) 4 il 64 12 
; Delhi - 155 (1) 15 49 69 13 
175 (3) 18 63 72 13 l 
200 (1) 23 91 73 13 
| 217 (1) 27 94 62 15 
: — “ 169 (8) 1 13 65 65 25 ' 
TY ve 175 (4) 0 6 53 99 17 
| 
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Knowing the uncertainties associated with the determination of the 
height distribution of ozone from the zenith sky measurements, not much 
accuracy can be attached to the figures but they are good enough for a 
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from Abu (24°-4 N.) to Troms (69°-7 N.) 


* In Fig. 4 (a), the results obtained from a V, rocket ascent have been added 
(F. S. Johnson, ef al., Jour. Geophys. Res., 1952, 57, 154). 
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general survey and at any rate, nothing better is available. The data have 
been piotted in two composite diagrams [Figs. 4 (a)* and 4 (6)] with total 
ozone amounts as abscisse and the amounts above the levels 36 km., 27 km., 
18 km. and 9 km. as ordinates. These diagrams enable us to make a few 
general statements :— 


(1) The ozone amount above 30 km. or thereabouts remains nearly the 
same irrespective of the total amount of ozone.°® 


(2) In low latitudes, more than half the ozone is concentrated between 
18 and 27km. This amount increases with increasing ozone amount, and 
the increase continues for the still larger amounts that occur over Arosa 
and Troms¢, tending to become constant for the highest value of ozone. 


There are substantially larger amounts of ozone over N. India in winter 
between 18 and 27 km. than over Arosa between the same levels for the same 
ozone amount. This may be partly due to the fact that low ozone values 
over Arosa refer to summer conditions and the data are therefore not strictly 
comparable. Fig. 4(b) also suggests that as we go from north Indian lati- 
tudes to the latitude of Arosa, part of the ozone in the layer 18-27 km. 
descends to lower levels and is distributed in a greater thickness of the . 
atmosphere. 


(3) Below 18 km., the ozone is present only in small quantities in low 
latitudes and probably has a maximum in the cold season. It increases with 
latitude from Kodaikanal (10° N.) to Delhi (284° N.). 


(4) In middle and high latitudes, there are substantial amounts of ozone 
below 18 km., and the larger the total amount, the larger is the proportion 
in the lowest levels below 9 km. For the largest ozone amounts, that occur in 
high latitudes (0-350 cm.), nearly half the quantity may be below 20 km. 


DISCUSSION 


6. One of the most remarkable features of the distribution of ozone 
which came out from the world survey organised by Dobson was the general 
increase of ozone from the Equator towards high latitudes and its strong 
annual variation in middle and high latitudes with maximum in spring and 
minimum in autumn. Ozone is a product of the photo-chemical action of 
ultra-violet solar radiation on the oxygen of the upper atmosphere, and if 
the atmosphere were stagnant, the rate of formation of ozone by light of 
wavelengths shorter than 2420 A° would balance the rate of its destruction 
by light of wavelengths < 11400 A° and an equilibrium amount with an 
appropriate vertical distribution would be established depending on the 
latitude and the season. The height of maximum ozone amount will rise 
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with increasing zenith distances of the sun. Such a photo-chemical theory 
first suggested by Chapman and worked out in greater detail by Wulf and 
Deming and others, explains generally the observed vertical distribution in 
the tropics, but does not explain either the actual ozone amounts or the 
dependence on latitude and season. 


K. Langlo in his recent memoir® ‘On the Amount of Atmospheric Ozone 
and its Relation to Meteorological Conditions” has given a comprehensive 
survey of the variations of ozone in the atmosphere. We can supplement it 
by the following remarks: 


(1) Fig. 4 (6) shows that approximately 0-100 cm. of ozone is contained 
in the atmosphere above 27 km. in all places where there is sunshine. Within 
the Arctic Circle during the polar night, this amount seems to disappear. 
Langlo has compared the five-day mean ozone values at Dombas (62° N.) 
and Troms¢ (70° N.) in the period December 1945 to January 1946. 
Troms¢é values are lower by about this amount, providing evidence for the 
view that the ozone above 27km. or thereabouts is in photo-chemical 
equilibrium with sunshine and if the sunshine is removed.for a few days, 
the ozone above this height disappears. It is interesting to note that if this 
view is accepted, the annual ozone variations at Troms¢ will fall into line 
with those at Oslo and Arosa and not show any anomaly. The ozone 
observations now being taken at Spitzbergen under International auspices 
will show whether this view is correct or not. 


(2) Most of the ozone which is situated below 27 km. is accumulated 
ozone and is in the stratosphere. Where the tropopause is permanently high 
as in the tropics (16-17 km.), the amount is permanently low. The high 
negative correlation between the height of the tropopause and the ozone 
amount in temperate and polar latitudes, which seems to apply even to the 
region of tropopause funnels elucidated by Palmen,!® shows that in general 
the lower the tropopause and the greater the accommodation available in the 
stratosphere, the larger is the ozone amount. The stratosphere below 27 km. 
seems to be a safe place for ozone. The extreme dryness of the stratosphere 
may have something to do with this.’ 


(3) The almost complete absence of ozone in the tropical troposphere 
is a fact of great significance. It would appear that water vapour and other 
material, oxidisable or catalytic, which are carried up into the upper atmo- 
sphere from the lower layers by convective currents are responsible for. 
destroying the ozone in the upper layers of the tropical troposphere. The 
destructive action may be aided by light filtered through ozone. This view 
is supported by the rapid disappearance of the occasional incursions of 
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ozone to the sub-tropical troposphere and the steady decline in ozone 


amount in late spring and summer in the atmospheres of middle and high 
latitudes. 


(4) Over “ cold sinks ”’, on earth such as the polar regions and large 
continental areas in winter where the tropopause is already low, some of 
the stratospheric air can sink into the troposphere and ride over inversions, 
retaining its ozone content for fairly long periods. Ehmert’s observations 
made on aircraft bring this out clearly.!2 The ozone will quickly disappear 
when it gets into regions of convective activity. 


(5) Langlo has noted that on the mean of the year, the ozone amount is 
nearly constant from Troms¢ to Arosa with a mean value of about 0-235 cm. 
and that it drops steeply to 0-190 cm. at 25° N. This rapid change runs 
parallel to the change in the average height of the tropopause. It is only in 
the season of western disturbances in N. India when composite tropopauses 
are frequent that marked day-to-day fluctuations in ozone occur. 


(6) The average tropopause has a height varying from 17 km. to 9 km. 
with a rather steep drop or break between 16 and 12km. and yet other 
breaks near polar fronts. Although the air between 9 and 16 km. is highly 
stratified, it is subject to large-scale movements, both vertical and horizontal, 
particularly between 20° N. and 60° N. associated with extra-tropical weather. 


A mechanism is thus available for a large-scale but slow stirring up of 
the air between 8 and 25 km. bringing down the ozone from higher levels 
and destroying it in the lower levels by mixing it with water-vapour and dust 
(perhaps also methane) contained in the tropospheric air of lower latitudes. 


The problem of atmospheric ozone and its variations thus becomes a 
problem of the large-scale turbulence and circulation of the upper tropo- 
sphere and lower stratosphere. 


SUMMARY 


The paper contains a summary of the results of the height distribution of 
ozone in the atmosphere obtained at Mount. Abu with a Dobson Spectro- 
photometer during the period October 1951 to April 1952, when western 
disturbances affect North Indian weather. The observations show that there 
is very little ozone in the atmosphere below 18 km., that the ozone amount 
above 27 km. show very little systematic change with total amount, and 
that most of the changes take place in the layer 18 to 27km. These 
results are compared with those obtained at Delhi, Kodaikanal, New Mexico 
(with V, rockets), Arosa (in Switzerland) and Troms¢ (Norway). In middle 
and high latitudes, there are substantial amounts of ozone below 18 km., 
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and the larger the total amount, the larger is the proportion in the lower 
levels. For high ozone amounts in high latitudes, nearly half the total 
amount may be below 20 km. 


These results are discussed in relation to the variations of ozone with 
latitude and season, connecting them with the photo-chemical formation of 
ozone in the atmosphere above 27 km., the existence of a shelter-region for 
ozone between 27 km. and the tropopause and the destruction of ozone by 
mixing it with water-vapour and oxidisable matter carried up in the tropo- 
sphere by convective activity. 


It is concluded that the problem of variations of atmospheric ozone is 
one of large-scale turbulence and circulation in the upper troposphere and 
lower stratosphere. 
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